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ABSTRACT 

We present the results of a search for all embedded protostars with internal 
luminosities < 1.0 L in the full sample of nearby, low-mass star-forming regions 
surveyed by the Spitzer Space Telescope Legacy Project "From Molecular Cores 
to Planet Forming Disks" (c2d). The internal luminosity of a source, Li nt , is the 
£j [ luminosity of the central source and excludes luminosity arising from external 

heating. On average, the Spitzer c2d data are sensitive to embedded protostars 
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with Li n t > 4 x 10 3 (<i/140pc) 2 L , a factor of 25 better than the sensitivity of 
the Infrared Astronomical Satellite (IRAS) to such objects. We present a set of 



selection criteria used to identify candidates from the Spitzer data and examine 
complementary data to decide whether each candidate is truly an embedded pro- 
tostar. We find a tight correlation between the 70 /zm flux and internal luminosity 
of a protostar, an empirical result based on both observations and detailed two- 
dimensional radiative transfer models of protostars. We identify 50 embedded 
protostars with L int < 1.0 L Q ; 15 have L int < 0.1 L Q . The intrinsic distribu- 
tion of source luminosities increases to lower luminosities. While we find sources 
down to the above sensitivity limit, indicating that the distribution may extend 
to luminosities lower than probed by these observations, we are able to rule out a 
continued rise in the distribution below L int = 0.1 L . Between 75 — 85% of cores 
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classified as starless prior to being observed by Spitzer remain starless to our lu- 
minosity sensitivity; the remaining 15 — 25% harbor low-luminosity, embedded 
protostars. We compile complete Spectral Energy Distributions for all 50 ob- 
jects and calculate standard evolutionary signatures (Lb i, T^, and L^/ L smm ), 
and argue that these objects are inconsistent with the simplest picture of star 
formation wherein mass accretes from the core onto the protostar at a constant 
rate. 

Subject headings: stars: formation - stars: low-mass, brown dwarfs 

1. Introduction 

Recently, the Spitzer Space Telescope Legacy Project "From Molecular Cores to Planet 
Forming Disks" (c2d; Evans et al. 2003) completed an extensive 3.6 — 160 fxm imaging survey 
of nearby, low-mass star-forming regions. One of the results to come out of this survey is 
the discovery of very low luminosity objects (VeLLOs; Young et al. 2004). If the internal 
luminosity of a source, Li nt , is the total luminosity of the central protostar and circumstellar 
disk (if present), a VeLLO is defined to be an object embedded within a dense core with 
Lint < 0.1 L Q (Di Francesco et al. 2007). The bolometric luminosity of an embedded 
protostar, an observable quantity that can be calculated by integrating over the full Spectral 
Energy Distribution (SED), is composed of both internal and external luminosity {L^i = 
L int +L ext ). The external luminosity is usually that arising from heating of the circumstellar 
envelope by the Interstellar Radiation Field (ISRF), and will add, on average, a few tenths of 
a solar luminosity to L hol (e.g., Evans et al. 2001). Thus, the distinction between L bo i and L int 
is most relevant for embedded protostars with L int < 1.0 L , where the external luminosity 
can be a significant fraction of the observed Lboi- For VeLLOs, the external luminosity can 
dominate the observed Lw Radiative transfer modeling of the SEDs of embedded protostars, 
including both the emission from the envelope at submillimeter and millimeter wavelengths 
and the emission from the central source itself at infrared wavelengths, is required to decouple 
internal and external luminosities (e.g., Shirley et al. 2002; Young et al. 2004; Dunham et 
al. 2006). 

Several VeLLOs have been discovered in cores that were previously classified as starless 
prior to being observed by Spitzer. In fact, the very first starless core observed by c2d, 
L1014, was found to harbor a VeLLO with L int ~ 0.09 L (Young et al. 2004). This 
discovery reinforces that the known sample of embedded protostars is not complete. This 
sample has been assembled primarily by two methods: (1) searching for IRAS sources that 
are associated with dense cores and have colors consistent with those expected for embedded 
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protostars (e.g., Myers et al. 1987), and (2) identifying molecular outflows and radio point 
sources associated with dense cores indicating the presence of protostars too deeply embedded 
to detect with IRAS (e.g., Andre et al. 1993). Myers et al. (1987) found that the IRAS 
data could detect protostars with L int > 0.1 (c?/140pc) 2 L , where d is the distance to the 
protostar, although this does not include the younger, more deeply embedded protostars that 
were only identified on a case-by-case basis by the second method. The regions surveyed by 
c2d with Spitzer are located at distances ranging from 125 — 500 pc. Even in the closest of 
these regions VeLLOs are likely to fall below the IRAS sensitivity limit. In the more distant 
regions, no protostars with L int < 1 L would be detected. While some of these protostars 
might have been identified on a case-by-case basis as described above, the full sample of 
embedded protostars with L int < 1 L Q is clearly incomplete. 

Constructing a complete sample of embedded protostars with L int < 1 L Q is important 
for studies of low-mass star formation. Despite substantial progress in recent decades, the 
details of the physical processes regulating mass accretion from the envelope to the proto- 
star remain poorly understood (see McKee & Ostriker [2007] for a recent review). Several 
authors have attempted to constrain evolutionary models of the formation of low-mass stars 
by determining the observational signatures of these models and comparing them to the 
properties of known protostars (e.g., Myers et al. 1998; Young & Evans 2005). A result 
common to all such studies is a substantial population of protostars with luminosities below 
model predictions. An idea proposed to explain this discrepancy is that the mass accretion is 
episodic in nature and the protostars with the lowest luminosities are those observed in qui- 
escent accretion states (e.g., Kenyon & Hartmann 1995; Young & Evans 2005; Enoch 2007; 
Enoch et al. 2008a, in preparation). Theoretical studies have provided several mechanisms 
by which such a process may occur, such as material piling up in a circumstellar disk until 
gravitational instabilities drive angular momentum outward and mass inward in short-lived 
bursts (Vorobyov & Basu 2005, 2006). Alternatively, quasi-periodic magnetically driven out- 
flows in the envelope can cause mass accretion onto the protostar to occur in magnetically 
controlled bursts (Tassis & Mouschovias 2005). Indeed, the evidence for non-steady mass 
accretion in young protostellar systems still in the embedded phase is steadily growing (e.g., 
Hartmann & Kenyon 1985; Dunham et al. 2006; Acosta-Pulido et al. 2007; Kospal et al. 
2007; Fedele et al. 2007). However, as the sample of embedded, low luminosity protostars 
is incomplete, the true nature of the discrepancy between evolutionary models and obser- 
vations of protostars is unknown. Future work devoted to assessing the validity of various 
models by comparing their predictions to the properties of known protostars depends on the 
existence of a sample that is as complete and unbiased as possible. 

The VeLLOs are a particularly interesting subset of embedded, low-luminosity proto- 
stars; in essence, they are an extreme case of the problem discussed above. To date, only 
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three VeLLOs have been studied in detail: L1014-IRS (L int ~ 0.09 L ; Young et al. 2004), 
L1521F-IRS (L^ ~ 0.06 L ; Bourke et al. 2006), and IRAM 04191-IRS (L int ~ 0.08 L ; 
Dunham et al. 2006). Despite the fact that all three have similar internal luminosities, 
they differ greatly in envelope and outflow properties, as discussed by Bourke et al. (2006). 
IRAM 04191 drives a large, bipolar molecular outflow, features bright molecular line and 
dust continuum emission, and shows evidence for infall, depletion, and deuteration (Andre et 
al. 1999; Belloche et al. 2002). L1521F is also bright in molecular line and dust continuum 
emission and also shows evidence for infall, depletion, and deuteration (Crapsi et al. 2004), 
but the envelope is not as centrally condensed as IRAM 04191 and the presence of an outflow 
is uncertain (Crapsi et al. 2004; Bourke et al. 2006). L1014 does not show significant evi- 
dence for infall, depletion, or deuteration (Crapsi et al. 2005a), but it does drive a compact, 
weak molecular outflow detected only in interferometer observations (Bourke et al. 2005; 
Crapsi et al. 2005a). A systematic search for all VeLLOs in the regions surveyed by c2d 
will allow their properties to be examined in detail both on a case-by-case basis and as a 
class of objects. Identifying the complete sample of VeLLOs will also allow us to determine 
how many cores classified as starless prior to being observed by Spitzer truly are starless, a 
question with important implications for estimates of the lifetime of starless cores (e.g., Kirk 
et al. 2005). 

In this paper, we present the results of a search for all embedded protostars with L int < 
1.0 L in the full c2d imaging dataset. Depending on the distance to each individual source, 
some will already have been detected by IRAS, while others will be new sources. We consider 
this work to be complementary to several related studies: A search by Kirk et al. (2007) for 
embedded protostars in 22 cores classified as starless prior to being observed by Spitzer; a 
search by J0rgensen et al. (2007; 2008, in preparation) for all embedded protostars in Perseus 
and Ophiuchus, regardless of luminosity, conducted by combining Spitzer and SCUBA 850 
fim dust continuum emission data; and a search by Enoch (2007) and Enoch et al. (2008a, 
in preparation) for all embedded objects in the Perseus, Ophiuchus, and Serpens molecular 
clouds, regardless of luminosity, conducted by combining Spitzer and Bolocam 1.1 mm dust 
continuum emission data. 

The key difference between the work presented here and the searches for embedded 
protostars listed above is that we do not start by identifying dense cores from their millimeter 
dust continuum emission and then look for associated Spitzer sources embedded within 
them. Instead, we develop a set of criteria to identify candidate embedded, low-luminosity 
protostars based on the 3.6 — 70 /im Spitzer data. This way, we are able to identify all of 
the candidates in the full c2d dataset, regardless of the availability and quality of millimeter 
wavelength observations for each region. Only after we identify all candidates based on 
Spitzer data alone do we turn to other observations to distinguish the objects of interest 
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from various contaminants masquerading in our sample. Our method identifies candidates 
for further examination once large-scale surveys of nearby star-forming regions are completed 
with SCUBA-2 (Ward-Thompson et al. 2007) and Herschel, and the method can easily be 
extended to search for embedded, low-luminosity protostars in the additional nearby, low- 
mass star-forming regions being surveyed by the Spitzer Gould Belt Legacy Project (L. Allen 
et al. 2008, in preparation). 

The organization of this paper is as follows: In §2], we provide a brief description of the 
c2d observations and data reduction, emphasizing those aspects relevant to this work. The 
criteria for identifying candidate embedded, low-luminosity protostars from the 3.6 — 70 /im 
Spitzer data are discussed in §|3l along with the possibilities for estimating source internal 
luminosities directly from observable quantities. A general proof-of-concept demonstrating 
the validity of these criteria is given in §H In £0, we discuss the contamination expected in 
the list of candidates, both from background extra-galactic sources and from more evolved 
Young Stellar Objects (YSOs) no longer embedded within their dense cores. We discuss the 
necessary requirements to prove that a candidate is truly an embedded protostar in §5.14 
we apply these requirements to our candidate list in §5.21 and we discuss the difficulties in 
including regions lacking good quality 70 /im data in §5.31 We discuss several general results 
of this work in Finally, we present our conclusions in £J7J 

2. Observations and Data Reduction 

The Spitzer c2d data have been published in several papers focusing on individual 
regions in the survey. We summarize here the observation strategy and data reduction 
method used by c2d; greater detail on individual regions can be found in the references 
given below. 

2.1. Observations 

2.1.1. Molecular Clouds 

The c2d project obtained complete 3.6 — 160 /im Spitzer maps of five nearby, large 
molecular clouds. The clouds (and approximate areas) surveyed include 1 deg. 2 of Serpens 
(Harvey et al. 2006; 2007a; 2007b), 2.5 deg. 2 total of the Lupus I, Lupus III, and Lupus 



- 6- 



IV cloudsE (Chapman et al. 2007; Merin et al. 2008), 1 deg. 2 of Chamaeleon II (Young 
et al. 2005; Porras et al. 2006; Alcala et al. 2008), 7 deg. 2 of Ophiuchus (Padgett et al. 
2007), and 4 deg. 2 of Perseus (J0rgensen et al. 2006; Rebull et al. 2007). These clouds were 
chosen to span a wide range of low-mass star- forming environments (Evans et al. 2003). We 
adopt distances of 260 ± 10, 150 ± 20, 200 ± 20, 150 ± 20, 178 ± 18, 125 ± 25, and 250 ± 50 
pc for Serpens, Lupus I, Lupus III, Lupus IV, Chamaeleon II, Ophiuchus, and Perseus, 
respectively (see discussion of distances in the references listed above). The observations 
of Serpens, Lupus, Chamaeleon II, Ophiuchus, and Perseus were obtained in Program IDs 
(PIDs) 174, 175, 176, 177, and 178, respectively. 

Observations at 3.6, 4.5, 5.8, and 8.0 /zm were obtained with the Infrared Array Cam- 
era (IRAC; Fazio et al. 2004). Each region was observed in two epochs to allow for the 
detection and removal of asteroids, and each region was observed twice per epoch with small 
dithers between the observations. The integration time per pointing was selected to be 12 s, 
providing an effective integration time per pixel of 10.4 s. This strategy results in four total 
observations per region (2 epochs of 2 dithers each), resulting in a total integration time 
of ~ 42 s per pixel. All observations were obtained in the "High Dynamic Range" mode, 
which provides an additional short (0.6 s) exposure to enable photometry of bright sources 
saturated in the longer exposures. 

Observations at 24, 70, and 160 fim were obtained with the Multiband Imaging Pho- 
tometer for Spitzer (MIPS; Rieke et al. 2004). The observations consist of two epochs of 
fast scan maps with a spacing between adjacent scan legs of 240". The second epoch was 
offset by 125" from the first in the cross-scan direction to fill in the 70 fim sky coverage 
that would have otherwise been missed due to detector problems, and by 80" in the scan 
direction to minimize missing 160 /im data. These mapping parameters result in two epochs 
of observations at 24 /im and one epoch at 70 and 160 fim, with total integration times per 
pixel of 30, 15, and 3 s, respectively. Small gaps of missing coverage remain in the 160 /im 
maps. 

For technical reasons relating to the fact that the length of a MIPS scan map leg must 
be chosen from a fixed list, the regions observed by MIPS are often significantly larger than 
those observed by IRAC in order to ensure full coverage. The approximate areas quoted 
above are for the IRAC+MIPS overlap regions; in general, the areas observed by MIPS are 
a factor of 2 — 3 larger. 



1 For simplicity, we refer to the ensemble of these three portions of the Lupus Molecular Cloud Complex 
as simply "Lupus" . 
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2.1.2. Small Dense Cores 

In addition to the five large molecular clouds, c2d also obtained small maps of 82 regions 
containing 95 small, dense cores. These cores are selected on the basis of being nearby (within 
500 pc), relatively small (generally less than 5', which corresponds to ~0.7 pc at a distance 
of 500 pc), relatively isolated, and showing evidence for dense gas and dust (Evans et al. 
2003). The list of all 82 regions is given in Evans et al. (2007). A discussion of the overall 
results of the Spitzer observations of these cores and estimates of the distance to each core 
are presented in T. Huard et al. (2008, in preparation). Although "relative isolation" was 
one of the criteria used to select these cores, we note here that many are not truly isolated 
in the sense that they are often loosely associated (both in projection on the sky and in 
velocity) with larger complexes. Approximately 70% of these cores were classified as starless 
prior to the launch of Spitzer, based primarily on showing no association with IRAS sources 
(Evans et al. 2003 and references therein). 

Observations at 3.6, 4.5, 5.8, and 8.0 fim were obtained with IRAC in PID 139 in a 
manner nearly identical to the clouds. The main differences are that the maps are much 
smaller (often only a single 5x5' frame or a small map < 30' on a side), the number of 
epochs observed (one or two) depends on the location of each core relative to the ecliptic, and 
short-exposure HDR images were only obtained for cores expected to contain bright sources 
based on previous surveys. Cores observed in only one epoch featured 4 dithers instead of 2 
so that all cores featured a total integration time per pixel identical to the clouds. 

Observations at 24 and 70 /im were obtained with MIPS in PID 139; 160 fxm observa- 
tions were not obtained. Unlike the clouds, the MIPS observations of the core regions were 
obtained in the pointed photometry mode. The exposure time per pointing was 3 seconds at 
24 jum and either 3 or 10 seconds at 70 /im, and small raster maps were obtained to match 
the IRAC maps. As a result, the core regions, unlike the clouds, do not feature significantly 
larger areas observed with MIPS than with IRAC. The exact integration time per pixel varies 
from region to region due to the differences in field and map size, but is generally 30 — 60 s 
at 24 /im and 50 — 100 s at 70 /im. 

2.2. Data Reduction 

The IRAC and MIPS images were processed by the Spitzer Science Center (SSC), using 
their standard pipeline, version S13, to produce Basic Calibrated Data (BCD) images. These 
images were then improved by the c2d Legacy Project to correct for artifacts remaining in 
the BCD images. A complete description of the improvements made, as well as the source 
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extraction and band-merging process summarized below, is given in Evans et al. (2007). 

After correcting for artifacts, mosaics were produced using the MOPEX software pro- 
vided by the SSC. Photometry at 3.6 — 24 /im was obtained using c2dphot, a modified 
version of DOPHOT (Schechter et al. 1993) that utilizes a digitized rather than analytic 
point source profile to better match the real Spitzer data, as well as incorporating several 
other changes (a complete description of which are given in Harvey et al. 2006). Sources 
that were not detected in all 5 bands from 3.6 — 24 /im were band-filled, a process whereby 
the c2dphot source extractor was forced to obtain fluxes at the source position (known from 
the band(s) in which it was detected) in the bands for which the source was not detected in 
the original source extraction (see Evans et al. (2007) for details). Finally, photometry at 
3.6 — 24 /im were band-merged into a final source catalog. 

Photometry at 70 /im was obtained using the SSC's MOPEX point-source fitting package 
on filtered BCD data for faint sources (< 2 Jy) and unfiltered BCD data for bright sources 
(> 2 Jy). Sources extracted at 70 were band-merged into the 3.6 — 24 /im catalogs 
described above, although this process is imperfect due to the significantly worse resolution 
at 70 /im (~18") compared to 24 /im (~6") and 3.6 - 8.0 /im (~2") (see Evans et al. 2007 
and ED - 

All of the above reduction and source extraction is part of the standard c2d pipeline. 
However, this pipeline does not include source extraction at 160 /im due to problems with 
diffuse emission and low spatial resolution. Thus, photometry at 160 /im for the five large 
clouds was obtained using the SSC's MOPEX point-source fitting package. The source 
extraction lists were then inspected and sources that appeared to be extended dust "clumps" 
with little or no central condensation were removed (T. Brooke, private communication). As 
a redundancy check, we also visually inspected all of the 160 /im images at the positions of 
the candidates identified in this work and performed aperture photometry on those that have 
associated 160 /im sources. The aperture sizes and sky annuli were chosen to allow us to use 
one of the standard aperture corrections for 160 /im point sources determined by the SSC. 
In general, we find very good agreement between the PSF and aperture photometry. The 
ratio of PSF to aperture flux for sources extracted by both methods has a mean and median 
of 1.07 and 1.02, respectively, with a standard deviation of 0.2. We thus add an overall 20% 
systematic uncertainty into our photometry at 160 /im. Both methods assume all 160 /im 
sources are point sources; fluxes of extended sources may thus be under-estimated. 
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3. Identification of Candidate Embedded, Low-Luminosity Protostars 

Generally speaking, we identify candidate embedded, low-luminosity protostars from the 
Spitzer observations by selecting sources with rising fluxes from shorter to longer wavelengths 
and infrared luminosities Lir, calculated over the 2MASS and Spitzer bands (1.25 — 70 fim), 
indicating L int < 1 L (see below). Throughout this paper we denote flux as F u (F v = uS u ). 
When we refer to the shape of the Spectral Energy Distribution (SED) of an object, it is 
always evaluated in terms of flux (F u ) rather than flux density (S u ). 

At this stage, we emphasize completeness over reliability; many objects identified as 
candidates are likely either more evolved Young Stellar Objects (YSOs) no longer embedded 
in their dense cores, or background galaxies. Separating these from the true embedded 
protostars is addressed in £j5l We have assembled a set of seven criteria that must be 
met for a source to be identified as a candidate. We list these criteria below, followed by 
an explanation of each one individually in §3.11 We provide a detailed discussion of the 
motivation for the third criterion in §3.21 and investigate the results of relaxing the first 
criterion in §3.31 Finally, we discuss the need to visually inspect each candidate in §3.41 and 
present a general proof-of-concept in §HJ 

The seven criteria that must be met are: 

1. Detected at 24 jttm with a signal-to-noise ratio (SNR) > 3. 

2. Located at a position observed by a minimum of two Spitzer photometric bands be- 
tween 3.6 and 70 fim. 

3. Detected at 70 jttm, unless confused or not observed at this wavelength. 

4. Rising (or flat) SED between the longest IRAC wavelength at which the source is 
detected and 24 fim. 

5. Rising (or flat) SED between 24 and 70 fim, unless confused or not observed at 70 /im. 

6. Observed luminosity calculated using all detections ranging from 1.25 — 70 fim of Ljr 
< 0.5 L . 

7. Not classified as a candidate galaxy in the classification method of Harvey et al. (2007b) 
unless Log(Pg a i) < —1.25, where P ga i is the un-normalized probability of being a galaxy. 
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3.1. Selection Criteria 

We impose the first criterion simply because 24 /mi observations are convenient for 
determining whether or not a source features a rising SED in the wavelength region covered 
by our Spitzer observations (3.6 — 70 /mi). First of all, 24 /mi is located approximately 
halfway between 3.6 and 70 /mi in log space, making it useful for determining the general 
shape of the SED over the wavelengths observed by Spitzer. Additionally, the diffraction 
limited resolution of Spitzer, 9 = 1.2-^, is a factor of 2.92 better at 24 /mi than 70 /mi. Thus, 
sources clustered close together (at least in projection on the sky) can often be resolved into 
individual sources at 24 /mi even when they are confused at 70 /mi, providing reliable flux 
information at a wavelength beyond the 3.6 — 8.0 /mi region covered by IRAC and better 
enabling us to evaluate the shape of the SED. We discuss the effects of relaxing this criterion 
in gH 

The second criterion ensures that we have flux information available at a minimum 
of one other Spitzer wavelength besides 24 /mi. The MIPS observations of the five large 
clouds were obtained in the scan map observing mode and often cover significantly larger 
areas than the IRAC observations of the same regions for technical reasons. Furthermore, 
there are minor differences between the exact areas observed at 24 and 70 /im, resulting 
in small regions that are observed at 24 /mi but not at any other wavelength. Since it is 
impossible to determine the shape of the SED of any source observed at only one wavelength, 
all sources in these "MIPS-24 only" regions are removed from consideration. This issue is 
not as significant in the regions with dense cores, where the MIPS data were obtained in the 
pointed observations mode and often cover very similar areas to the IRAC observations. 

The third criterion is imposed because 70 /im is a crucial wavelength for determining 
Li n t of embedded, low-luminosity protostars. We will return to this point in greater detail 
below. 

The fourth and fifth criteria select sources with rising or flat SEDs from shorter to 
longer wavelengths. Simple, 1-D models of protostars still embedded within their cores show 
that they should essentially always feature rising SEDs both between the IRAC wavelength 
regime (3.6 — 8.0 /mi) and 24 /im and between 24 and 70 /ma (e.g., Young et al. 2005). Two- 
dimensional models that include the effects of outflow cavities show that the exact shape of 
the 3.6 — 70 /mi SED depends on the inclination of the source (e.g., Whitney et al. 2003a; 
Whitney et al. 2003b; Robitaille et al. 2006). In particular, embedded protostars observed 
at nearly pole-on inclinations exhibit much flatter SEDs over these wavelengths since more of 
the protostellar emission can escape through the outflow cavities instead of being reprocessed 
to longer-wavelength emission by the envelope. While this flattening of the SED can be very 
significant, even at the most extreme inclinations the models predict essentially flat rather 
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than falling SEDs at these wavelengths. We return to this point in §4.21 

We thus require each source to have an increasing or constant flux between the longest 
IRAC wavelength at which the source is detected and 24 /mi, and between 24 and 70 /im. If 
a source is not detected at any IRAC wavelength, it is considered to pass the fourth criterion. 
We take the photometric uncertainties into account; a source may actually have a decreasing 
flux between the two wavelengths considered as long as it is consistent with being at least 
flat given the uncertainties. Not all sources have 70 /im detections; we will return to these 
sources at the end of this section, but for now, we keep any source that meets the first rising 
criterion but has no flux information at 70 /mi. We do not place any rising requirements on 
the IRAC fluxes themselves since the effects of geometry (Whitney et al. 2003a; Whitney et 
al. 2003b), solid-state ice features (Boogert et al. 2004), and shocked emission from outflows 
(e.g., Noriega-Crespo et al. 2004) greatly complicate the exact shape of the 3.6 — 8.0 /im 
SED. 

The sixth criterion separates out the low-luminosity sources, those of interest to this 
study, from the more luminous sources. As described in £[TJ the internal luminosity of an 
embedded, low-luminosity protostar is not an observable quantity and is usually derived from 
radiative transfer modeling. Here we investigate whether or not we can estimate Li nt without 
the need to construct detailed models of each source. We have identified 11 embedded 
protostars from the literature that have been observed with Spitzer by c2d and have accurate 
internal luminosity determinations from radiative transfer models. We list these 11 objects 
in Table (TJ along with Lj R calculated from the observations and Lf, i and L int derived from 
the published models. 

Figured^ plots L int vs. Lm in log-log space for these 11 objects; a linear correlation is 
clearly observed. The solid line shows the results of a linear least-squares fit; it has a slope 
of 0.88 and a y-intercept of 0.32. However, any linear correlation in a luminosity-luminosity 
plot is suspect since it may simply arise from plotting the square of the distance to the object 
versus itself. Thus, Figure [lb plots the ratio of L int /LiR, a quantity that is independent of 
distance, vs. the log of Ljr. An approximately constant ratio is observed over several orders 
of magnitude of Lj R , as would be expected if a true linear correlation exists between these 
quantities. The solid line shows the average ratio of 1.7, weighted by the uncertainties in L int . 
Inverting this quantity, we find that the luminosity calculated between 1.25 — 70 /im gives 
a result that is approximately half the true value of the internal luminosity. Thus, to select 
sources with L int < 1.0 L , we impose the requirement that L IR < 0.5 L . We emphasize 
that this relationship between L int and L IR is only an approximation; Figure [Tb clearly 
shows there is some variation from source to source. In fact, the lowest-luminosity sources 
may have a larger ratio of Li nt /Lm than more luminous sources, although the uncertainties 
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are large. However, even if the ratio of Li nt /LiR does increase for lower-luminosity sources, 
they will have Ljr far below the cut of 0.5 L and will thus be included in our sample. 

The seventh criterion is related to the fact that background galaxies are difficult to 
separate from YSOs when analyzing Spitzer data of star-forming regions. A number of 
authors have recently shown that such galaxies often exhibit similar mid-infrared colors and 
SEDs to YSOs (e.g., Harvey et al. 2006; J0rgensen et al. 2006). In an effort to accurately 
separate out background galaxies from YSOs, Harvey et al. (2007b) developed a classification 
method that assigns to each object that can not be fitted by an extincted stellar photosphere 
an un-normalized "probability" of being a galaxy, P ga i- The probabilities are assigned based 
on the position of each source in various color-color and color-magnitude diagrams, along 
with information about the nature of the source (point-like vs. extended) given by the c2d 
source extraction software (see $2] and references therein). Based on its value of P ga i, Harvey 
et al. classified each source as a candidate YSO or candidate Galaxy. 

Figure [2] shows a histogram of P ga \ for all 851 sources assigned a probability in the 
ensemble of 82 regions with dense cores observed by c2d. The dotted line shows Log(P ga i) = 
—1.47 while the dashed line shows Log(P ga i) = —1.25. This figure is similar to Figure 4 
of Harvey et al. (2007b), which plots the distribution of P ga i for all sources assigned a 
probability in the c2d observations of the Serpens Molecular Cloud. As noted by Harvey et 
al., there are two distinct peaks in probability space: one at Log{P ga {) = —5.00 and one at 
Log(P ga i) ~ —0.4, with a tail connecting the two. In a similar analysis on a Spitzer dataset 
of extragalactic objects, Harvey et al. found only the peak at Log(P ga i) ~ —0.4 and not the 
other peak or the tail connecting the two. Based on this, they set the boundary between 
candidate YSO and candidate galaxy at Log(P ga i) = —1.47, such that objects in the peak 
at Log(P ga i) = —5.00 and the tail between the two peaks are considered candidate YSOs 
while objects in the peak at Log(P ga i) ~ —0.4 are considered candidate Galaxies. However, 
as is evident both from Figure [2] of this paper and Figure 4 of Harvey et. al, there is likely 
some overlap in probability space between the tail of objects considered candidate YSOs and 
the peak of objects considered candidate galaxies. Some of the objects near the boundary 
of Log(P ga i) = —1.47, but slightly into the regime of galaxies, may in fact still be YSOs. 
Furthermore, there was no galaxy in the Spitzer dataset of extragalactic objects considered 
by Harvey et al. with Log(P ga i) < —1.25. Thus, for all the objects that pass the first 
six criteria, we reject all those classified as candidate galaxies based on this classification 
method except for those with —1.47 < Log{P ga \) < —1.25. Many of these intermediate 
objects are likely to be galaxies, but we include them as candidates at this stage for the sake 
of completeness. Separating out these false candidates is the subject of §0 
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3.2. Importance of 70 jim Data 



We now return to the issue of detection at 70 /xm. This is a crucial wavelength for 
determining Lj nt for embedded, low-luminosity protostars. Radiative transfer models of 
these objects are strongly constrained by this wavelength since the flux at 70 /zm, F 70 , 
is mostly unaffected by the details of the source geometry and presence or absence of a 
circumstellar disk that significantly affect the 3.6 — 24 /zm SED. F-?o is also quite unaffected 
by the amount of external heating from the ISRF that is important at wavelengths > 100 
/iin (e.g., Dunham et al. 2006). To examine this more fully, we again use the sources listed in 
Tabled] observed with Spitzer by c2d and possessing accurate determinations of L int through 
radiative transfer models. Figure [3] plots F u , normalized to 140 pc, for all six wavelengths 
observed by Spitzer vs. L int for these 11 sources, along with the results of linear least-squares 
fits to the data in log-log space. These fits show that a correlation exists between the flux at 
each wavelength and Li nt \ in general terms, more luminous sources emit more energy at all 
wavelengths. However, the correlation is clearly strongest at 70 fim, weaker at 24 /zm, and 
very poor at 3.6 — 8.0 /zm. To quantify this, at each wavelength we calculate x% a reduced 
chi-squared for the fit: 



where n is the number of datapoints (n = 11), p is the number of free parameters 
(p = 2), yi and Xi are the log of the x and y values for each datapoint, a. t is the size of the 
uncertainty in log space, and m and b are the slope and y-intercept, respectively, derived 
from the fit in log-log space. The values for m and b, along with their statistical uncertainties 
from the fit and the results of the xl calculations, are presented in Table [2j We find xl — 3 
at 70 jum, Xr — 85 at 24 /im, and xl > 100 for all four IRAC wavelengths, in agreement with 
the qualitative conclusions from above. 

These results suggest that a direct correlation exists between F 70 and L int . However, 
this analysis is based on only 11 data points, and furthermore, many of the radiative transfer 
model determinations of L int listed in Table [1] are based on the same simple, 1-D modeling 
setup that, among other things, uses a very simple method of including a disk in a 1-D 
model (Butner et al. 1994). We do not expect this to be significant, as previous work has 
concluded that even these simple, 1-D models accurately constrain L int (e.g., Dunham et al. 
2006), but to test the validity of these results, we attempt to reproduce them by running a 
grid of two-dimensional models that are more physically realistic than those on which many 
of the values of L int in Table [T] are based. 




(1) 
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We used the two-dimensional, axisymmetric, Monte Carlo dust radiative transfer code 
RADMC (Dullemond & Dominik 2004) to compute the radiative transfer and calculate the 
emergent SEDs of embedded, low-luminosity protostars and look for correlations between 
L int and F v in the Spitzer bands. The density structure assumed for the models consists of a 
protostellar disk and a rotationally flattened, infalling envelope with an outflow cavity, and 
is described in greater detail in Crapsi et al. (2008). The details of the models are nearly 
identical to those of Crapsi et al.; we do not describe them here except to note that they 
have been tuned to their specific purpose within this work in the following ways: 

• The protostar has a fixed temperature of 3000 K and a random luminosity in the range 
of 0.03 - 10 L . 

• The disk has a fixed size of 100 AU, variable mass between 10~ 3 and 10~ 5 M Q , and fixed 
flaring parameters (a scale height H of 20 AU at the outer radius and then decreasing 
as H oc r 9 / 7 , corresponding to the self-irradiated passive disk of Chiang & Goldreich 
[1997]). 

• The envelope has a fixed radius of 14,000 AU, a variable mass in the range of 1 — 10 
M , and a variable centrifugal radius in the range of 100 — 900 AU. 

• The envelope is externally heated by the ISRF, for which we adopt that of Black (1994), 
modified at the ultraviolet wavelengths to reproduce the ISRF of Draine (1978). This 
Black-Draine ISRF is then attenuated by dust with properties given by Draine & Lee 
(1984) to simulate being embedded in a parent cloud and multiplied by a scale factor 
to account for the fact that the ISRF is not uniform everywhere, where the number of 
magnitudes of visual extinction and scale factor are chosen randomly. 

The observed SED was obtained by raytracing the density, temperature, and scattering 
structures calculated by RADMC at five different inclinations: 20, 35, 50, 65, and 80° from 
the axis of symmetry. Fluxes were integrated inside apertures comparable to the resolution 
of Spitzer: 2" for wavelengths shorter than 10 /im, 6" for wavelengths in the range of 10 — 40 
/im, and 20" in the wavelength range of 40 — 100 /im. 

Following this description, a grid of 292 models were run and 1460 SEDs were obtained 
(one at each inclination). The results are shown in Figure H] and, qualitatively, appear to 
confirm the correlation between Fjq and Lint- The parameters derived from linear least- 
squares fits in log-log space for F u vs. L int for the models are given in Table [2J It is clear 
from these models that there is in fact a tight correlation between F70 and L int (most of the 
scatter seen in the MIPS2 panel of Figure H] is due to inclination effects). In fact, not only do 
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the observations and models both show this correlation, but they also agree remarkably well 
on the quantitative relationship (within 3% for the slope and within 1% for the y-intercept) 
considering that, to date, only 11 sources have been modeled by us and are thus available 
to examine this from an observational standpoint. Since F 7Q is an observable quantity, the 
correlation we find between F m and L int can be used to estimate L int for sources which either 
lack sufficient data to constrain radiative transfer models or have such data but have not yet 
been modeled. Using the results of the linear least-square fit to the modeled F 70 vs. L int in 
log-log space, we derive this estimate to be 

L int =3.3x lO 8 F 7O a94 L , (2) 

where F 7 o is in cgs units (erg cm -2 s" 1 ) and is normalized to 140 pc. 

Since a detection at 70 fim is crucial for constraining radiative transfer model determi- 
nations of L int and gives an immediate, direct estimate of L int , we require each source to 
be detected at 70 /zm unless it is not observed at this wavelength or is located too close to 
another source detected at 70 //m to be resolved into a separate source (see §5.3|) . Thus, with 
these two exceptions, we effectively require candidates to be above the detection threshold 
at 70 jtim. This requirement sets the fundamental limit for our luminosity sensitivity. Fig- 
ure El plots a histogram showing the 70 /im 3a point source sensitivity for each of the 82 
regions with dense cores observed by c2d. This sensitivity was derived by first calculating 
the standard deviation of the background intensity for each region, a s k y , in the image units 
of MJy sr _1 . Approximating the 70 /im Spitzer PSF as a Gaussian with #fwhm equal to 
the diffraction limited resolution of ~ 18", this background intensity is then converted into 
units of mJy beam _1 as follows: 

a sky (mJybeam)- 1 = 1 x 10 9 ( ^f™ ) ^s ky (MJy sr)' 1 (3) 

This gives the 3a point source sensitivity for each region, since the total flux density of a 
point source is equivalent to its flux density per beam. 

This 70 /j,m point source sensitivity can be directly translated into a luminosity sensi- 
tivity using Equation [21 From this relationship, the mean 70 /im 3a point source sensitivity 
of 38.6 mJy translates into a luminosity sensitivity of ~ 4 x 10~ 3 L Q at 140 pc. Figure [6] 
plots the luminosity sensitivity limit for each of the 82 regions vs. the distance to the region, 
calculated by translating the 70 /xm 3a point source sensitivity into a luminosity sensitivity 
using Equation [2] and then scaling from 140 pc to the distance to the region. The solid 
line shows the relation Li nt = 4 x 10~ 3 (d/140pc) 2 L . There is some scatter between the 
line and the actual data, caused by the variation in 70 jum point-source sensitivity from one 
region to the next, but this relation clearly provides an accurate estimate of the luminosity 
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sensitivity. On average, the c2d observations of regions with dense cores are sensitive to 
embedded protostars with internal luminosities L int > 4 x 1CT 3 (<i/140pc) 2 L . This same 
argument holds true for the c2d observations of the clouds, as the 70 /zm cloud maps reach 
similar depths as the cores. This is a factor of 25 better than the sensitivity of IRAS quoted 
in 

Applying these criteria results in the identification of 673 candidate embedded, low- 
luminosity protostars: 106 in the ensemble of dense core regions, 196 in Perseus, 112 in 
Chamaeleon II, 153 in Lupus, 57 in Ophiuchus, and 49 in Serpens. 

3.3. Relaxing the First Criterion 

The first of the seven criteria described above for selecting candidates from the Spitzer 
c2d observations is a detection at 24 /xm. As discussed in §3.11 this is an important wavelength 
for determining whether or not a source features a rising (or flat) SED from 3.6 — 70 /xm. 
However, it is possible that a source of interest has no detection at 24 /im. For example, 
it might simply not be covered at 24 fim if it is located close to the edge of an observed 
region where the coverage between bands is not uniform, or it might be so deeply embedded 
that it is only detected at 70 /xm. The latter could be particularly interesting sources, either 
very young, deeply embedded protostars or perhaps even first hydrostatic cores, which are 
short-lived, hydrostatic objects in a phase between the initial collapse of a dense core and 
the formation of a Class protostar predicted by theory but not yet found by observations 
(e.g., Boss & Yorke 1995). 

Thus, to make sure no sources are missed by imposing this first criteron, we search for 
any source not detected at 24 fim but detected at 70 /xm. We then apply the sixth criterion, 
Lir < 0.5 L , to select the low luminosity sources that are the focus of this work. This 
gives an additional 53 candidates: 13 in the ensemble of dense core regions, 7 in Perseus, 
17 in Chamaeleon II, 10 in Lupus, 2 in Ophiuchus, and 4 in Serpens. We add these to 
the candidates identified above for a total of 726 candidates. However, only 4 of these 
53 additional candidates survive the visual inspection described in §3.4[ and none of these 
surviving 4 are particularly strong candidates (see §3.41) . 

3.4. Visual Inspection 

The final step in the identification of candidates is visual inspection. For each of the 
726 sources identified above, we examined the images at the six wavelengths observed by 
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Spitzer and the source SED. Throughout the course of this inspection, we found five reasons 
to reject candidates: (1) The source is obviously a resolved galaxy; (2) The 70 /im detection 
is not real; (3) The source SED is not consistent with that of an embedded protostar; (4) 
The source is located in a region of nebulosity that produces a false infrared excess at the 
longer wavelengths; and (5) There is no real source. 

The first reason is due to the fact that a source must meet certain requirements to be 
classified as a candidate YSO or candidate galaxy in the method of Harvey et al. (2007b). 
Specifically a source must be detected in both epochs of observations (or in the single epoch 
if only one epoch was observed) with SNR > 3 in all 4 IRAC bands and the 24 /im MIPS 
band. Any source outside the area of overlap between all 5 bands, or simply not detected 
in any one of these bands, does not get evaluated by the method of Harvey et al. and thus 
can not be classified as a candidate YSO or galaxy. This results in many galaxies making 
the candidate list. Visual inspection removes those large and/or close enough to be resolved 
by IRAC; all others remain candidates at this stage. 

The second reason relates to the fact that the exact sensitivity of the 70 /im c2d ob- 
servations is a complicated function of the position within the scan map (clouds) or pointed 
observation (core regions), since the exact number of frames at a given position depends 
on the technical details of how Spitzer executes observations. In regions with fewer frames, 
often near the edges of maps, noise spikes are sometimes falsely identified as sources; we 
remove these from consideration. Of the 53 candidates identified by searching for sources 
detected at 70 /im but not at 24 /im, 49 were removed for this reason, and the 4 that remain 
are all sources with questionable 70 /im detections that simply weren't obvious enough false 
detections to remove. 

The third reason relates to the details of how rising sources were selected. Figure [7] 
shows the SED of SSTc2d J032856. 64+311835. 6, a source in Perseus representative of those 
removed for this reason. It is detected at all 6 Spitzer wavelengths, features a rising flux 
between 8 and 24 /im, features a flat, but slightly rising, flux between 24 and 70 /im, has 
Lir = 0.37 L Q , and is classified as a candidate YSO by the classification method of Harvey 
et al. (2007b). Thus, it fulfills all the criteria for selecting candidates discussed in §3.11 
However, it clearly emits more energy in the near-infrared (1.25 — 2.17 /im) than at 24 and 
70 /xm. Comparison with 2-D models of sources (e.g., Whitney et al. 2003a; Whitney et al. 
2003b; Robitaille et al. 2006) shows that, even at the most extreme pole-on inclinations, the 
effects of outflow cavities are unlikely to result in an observed SED more luminous in the 
near-infrared than at 24 and 70 /xm. A more likely explanation for objects featuring such 
SEDs is that they are more evolved YSOs surrounded by circumstellar disks but no longer 
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embedded within dense cores. Indeed, a search of the SIMBAD^ database shows that this 
is SSS 108, which Aspin et al. (1994) and Aspin (2003) conclude is a pre-main sequence 
star with a small thermal excess at 2 /im (i.e., from a disk). The nature of these objects 
with large infrared excesses at 24 and 70 /im compared to 3.6 — 8.0 fim could potentially 
be very interesting; indeed, these SEDs appear similar to those of a sample of YSOs with 
disks featuring large inner holes presented by Brown et al. (2007). Sources that meet our 
selection criteria but are more luminous in the near-infrared than at 24 and 70 /im may be 
interesting objects, but they are not relevant to this study and are thus removed from our 
sample. 

The final two reasons for rejecting sources upon visual inspection both relate to the 
bandfilling process described in §2J The first reason is that stars are often detected only in 
the first two IRAC bands, due both to the decreasing Spitzer sensitivity at longer wavelengths 
and the fact that Spitzer observes the Rayleigh- Jeans portion of stellar SEDs. Some stars 
detected only in these two bands are located, at least in projection, in regions of nebulosity. 
For these objects, the bandfilling process can assign source fluxes in bands beyond IRAC2 
that are contaminated by this extended nebulosity, creating false excesses and thus falsely 
rising SEDs. These are obvious in that there is no real point-source in the longer wavelength 
images. The second reason is that Spitzer images of star-forming regions often exhibit copious 
amounts of nebulous, diffuse emission due to scattered light, shocked emisison from outflows, 
thermal emission from hot dust, PAH emission, etc. In regions where this diffuse emission 
becomes clumpy, the source extraction may falsely detect a point-source at one wavelength 
that then gets band-filled at all other wavelengths. Some of these fake sources may pass all 
7 criteria for selection, but they are obvious upon visual inspection in that there is no point 
source present in the images. Both types of sources are removed from consideration. 

After removing sources for the reasons described above, we are left with 218 candidates: 
49 in the ensemble of dense core regions, 70 in Perseus, 16 in Chamaeleon II, 42 in Lupus, 
22 in Ophiuchus, and 19 in Serpens. The source number from this work, group number 
indicating likelihood of being an embedded object (see §5]), Spitzer source name, position, 
L IR , and distance to and name of the nearest source in the IRAS point-source catalog are 
given for each candidate in Tables [3] — [HI For the candidates in the dense cores, Table [3] 
also list the name of the core region in which each candidate is located and the assumed 
distance to this region; these distances have been derived by reddening or association with 
regions with known distances (T. Bourke, private communication). The full list of distances 
to core regions will be presented by T. Huard et al. (2008, in preparation). We emphasize 



2 Available at: 



http : / /simbad . u-strasbg . fr /simbad / 
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here that these objects are only candidates. We have chosen to emphasize completeness over 
reliability, with the result being that many of these 218 candidates are in fact galaxies or 
more evolved YSOs no longer embedded in their dense cores. Separating the wheat from the 
chaff is the focus of §0 



4. Proof of Concept 

In §3.1[ we presented a set of criteria to select embedded, low-luminosity protostars from 
Spitzer c2d observations. We imposed certain restrictions on the shape of the 3.6 — 70 
SED in order to pick out the embedded objects (the fourth and fifth criteria), and justified 
their validity by comparing to both 1-D and 2-D models of embedded sources. Here we 
further examine the ability of these criteria to identify embedded protostars. 



4.1. Comparison to Known Class Objects 

Class objects were added to the original infrared spectral slope classification system 
of Lada (1987) by Andre et al. (1993) as objects younger and more deeply embedded than 
Class I objects and thus not detected in the infrared. However, with the advent of sensitive 
infrared facilities such as Spitzer, both Class and I objects are detected in the infrared and 
are often indistinguishable from one another based on their infrared spectral slope a (e.g., 
Enoch 2007; Enoch et al. 2008a, in preparation; Kauffmann et al. 2008). Robitaille et al. 
(2006) discussed the distinction between SED "class" and physical "stage". Physically, the 
distinction between a Stage and Stage I object is that a Stage object still has more than 
half of its total mass in the circumstellar envelope (Andre et al. 1993). Various evolutionary 
indicators can be used to separate embedded objects into Class and Class I objects, which 
should correlate with the true physical stage (see §6.2.2p . 

Froebrich (2005) searched the literature and compiled a database of photometry of all 
known Class objects. This database is available online§ and is updated as new Class 
objects are discovered and new photometry is published. If we remove our sixth criterion, 
which selects only those embedded protostars with L int < 1 L , our criteria should select all 
of the objects in this database with available photometry between 3.6 — 70 fim. 

Since both IRAS and the Infrared Space Observatory (ISO) obtained photometry at 25 
and 60 /im, there are many objects in this database with photometry at these two wavelengths 



3 Available at |http : / / astro . kent . ac . uk/protostars / index . html 
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but not at 24 and 70 /im (in other words, observed with IRAS and/or ISO, but not with 
Spitzer). Thus, we first selected all sources in the database that had a detection at 60 or 
70 /im; 80 of the 135 sources in the database meet this criterion. We then imposed the 
same criteria as described in §3.11 and 13.3} omitting the sixth criterion as described above, 
repacing 24 and 70 /im with 25 and 60 /im, respectively, when necessary, and omitting the 
second, third, and seventh criteria because these require information not available from the 
database. 

Of the 80 objects evaluated, 79 pass our criteria and are selected as embedded objects. 
The one source that does not is IRAS 12500—7658, located in the Chamaeleon II molecular 
cloud. According to the Class database, this source has F 24 = 8.58 x 10~ n ± 0.05 x 10~ n 
erg s _1 cm -2 and F 7Q = 6.98 x 10~ n ± 0.43 x 10 -11 erg s" 1 cm -2 , and thus does not have 
a rising (or flat) flux between 24 or 70 /im, even when the uncertainties are taken into ac- 
count. This source was observed by c2d, and in fact these values are taken from Young 
et al. (2005), who presented MIPS data on the Chamaeleon II cloud using a preliminary 
version of the c2d pipeline. In the final c2d source catalog for Chamaeleon II, which in- 
cludes, among other things, improved source photometry and more accurate photometric 
uncertainties (Evans et al. 2007), this source has F 24 = 7.64 x 10~ n ± 0.71 x 10~ n erg s _1 
cm -2 and F 70 = 6.60 x 10~ n ± 0.64 x 10 -11 erg s -1 cm -2 . Using these new, more accurate 
values, IRAS 12500—7658 does indeed have a rising (or flat) SED between 24 or 70 /im when 
the uncertainties are taken into account. Thus, all 80 sources are recovered; we identify all 
of the Class objects. 

4.2. Comparison to Embedded Objects in Perseus, Serpens, and Ophiuchus 

The above exercise demonstrates that our criteria successfully identify the youngest, 
most heavily embedded protostars. However, we are interested in the full population of 
embedded objects, not just the Class objects. To test our method for Class I objects, and 
to verify the above results, we compare our objects to those identified by Enoch (2007) and 
Enoch et al. (2008a, in preparation). They searched for Spitzer sources associated with 
dense cores identified by large-scale, uniform 1.1 mm Bolocam dust continuum emission 
surveys of Perseus, Serpens, and Ophiuchus. They did not restrict themselves to the low- 
luminosity population, and since they had good-quality millimeter data available for all three 
clouds, they started by searching for sources associated with known dense cores. Our work 
is complementary in that all sources identified here should be identified by Enoch et al., but 
we can extend our analysis to regions lacking such millimeter data. 

Comparing the results of both searches shows, in general, good agreement. Enoch et al. 
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separated their objects into Class and Class I objects based on the bolometric temperature 
of each source (see §6.2p . After removing sources with Ljr > 0.5 L from their sample, we 
identify all of their Class sources. This confirms the above result that we identify all Class 
objects. We also identify most of their Class I sources, but find a population of objects 
in their sample not identified as candidates in our work. Inspection of this population of 
objects shows that they are not selected by our criteria because they do not have rising 
fluxes from 24 to 70 /im. Out of the 37 Class I objects with Ljr < 0.5 L in the Enoch 
(2007) and Enoch et al. (2008a) sample, 15 (40%) are not selected by our criteria for this 
reason. These 15 objects range in luminosities from 0.03 < Lm < 0.4 L . Figure [8] shows 
the average SEDs, weighted by l/L^i so that the average is not dominated by the most 
luminous sources, for all sources identified by Enoch (2007) and Enoch et al. (2008a, in 
preparation), divided into four bins: Class 0, Class I with rising (or flat) fluxes from 24 to 
70 /im, Class I with decreasing fluxes from 24 to 70 /im, and Class II. 

Following Whitney et al. (2003b), Enoch et al. (2008a, in preparation) use the bolo- 
metric temperature of each object to go beyond the simple Class 0/1 division and divide 
their sources into Early Class 0, Late Class 0, Early Class I, and Late Class I objects. They 
present the average SED for each of these four classes of objects; it is the Late Class I ob- 
jects (300 K < Tb i < 650 K) that feature decreasing fluxes between 24 and 70 /im. Thus, 
a possible conclusion is that our selection criteria miss the older Class I objects nearing the 
end of the embedded phase. 

However, as described in §3.1} we developed our selection criteria based on comparison 
to models of protostars still embedded within their dense cores. This includes both Class 
and Class I protostars. Two-dimensional models featuring outflow cavities predict rising or 
flat fluxes from 24 to 70 /im at even the most extreme pole-on inclinations, (e.g., Whitney 
et al. 2003a; Whitney et al. 2003b; Robitaille et al. 2006). Indeed, Enoch et al. (2008a, 
in preparation) compare their average SEDs for Early and Late Class and I objects to 
the two-dimensional models of each type of object presented by Whitney et al. (2003b) and 
note that their Late Class I objects differ from the models in that, among other things, they 
decrease in flux from 24 to 70 /im. 

The nature of these objects is uncertain. The overall similarities between the average 
SEDs of the Class I objects rising and falling between 24 and 70 /im (see Figure [8]) and the 
fact that these objects are associated with dense cores suggests they truly are embedded 
objects. On the other hand, no published models of embedded protostars predict fluxes that 
decrease between these two wavelengths. Robitaille et al. (2006) noted that the SED class 
and physical stage of an object do not always agree. For example, a Stage II object located 
behind a large amount of material (such as a dense core) could feature a Class I SED. It is 
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possible that these objects with decreasing fluxes from 24 to 70 /im are not truly embedded 
protostars but heavily extincted Stage II objects located behind dense cores. A future, 
detailed study of these objects is needed, but at present, their true physical stage remains 
unknown. Our criteria will identify the majority of embedded, low-luminosity objects (100% 
of the Class objects; > 60% of the Class I objects). Even if the objects with decreasing 
fluxes from 24 to 70 /im are in fact embedded objects, including them in our sample would 
not significantly change any of our results. 



5. Confirming the Candidates 

Many of the 218 candidates identified in $3] are not truly embedded protostars. The 
visual inspection described above removed a large number of false candidates, cutting our 
list down from 726 candidates to 218. However, we only removed sources from consideration 
if it was obvious that they fell into one of the five types of false candidates discussed in §3.41 
Candidates for which there was any ambiguity were kept, and at least some of these are 
expected to be false candidates. Furthermore, as discussed above, a source is only evaluated 
by the classification method of Harvey et al. (2007b) for identifying candidate YSOs and 
candidate galaxies if it is detected in both epochs of observations (or in the single epoch 
if only one epoch was observed) with SNR > 3 in all 4 IRAC bands and the 24 /im MIPS 
band. Any source not covered by all 5 of these bands is ineligible for classification as either 
candidate galaxy or candidate YSO, and instead receives one of several general source types 
based on the shape of the SED (Evans et al. 2007). Many galaxies thus slip into our 
candidate list because they pass criterion 7 of §3.11 whereas they would have been classified 
as candidate galaxies and rejected according to this criterion if they had been observed in all 
5 bands. Only those that are large and/or close enough to resolve are removed by the visual 
inspection described in §3.41 This is particularly relevant for the five large clouds, since the 
MIPS coverage is substantially larger than the IRAC coverage. 

To demonstrate the extent of extra-galactic contamination expected, we note that 148 
of the 218 candidates have Lir < 0.05 L , which roughly corresponds to L int < 0.1 L using 
the relation between the two derived in §3.11 This suggests a very large number of VeLLOs 
exist; however, out of these 148 sources, 114 are not considered by the Harvey et al. (2007b) 
classification method because they do not meet the requirement of SNR > 3 detections in 
all 5 bands. As Figure M demonstrates, the majority of these low-luminosity sources are 
expected to be galaxies. Out of 851 sources in the 82 regions with dense cores observed by 
c2d that are classified as either candidate YSOs or candidate galaxies, 604 have Lm < 0.05 
L Q , and 518 of these 604 (~ 86%) are classified as candidate galaxies. Thus, as many as 
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98 of the 114 candidates not considered by the Harvey et al. (2007b) classification method 
because they do not meet the requirement of SNR > 3 detections in all 5 bands may be 
extra-galactic in nature. Clearly, further effort is required in order to determine which of the 
218 candidates are in fact embedded, low-luminosity protostars. 

5.1. How to Prove a Source is Embedded 

We examined both the Spitzer data and complementary data for each of the 218 can- 
didates to search for evidence proving they are truly embedded objects. We divide this 
evidence into two "levels of certainty" : 

1. Evidence exists proving that a source is associated, in projection on the sky, with a 
region of high volume density. 

2. Evidence exists proving that this soure is actually embedded within the dense region 
(to remove the possibility of chance alignment). 

As discussed below, data on molecular line and dust continuum emission, molecular outflows, 
extinction, and infrared nebulosity are used to evaluate these conditions. We note here that 
it is much easier to prove the first level than the second. 

5.1.1. Evaluating the First Level 

The first level is evaluated by observing the region in a tracer of dense material. Such 
tracers include millimeter dust continuum emission, which typically traces regions with vol- 
ume densities of n > 10 4 cm~ 3 (e.g., Enoch et al. 2007), inversion transition lines of NH 3 , 
which also trace regions with volume densities of n > 10 4 cm -3 (e.g., Ho & Townes 1983; 
Benson & Myers 1989), and rotational transition lines of N2H + , which trace similar material 
as the NH 3 lines (Caselli et al. 2002). We thus looked for associations with millimeter dust 
continuum, NH 3 , or N 2 H + sources using a number of different surveys of the regions covered 
by c2d in these tracers, some of which were specifically designed as complementary surveys 
to c2d; we list these surveys below. In addition, we searched the SIMBAD database for 
other such observations besides those listed below. In all candidate is said to be 

associated (in projection) with a millimeter dust continuum, NH 3 , or N 2 H + source if it is 
located within one beam of the peak position of the source. 
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The surveys we used for the candidates in the 82 regions with dense cores are as follows: 
An on-going 350 /im dust continuum survey of cores with the Submillimeter High Angular 
Resolution Camera II (SHARC-II) at the Caltech Submillimeter Observatory (CSO; Wu 
et al. 2007; M. Dunham et al. 2008, in preparation), a 450 and 850 fim dust continuum 
survey of 38 cores with the Submillimeter Common-User Bolometer Array (SCUBA) at the 
James Clerk Maxwell Telescope (JCMT; Young et al. 2006a), a 1.2 mm dust continuum 
survey of 37 cores with the Max Planck Millimeter Bolometer (MAMBO) at the IRAM 30 
m telescope (Kauffmann et al. 2008), a 1.2 mm dust continuum survey of 151 southern 
cores with the SEST Imaging Bolometer Array (SIMBA) at the Swedish ESO Submillimeter 
Telescope (SEST; K. Brede et al. 2008, in preparation), an N 2 H + (1-0) survey of 38 cores 
observed with the Five College Radio Astronomy Observatory (FCRAO; C. De Vries et al. 
2008, in preparation), an N 2 H + (1-0) survey of 59 cores also observed with the FCRAO 
(Casein et al. 2002), and observations of the (J, K) = (1, 1) and (2,2) lines of NH 3 for 264 
cores compiled from the literature by Jijina et al. (1999). 

The surveys used for Perseus are a 1.1 mm dust continuum survey of 7.5 deg 2 with 
Bolocam at the CSO (Enoch et al. 2006), an 850 /im dust continuum survey of 3.5 deg 2 
with SCUBA at the JCMT (Kirk et al. 2006), a 450 and 850 /im dust continuum survey of 
3.0 deg 2 with SCUBA at the JCMT (Hatchell et al. 2005; Hatchell et al. 2007a), and the 
SHARC-II survey described above, which includes several cores in Perseus. For Chamaeleon 
II, we used a 1.3 mm dust continuum survey of 36 YSOs in Chamaeleon I and II with the 3 He- 
cooled bolometer system at the SEST (Henning et al. 1993), and the SIMBA survey of 151 
southern cores described above, which includes several cores in Chamaeleon II. For Lupus, 
a 1.2 mm dust continuum survey of ~625 arcmin 2 with SIMBA at the SEST (Tachihara et 
al. 2007) was used. For Ophiuchus, we used a 1.1 mm dust continuum survey of 10.8 deg 2 
with Bolocam at the CSO (Young et al. 2006b), a 1.2 mm dust continuum survey of 1.0 
deg 2 with SIMBA at the SEST (Stanke et al. 2006), and an 850 fim dust continuum survey 
of 700 arcmin 2 with SCUBA at the JCMT (Johnstone et al. 2000). Finally, for Serpens we 
used a 1.1 mm dust continuum survey of 1.5 deg 2 with Bolocam at the CSO (Enoch et al. 
2007), a 450 and 850 ^m dust continuum map of 120 arcmin 2 with SCUBA at the JCMT 
(Davis et al. 1999), and a 1.2 mm dust continuum map of ~150 arcmin 2 with MAMBO at 
the IRAM 30 m telescope (Djupvik et al. 2006). 

Observations tracing dense material are not always available at the positions of the 
candidates. This is especially true for candidates in the southern clouds and cores. As an 
alternative, observations indicating high column density can be used to evaluate the first 
level. A positive result does not guarantee this level is satisfied, as high column densities do 
not necessarily indicate high volume densities, but it does increase the likelihood that the 
source is an embedded object. 
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We utilize two indicators of high column density: extinction maps and absorption against 
the mid-infrared background. For the extinction maps, we used the maps created by the 
c2d team. We give a brief description of them here; a more complete description can be 
found in Evans et al. (2007). A line-of-sight extinction estimate is produced by the c2d 
pipeline for each source with a 1.25 — 24 //m SED consistent with that of an extincted stellar 
photosphere, and these line-of-sight estimates were then convolved with uniformly spaced 
Gaussian beams. Maps of different resolution were produced by using Gaussian beams with 
different FWHM. We used the highest resolution maps available for each cloud (180" for 
Perseus, 120" for Chamaeleon, 120" for Lupus I, 90" for Lupus III, 90" for Lupus IV, 240" 
for Ophiuchus, and 90" for Serpens); extinction maps are not yet available for the cores. We 
did not attempt to derive quantitative estimates of the column density from the extinction 
maps; instead, we examined the maps at the positions of the candidates and identified any 
candidate located in a localized region of higher extinction as showing evidence for high 
column density. 

For the absorption against the mid-infrared background, we looked for "dark cores" 
or "shadows" in the 8 and 24 /zm images. We examined both the images themselves and 
radial profiles of the background intensity. Figure [TU] shows an example of a dark core at 
the position of the candidate in L673-7 (source 031). While a quantitative analysis of the 
extinction profile and column density through the core is possible (e.g., Stutz et al. 2007), it 
is beyond the scope of this paper. Here we simply identified any candidates that appeared 
to be located within such a dark core as showing evidence for high column density. An 
important note is that, in both techniques, negative results do not rule out the possibility 
of high column density. Very small, localized column density enhancements may not be 
seen in the extinction maps due to beam dilution, and variations in the strength of both 
the background and foreground emission will have a significant impact on the presence or 
absence of absorption against the mid- infrared background. As a result, we consider cases 
where data are available to evaluate the existence of regions of high column density but 
return a negative result to be equivalent to cases where no such data are available. 



5.1.2. Evaluating the Second Level 

The second level is primarily evaluated by searching the region around a source for a 
molecular outflow that is centered both spatially on the source and kinematically at the 
velocity of the dense core with which the source is associated. Large-scale surveys of 12 CO, 
the primary tracer of molecular outflows, with the necessary sensitivity, spatial resolution, 
and velocity resolution to detect molecular outflows are not as common as surveys of dust 
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continuum. In fact, the only such survey for the regions observed by c2d was a search for 
outflows in Perseus by Hatchell et al. (2007b), and even this was not an unbiased survey 
as they targeted known cores from a previous 850 /xm dust continuum survey (Hatchell et 
al. 2005). Thus, we searched the literature for each of the 218 candidates to find those 
with known molecular outflows. Many candidates had no published 12 CO observations, 
including several candidates in Perseus that were not part of the Hatchell et al. (2007b) 
survey. Some of these were observed at the CSO as part of a search for molecular outflows, 
resulting in detections of outflows by us in L673-7 and L1251 (Sources 031, 044, and 045 from 
this work; M. Dunham et al. 2008, in preparation; J.-E. Lee et al. 2008, in preparation). 
While a positive result will prove a source is embedded in a dense core, a negative result 
will not disprove this since the answer can depend on the available data. As an example, 
L1014 (Source 038 from this work) showed no signs of driving an outflow from single-dish 
molecular line observations (Crapsi et al. 2005a), but a compact, low- velocity 12 CO (2- 
1) outflow was discovered by Bourke et al. (2005) with the Submillimeter Array. Thus, in 
evaluating whether each source shows evidence for being an embedded protostar, we consider 
cases where observations tracing molecular outflows exist but no outflow is detected to be 
equivalent to cases where no such observations exist. 

We also use a 350 /iin SHARC-II detection to satisfy the second level. As with all 
submillimeter and millimeter observations of dust continuum emission, a detection with 
SHARC-II indicates high volume density. However, unlike other observations of the dust 
continuum, Wu et al. (2007) concluded that, through a combination of temperature and 
instrumental effects, a SHARC-II 350 fxm detection effectively always indicates that a core 
has a protostar embedded within it. 

Many of the 218 candidates have not been observed either in 12 CO or with SHARC- 
II at 350 /im. Both a search for molecular outflows and an extension of the SHARC-II 
survey of dense cores presented by Wu et al. (2007) are on-going at the CSO (M. Dunham 
et al. 2008, in preparation), but in the meantime, this leaves us unable to evaluate the 
second level for many candidates. Thus, we also examine the IRAC 2 (4.5 /jm) image for 
each candidate to look for extended nebulosity or jets, suggesting the presence of molecular 
outflows and outflow cavities. We specifically focus on 4.5 /im for two reasons: the wavelength 
is short enough for scattered light off the edges of outflow cavities to often be visible, and the 
photometric band overlaps with emission lines of molecular hydrogen shocked by outflows 
(e.g., Noriega-Crespo et al. 2004). Examples of this are shown for Sources 001 and 004 
from this work in Dunham et al. (2006) and Bourke et al. (2006), respectively. We do 
not separately identify "extended nebulosity" and "jets"; instead we consider evidence for 
either to be the same thing. While this is not strictly correct, it avoids the ambiguity 
sometimes inherent in assigning extended structure one label or the other. Furthermore, 
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both are likely indicators of outflow activity, the actual source property of relevance for this 
study. A positive result does not guarantee that the second level is satisfied, but it does add 
to the likelihood that the source is an embedded object. Once again, however, a negative 
result does not prove that the source is not embedded. To again use L1014 (Source 038) as 
an example, it does not show any extended nebulosity or jet features in the 4.5 /zm image 
(or any other Spitzer image); only in deep near-infrared images is scattered light from an 
outflow cavity seen (Huard et al. 2006). Thus, as above, we consider cases where 4.5 jttm 
Spitzer images are available but show no extended nebulosity and/or jets to be equivalent 
to cases where no such images are available. 

5.2. Assessing the Likelihood of Being Embedded 

Based on the results of evaluating the two levels of certainty described above, we divide 
the candidates into 6 groups of descending likelihood of being embedded protostars. We list 
the groups in Table M and describe them below: 

1. Group 1 consists of candidates that are confirmed as embedded, low-luminosity pro- 
tostars. Existing observations confirm that sources in this group are embedded within 
regions of high volume density. 

2. Group 2 consists of candidates that have a very high likelihood of being embedded, 
low-luminosity protostars. Sources in this group are associated with regions of high 
volume density. They show extended nebulosity and/or jets in 4.5 fim Spitzer images, 
but they are not confirmed to be embedded within the high-density regions, either 
because no observations are available to evaluate this condition or such observations 
are available but give a negative result. 

3. Group 3 consists of candidates that have a high likelihood of being embedded, low- 
luminosity protostars, although not as high as group 2 candidates. Similar to group 
2, group 3 candidates are associated with regions of high volume density but are not 
confirmed to be embedded within these regions. Unlike group 2, however, group 3 
candidates are not associated with extended nebulosity and/or jets at 4.5 //m, either 
because no such extended structure is detected or because 4.5 fim images are not 
available. 

4. Group 4 consists of candidates that show evidence for being embedded, low-luminosity 
protostars but lack confirmation that either level of certainty is fulfilled. Sources in this 
group have no available observations to evaluate whether or not they are associated 
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with regions of high volume density and either no available observations or observations 
giving a negative result on whether or not they are embedded within regions of high 
volume density. These sources either are associated with regions of high column density 
and show extended nebulosity and/or jets in 4.5 /zm Spitzer images (group 4a), are 
only associated with regions of high column density (group 4b), or only show extended 
nebulosity and/or jets in 4.5 /j,m Spitzer images (group 4c). 

5. Group 5 consists of candidates that might be embedded, low-luminosity protostars, 
but show no sign of being such objects in their limited available data. Sources in this 
group have no available observations to evaluate whether or not they are associated 
with regions of high volume density and either no available observations or observations 
returning a negative result for proving they are embedded within dense regions, evalu- 
ating their association with regions of high column density, and identifying nebulosity 
or jets at 4.5 /im. 

6. Group 6 consists of candidates that are most likely not embedded, low- luminosity 
protostars. Observations show that sources in this group are not associated with regions 
of high volume density, with the important caveat that most dust continuum emission 
surveys used to search for regions of high volume density are only complete to cores 
with masses > 0.1 - 1.0 M (e.g., Enoch et al. 2007). 

The last column of Table [9] lists the total number of the 218 candidates in each of the 
groups. Table [10] lists, for all 218 candidates, sorted by group, the source number, cloud 
or core region in which the source is located, Lm, the status of 160 fim observations for 
this source, and the SIMBAD name of the source, where applicable. If a candidate is not 
associated with an infrared or radio point source in SIMBAD, we give the name of the 
nearest extended dense core within a search radius of 1'. Note that this does not guarantee 
the candidate is embedded within this core. 

We find that 149 of the 218 candidates are in either group 5 or group 6 and thus show no 
evidence for being embedded, low-luminosity protostars. A fraction of these 149 candidates 
are likely real embedded sources that simply lack complementary data of sufficient quality 
to prove this; these sources should be re-evaluated as additional surveys of these regions 
become available. Most, however, are likely either galaxies or more evolved sources no 
longer embedded within dense cores, consistent with our earlier claim that as many as 98 of 
the candidates may be extra-galactic in nature ($HJ). This leaves 69/218 candidates that, at 
minimum, show some evidence for being embedded objects. We are able to obtain complete 
SEDs for 50 of these 69 candidates; these SEDs are listed in Table [TTJ We consider only 
these 50 sources in the rest of the paper. Of the 19 remaining candidates, 16 are located 
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close enough to other sources that they cannot be separated into individual sources at 70 
/im; these are discussed below. Two sources, 027 and 202, are not covered by the 70 jLtm 
c2d observations and have no obvious association with IRAS 60 /mi sources. As 70 jum is 
an essential wavelength for this study (see below), we do not list these objects in Table [TT1 
The last remaining source is Source 132, which is identical to Source 017 as it was observed 
by both IRAC and MIPS in the map of the dense core DC303.8-14.2 and by MIPS in the 
map of Chamaeleon II. We thus remove Source 132 from consideration to avoid duplication. 

5.3. Effects of Confused 70 /im Data 

Good quality 70 /mi data are essential to this study. Photometry at this wavelength is 
important for distinguishing between SEDs of embedded protostars and more evolved YSOs 
no longer embedded in dense cores ( §3. It §4.2j) . and for deriving accurate protostellar internal 
luminosities ( §3.21) . However, the MIPS 70 /mi PSF is ~ 18", compared to 6" at 24 /mi and 
2" at 3.6 — 8.0 /mi. Embedded protostars spaced by more than 6" but less than 18" will 
be resolved into individual sources at 24 /mi and thus identified by our selection criteria, 
but confused at 70 /im. Various methods exist for dividing the total 70 /mi flux among the 
confused sources, such as by the ratio of 24 /im fluxes of each source (e.g., Lee et al. 2006), 
but all such methods introduce large uncertainties. 

Of the 69 candidates that show at least some evidence for being embedded objects, 16 
are located too close to other objects to resolve into individual sources at 70 /mti Many 
likely have L int > 1.0 L and are thus not relevant to this study, since most have values of 
L IR , integrated only to 24 /mi, greater than 0.1 L . However, a few might be embedded, 
low-luminosity protostars. We do not list any of these 16 candidates in Table [TT] and do 
not consider them in the analysis discussed below. Higher-resolution far-infrared data (e.g., 
Herschel) will be necessary to evaluate the nature of these objects. 

6. Discussion 

6.1. Estimating Internal Luminosities and the Luminosity Distribution 

The internal luminosity is a key parameter in understanding the evolution of an object 
from dense core to star, but for low-luminosity, embedded protostars, it is not a directly 



4 The source numbers are 002, 006, 027, 048, 049, 055, 056, 073, 208, 210, 213, 214, 215, 108, 211, 212. 
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observable quantity since a significant component of the bolometric luminosity will arise 
from external heating by the ISRF. In §31 we presented two methods of estimating L int based 
on observable quantities. The first estimate is based on the result that an approximately 
constant ratio exists between L int and L IR , where L IR is an observable quantity. We call 
this estimate of the internal luminosity L\^ t . The second estimate is based on the tight 
correlation observed between the 70 /im flux, scaled to the value it would have if the source 
were located at 140 pc, and the internal luminosity. This estimate, which we call L™ f , is 
given by Equation [2J 

Figure [HI which plots L\^ t vs. Lj° t , demonstrates that the two estimates agree to within 
a factor of two of each other for nearly all of the sources. We consider L™ t to be a more 
accurate estimate of L int , both because nearly identical relationships between L int and 70 
/im flux normalized to 140 pc were derived separately from observations of protostars and a 
large grid of models and because L\^ t is based on a relationship between L int and Lj R that 
shows some variation from source to source. Throughout the rest of this paper we use Lj® t 
when we refer to internal luminosities. Detailed radiative transfer modeling of each source 
is necessary to obtain more accurate internal luminosities. 

Figure [121 shows the distribution of Log(L ini ) for the 50 objects listed in Table [TT1 The 
decrease in the number of protostars near Li nt ~ 1 L Q is an artifact introduced by requiring 
that all sources have Lm < 0.5 L , which is roughly equivalent, but not identical, to requiring 
that all sources have L int < 1.0 L given that the ratio of L int /L IR has a weighted average 
of 1.7 but varies from source to source. The full distribution of protostellar luminosities will 
be examined in a future paper (Evans et al. 2008, in preparation) that combines this work 
with other, complementary searches for embedded protostars in the c2d sample (e.g., Enoch 
2007; Enoch et al. 2008a, in preparation; J0rgensen et al. 2007; J0rgensen et al. 2008, in 
preparation). Here we only focus on the low end of the distribution. 

We find that 15/50 (30%) of the objects identified here have L int < 0.1 L and are thus 
classified as VeLLOs. Qualitatively speaking, this indicates that there are more embedded 
protostars at lower luminosities than at higher luminosities (if they were distributed evenly 
in luminosity, only 10% would have L int < 0.1 L ). To quantify this, Figure [TBI shows the 
distribution of source internal luminosities, cut off at L int = 0.5 L since our sample is likely 
incomplete near 1 L . A visual inspection of this figure suggests that the distribution is not 
uniform but increases with decreasing luminosity. Applying a K-S test, we find that there is 
a 94.2% probability that these sources are not drawn from a uniform luminosity distribution. 
The implications of this increase in number of protostars with decreasing luminosity will be 
explored in §6.41 



Six VeLLOs have internal luminosities lower than our sensitivity limit for the most 
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distant regions (~ 5 x 10~ 2 L Q ), raising the question of how many such sources we miss in 
the more distant regions of our sample. A discussion of correction for completeness due to 
the non-uniform distances to observed regions is presented in Appendix A. Based on this 
discussion, we conclude that only ~ 2 sources with luminosities between our sensitivity limit 
for the closest regions (L irit = 2.5 x 10~ 3 L ) and the most distant regions (L int = 5x 10~ 2 L ) 
are missed. Also based on Appendix A, we conclude that the intrinsic luminosity distribution 
does not continue to increase below L int = 0.1 L , although there are sources present all the 
way down to the sensitivity limit. Finally, we note that we can not draw any conclusions 
about the possible presence of objects with extremely low luminosities below ~ 1CT 3 L . 
Future infrared facilities with sensitivities much greater than that of Spitzer will be required 
to search for such objects and determine the true lower limit to the intrinsic luminosity 
distribution. 

The above conclusions are based on the analysis presented in Appendix A. This analysis 
assumes that all sources with luminosities above our sensitivity limit of Li nt = 4 x 10~ 3 
(d/140pc) 2 L are in fact identified. To verify that this is indeed the case, we re-examine 
the criteria described in £j3l The only criterion that could potentially filter out sources with 
very low luminosities is the requirement that a source not be classified as a candidate galaxy 
in the classification method of Harvey et al. (2007b) unless Log(P ga i) < —1.25, since several 
of the steps in this classification method are based on the general fact that galaxies are faint 
(see Harvey et al. for details). 

Using Perseus as an example, we re-applied our selection criteria, leaving out the step 
that filters out candidate galaxies with Log(P ga i) > —1.25. This results in only an additional 
36 sources, since most of the additional sources that would otherwise be selected are rejected 
because they are not detected at 70 /zm. We visually inspected these sources and rejected 
14 based on either showing obvious galaxy morpologies or not being real sources (see §3.41 
for details). This left us with 22 sources. The positions of all 22 sources were covered by at 
least one of the three sub millimeter/millimeter dust continuum emission surveys of Perseus 
listed in §5.1.11 but none were associated with dense cores. All 22 sources would thus be 
placed into group 6 in the terminology of §5.21 We conclude that we do indeed identify all 
sources above our sensitivity limit. 

6.2. Evolutionary Indicators 

For the 50 objects with complete SEDs listed in Table [TTJ we calculate the following 
quantities: the bolometric luminosity (Lboi), the ratio of bolometric to submillimeter lumi- 
nosity {Lboil Lsmm), and the bolometric temperature (T bo i). L boi is calculated by intergrating 
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over the full observed SED, 

P oo 

L bol = And 2 / S v dv , (4) 
Jo 

while the submillimeter luminosity is calculated by integrating over the observed SED for A 
> 350 /im, 

rv=c/350 /im 

L S mm = 47T(i 2 / S u du . (5) 

Jo 

The bolometric temperature is defined to be the temperature of a blackbody with the same 
flux- weighted mean frequency as the source (Myers & Ladd 1993). Following Myers Sz Ladd, 
T bo i is calculated as 

T bol = 1.25 x 10- 11 Jo K . (6) 

Jo S » dv 

We list L^i, L bo i/L smm , and T bo [ for these 50 objects in the last column of Table dH The 
stated uncertainties reflect only those arising from measurement error. Since we are inte- 
grating over finitely sampled SEDs, the values calculated for all three quantities will depend 
on how well each SED is sampled. We discuss the errors introduced by incomplete sampling 
in Appendix B. 



6.2.1. Bolometric Luminosity-Temperature (BLT) Diagram 

Figure [TH places the 50 objects listed in Table [11] on Figure 19 from Young & Evans 
(2005), which plots L boi vs. T boi . As discussed by Myers et al. (1998), who refer to this 
as a Bolometric Luminosity- Temperature (BLT) diagram, this is effectively a Hertzsprung- 
Russell diagram for protostars. The lines show the evolutionary tracks followed by various 
models from Young & Evans (2005) and Myers et al. (1998). In general, models predict 
that protostars form at low values of L bo i and T bo i and move up and to the left (increasing 
L bo i, increasing T bo i) as the protostar grows. Whether or not L bo \ steadily increases or 
begins to decrease depends on the details of the mass accretion process. Young & Evans 
considered the collapse of a singular isothermal sphere in the "standard model" of star 
formation (Shu, Adams, & Lizano 1987), which features a constant mass accretion rate over 
the entire protostellar phase. Myers et al. (1998), on the other hand, assume that the 
accretion rate decreases exponentially with time and mass is lost from the system (from 
outflows, winds, etc.), explaining the different trends shown by their evolutionary tracks. 

The colors/symbols of the points in Figure [T41 for the 50 objects from this work reflect 
how well the source SED is sampled. All four colors/symbols have SEDs that are sampled 
at 1.25 — 70 /im by 2MASS and Spitzer and at millimeter wavelengths. Black circles are also 
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sampled at, at minimum, one wavelength in both the 100 — 200 and 350 — 450 /im wavelength 
ranges (category 1 SEDs in the notation of Appendix B). Blue triangles are sampled between 
350 — 450 fim but not 100 — 200 fim (category 2 SEDs). Green diamonds are sampled between 
100 — 200 fj,m but not 350 — 450 /im (category 3 SEDs). Red squares are not sampled in either 
wavelength range (category 4 SEDs). A detailed analysis of the errors introduced by the 
incomplete sampling in these various categories is presented in Appendix B; here we simply 
note that, depending on the category, the error in L boi and T bot is between approximately 
20 - 60%. 

There are many Class I sources (Tt> i > 70 K; Chen et al. [1995]) with Lboi below 
any evolutionary track by up to an order of magnitude. The most extreme cases are all 
sources with category 4 SEDs, where L ho \ is underestimated by ~ 50%, on average. However, 
even accounting for this, many are still substantially below any of the evolutionary tracks. 
A similar result was found by Enoch (2007) and Enoch et al. (2008a, in preparation), 
who suggest episodic accretion as a potential solution. Another potential explanation for 
this discrepancy is that the sources with the lowest L^i are simply those with envelope 
masses below those considered in the models. However, this is ruled out by Enoch et al., 
who find good agreement between observations and models when comparing envelope mass 
vs. bolometric temperature despite the scatter seen in BLT diagrams, as well as several 
individual studies of VeLLOs that conclude at least some of the lowest-luminosity protostars 
are embedded in envelopes with masses similar to those of more luminous protostars (a few 
M ; Young et al. 2004; Bourke et al. 2006; Dunham et al. 2006). 

A similar population of Class objects with luminosities below model predictions is 
not as apparent from Figure [TH and is also not seen by Enoch et al. However, on average, 
external heating from the ISRF will add a larger contribution to Lf, i for objects with more 
massive envelopes, thus such heating is likely to be more significant for Class objects than 
Class I objects and can add up to several tenths of a solar luminosity to L bot (e.g., Evans et 
al. 2001). Since, as described in JQ, this distinction between Lu and L int is most relevant for 
very low-luminosity sources where external heating from the ISRF can dominate the total 
observed L hoh a similar population of Class objects with very low luminosities may be 
present but not readily apparent when examining a plot of L hoi vs. T ho \. We return to this 
point in §6.41 

6.2.2. Comparison of two Indicators of Evolutionary Status 

Both Tboi and Lboi/L smm are indicators of the evolutionary status of a protostar. As 
the envelope mass decreases through a combination of accretion and dissipation processes, 
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both quantities should increase. Eventually, once the star has fully formed, T bo [ will become 
equal to T e ff and L bo i/L smm will become very large. We show a plot of L bo i/L smm vs. T bo[ 
for the objects considered here in Figure [TBI The division between Class 0/1 objects in T bol 
is from Chen et al. (1995), while the division between Class 0/1 and Class I /II objects in 
L bo i/L smm is from Young & Evans (2005). The colors/symbols of the data points hold the 
same meaning as in Figure [TH 

As stated in §4.11 the distinction between a Stage and Stage I object is that a Stage 
object still has more than half of the total mass of the system in the envelope, whereas less 
than half the total mass remains in the envelope of a Stage I object. Both T bo i and L bo i/L smm 
are used to separate embedded objects into Class and Class I objects; in principle, these 
observational classes should correspond to Stage and Stage I objects. Young & Evans 
(2005) concluded that L bo i/L smm is a better measure of physical stage than T boh since their 
evolutionary models of different initial mass cores featured similar values of L bo i/L smm but 
different values of T bo i when they reached the point at which half of the total mass had 
accreted onto the protostar. They also found that embedded objects cross the Class 0/1 
boundary in T bo \ {T bo \ = 70K; Chen et al. 1995) when they are still Stage objects. We 
see a similar trend in Figure [TBI in that several objects are Class I objects according to T bol 
but Class objects according to L bol /L smm . However, we caution that the details of where 
the observed SED is sampled affect the calculated value of L bo i/L smm significantly more 
than the calculated value of T bo i; this explains the objects seen in Figure [15] with extremely 
high values of L bo i/L smm (see Appendix B). Great care must be taken to consider the errors 
introduced by incomplete sampling before using L bo i/L smm as an evolutionary indicator for 
any particular object. 



6.3. Fraction of Starless Cores that Remain Starless 

Out of the 95 individual dense cores observed by c2c§|, 67 (~70%) were classified as 
starless prior to being observed by Spitzer. The very first starless core observed, L1014, was 
found to actually harbor a very low-luminosity, embedded protostar [L int ~ 0.09 L ; Young 
et al. 2004; Bourke et al. 2005; Huard et al. 2005). This naturally gives rise to the question 
of how many cores classified as starless prior to the launch of Spitzer remain so after being 



5 Evans et al. (2007) list the 82 core regions surveyed by c2d. Those regions that covered multiple cores are 
as follows (see also Evans et al. 2003 and references therein): CB130-3 (CB130-1, CB130-2, and CB130-3), 
CG30-31 (CG30, CG31A, CG31B, CG31C), IRAM04191+1522 (IRAM04191+1522 and IRAS04191+1523), 
L1251 (L1251A, L1251C, and L1251E), L43 (L43-East and L43-RN091), L673 (L673-SMM1 and L673- 
SMM2), and TMC1 (TMC1-1, TMC1-2, TMC1-1C-1, and TMC1-1C-2). 
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observed with Spitzer! A substantial number of protostellar cores misclassified as starless 
cores would result in an incorrect estimate of the lifetime of the starless core phased, as such 
estimates are derived by comparing the numbers of starless and protostellar cores (e.g., Kirk 
et al. 2005; Enoch et al. 2008b). Kirk et al. (2007) present the results of a search for 
embedded, low-luminosity protostars in Spitzer observations of 22 starless cores previously 
studied using submillimeter data (Kirk et al. 2005); they find only one such source, L1521F- 
IRS (Bourke et al. 2006). This suggests that ~95% of starless cores remain starless, leading 
them to conclude that the errors in starless core lifetime estimates introduced by previously 
undetected, low-luminosity protostars are smaller than the uncertainties in such estimates. 

Our sample of 67 starless cores allows us to test these conclusions with a dataset approx- 
imately three times larger than that of Kirk et al. (2007). Searching all 95 cores considered 
here, we have identified a total of 49 candidate embedded, low-luminosity protostars in 33 
different cores (Table [3]). Seventeen of these cores were classified as starless prior to being 
observed by Spitzer. If we restrict ourselves only to those sources in Groups 1 — 4, those 
that, at minimum, show some additional evidence of being embedded protostars besides the 
shapes of their SEDs, we identify 29 candidates in 21 different cores, 9 of which were classified 
as starless prior to being observed by Spitzer. Thus, between 9 and 17 out of the 67 starless 
cores actually harbor embedded, low-luminosity protostars. Approximately 75 — 85% of star- 
less cores remain starless down to our luminosity sensitivity of L int > 4 x 10~ 3 (d/140pc) 2 . 
Most starless cores do indeed remain starless after being observed by Spitzer, in agreement 
with Kirk et al. (2007), although the fraction is smaller than they found with their smaller 
sample. 



6.4. Non-Steady Mass Accretion 

In this work, we have identified a sample of 50 embedded protostars with Li nt < 1.0 
L ; 15 of these objects have L int < 0.1 L Q and are thus classified as VeLLOs. Assuming 
spherical mass accretion at the rate predicted by the standard model (M acc ~ 2 x 10~ 6 
M Q yr _1 ; Shu, Adams, & Lizano 1987) onto an object with a typical protostellar radius of 
R ~ 3 R , a protostar located on the stellar /substellar boundary (M = 0.08 M ) would 
have an accretion luminosity, L acc = GM ^ facc , of L ~ 1.6 L Q . The objects identified here 
are difficult to understand in the context of the standard model of star formation. VeLLOs, 
with luminosities more than an order of magnitude lower than the above calculation, are 



6 In actuality, this is an estimate of the lifetime of the portion of the starless core phase detectable by 
(sub) millimeter observations. This is often referred to as the pre-stellar core phase (Kirk et al. 2005). 
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an extreme example of this problem. If the objects identified here were observed edge-on 
through circumstellar disks their luminosities could possibly be underestimated, but, for at 
least some of them, this possibility can be eliminated (e.g., IRAM04191, Dunham et al. 
2006). Thus, they must either feature mass accretion rates far below that predicted by the 
standard model, masses far below the stellar/substellar boundary, or some combination of 
the two. Several authors have invoked episodic accretion to explain the existence of these 
low-luminosity objects (e.g., Kenyon & Hartmann 1995; Young Sz Evans 2005; Enoch 2007; 
Enoch et al. 2008a, in preparation). 

Strong evidence for non-steady mass accretion exists for the VeLLO IRAM 04191-IRS. 
This object drives a well-defined, bipolar molecular outflow (Andre et al. 1999). Taking the 
average mass accretion rate implied by the outflow and the dynamical time of this outflow 
((M acc ) ~ 5 x 10~ 6 M yr _1 and t d ~ 10 4 years), Andre et al. (1999) calculate that a 
protostellar mass of 0.05 M will have accreted over the lifetime of the outflow. Accretion at 
this rate onto an object with a mass of 0.05 M and a radius of 3 R would give rise to an 
accretion luminosity of L acc ~ 2 L Q . This is clearly inconsistent with the internal luminosity 
of L int ~ 0.08 L derived by Dunham et al. (2006). Dunham et al. suggested non-steady 
mass accretion to resolve this inconsistency: the current mass accretion rate must be much 
lower than the average rate over the lifetime of the outflow. 

We show the location of the 15 VeLLOs in our sample on the BLT diagram in Figure 
[TBI Most are located at T bo i < 70 K and L bol > 0.1 L Q and are thus Class objects with very 
low internal luminosities but bolometric luminosities that generally appear consistent with 
evolutionary tracks. If there is a population of Class objects with luminosities lower than 
predicted by evolutionary tracks, as is the case for Class I objects (see §6.2. ip . the inconsis- 
tency between the luminosities of these Class objects and the predictions of evolutionary 
tracks will likely be masked on a BLT diagram since, as described in §6.2.11 the contribution 
to Lboi from external heating will be larger for younger objects with more massive envelopes. 
In addition to L 0O i and Tf, i, considering time as a third dimension relevant to protostellar 
evolution can help illuminate whether or not these Class objects with Li nt < 0.1 L but 
Lboi > 0.1 L are truly as consistent with constant mass-accretion evolutionary tracks as 
they appear to be on a BLT diagram. 

Assuming a constant rate of mass accretion, the protostellar mass M will grow linearly 
with time. If we assume the protostellar radius remains fixed, L acc oc M and thus also 
grows linearly with time. Adding in the assumption that star formation has been occuring 
continuously over a time much longer than the lifetime of the embedded phase, the number 
of protostars as a function of luminosity should be constant. However, as discussed in §6.11 
the intrinsic luminosity distribution is not uniform and instead increases with decreasing 
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luminosity. Even though many of the lowest luminosity objects appear consistent with the 
location of evolutionary tracks calculated assuming constant mass accretion, there are too 
many relative to higher luminosity sources to explain by constant mass accretion. 

We have presented evidence for non-steady mass accretion in low-mass protostars, both 
from the above statistical argument and from a detailed study of one individual source (Dun- 
ham et al. 2006). While we note here that non-steady mass accretion does not necessarily 
imply episodic mass accretion, these results, combined with other studies that find clear evi- 
dence for episodicity in outflow activity (e.g., HH211; Lee et al. 2007), strongly suggest that 
protostellar mass accretion is episodic in nature. Future work is needed on both observational 
and theoretical fronts to better understand the implications of the objects identified in this 
study on the mass accretion process in low-mass star formation. On the observational front, 
the detailed study of the outflow driven by IRAM04191 must be repeated for other sources 
known to drive outflows, and dedicated outflow searches are needed towards objects lacking 
such data. Additionally, detailed chemical and physical studies of the objects identified in 
this work are needed to explore the nature of these sources beyond the simple evolutionary 
indicators calculated here. Such studies are of particular relevance to understanding the 
evolution of starless cores and identifying those cores on the verge of collapse since VeLLOs 
have been discovered in cores that were classified as starless prior to Spitzer observations 
but were not always those cores believed to be the most evolved and nearest to the onset 
of collapse (Crapsi et al. 2005b; Bourke et al. 2006). Chemical studies can also be used 
to distinguish between low-luminosity protostars that have featured higher past luminosities 
(presumably through higher past mass accretion rates) and those that have always featured 
such low luminosities (e.g., Lee 2008). Candidates identified in this study that show no evi- 
dence for being embedded in dense cores should be re-examined as additional data on each 
source becomes available (e.g., through large-scale surveys with Herschel, SCUBA-2, etc.) 
to ensure as complete a sample as possible. Finally, the full sample of embedded, low-mass 
protostars must be assembled from this work and the various other searches for such objects 
outlined in JO 

On the theoretical front, evolutionary models predicting the observational signatures of 
dense cores forming stars must be revisited. These models must be able to explain both the 
existence of embedded protostars with luminosities much lower than predicted by current 
models and the very large scatter in observed at each value of T^i in Figure [TH Models 
that incorporate episodic mass accretion and the effects of source geometry may provide 
such explanations, but will require higher dimensions than the 1-D models considered by us 
to date. 
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7. Conclusions 

We have conducted a search for all embedded protostars with Lj ni < 1.0 L Q in the 
c2d dataset of nearby, low-mass star-forming regions. We identify 218 candidates from the 
Spitzer data alone; examining all available complementary data for each candidate results 
in a sample of 50 objects that show at least some evidence that they are indeed embedded 
within dense cores. A summary of our major results is as follows: 

• On average, the Spitzer c2d data are sensitive to embedded protostars with L int > 
4 x 10~ 3 (rf/140pc) 2 L , a factor of 25 better than the sensitivity of the Infrared 
Astronomical Satellite (IRAS) to such objects. 

• The 70 /iin flux and internal luminosity of a protostar are tightly correlated. As the 
former is a directly observable quantity but the latter is not, this correlation gives a 
powerful method for estimating protostellar internal luminosities when detailed radia- 
tive transfer models for each source are lacking. 

• Of the 50 objects in our sample, 15 (30%) have L int < 0.1 L and are thus classified as 
VeLLOs. The distribution of source luminosities is not uniform and instead increases 
with decreasing luminosity. Accounting for incompleteness arising from non-uniform 
distances to the observed regions, we find sources down to the above sensitivity limit, 
indicating that the intrinsic luminosity distribution may extend to lower luminosities 
than probed by these observations. Despite this, we are able to rule out a continued 
rise in the distribution below L int =0.1 L . 

• Between 75 — 85% of cores classified as starless prior to being observed by Spitzer 
remain starless down to the above luminosity sensitivity; the remaining 15 — 25% 
harbor low-luminosity, embedded protostars. This is in general agreement with Kirk 
et al. (2007), who examined archival Spitzer data of 22 starless cores and found only 
one to be harboring a low-luminosity protostar. However, with our larger sample size, 
we are able to better constrain the fraction of cores previously classified as starless that 
in fact harbor such objects. We confirm that recent estimates of starless core lifetimes 
(e.g., Kirk et al. 2005; Enoch et al. 2008b) do not feature large errors introduced by 
previously undetected, low-luminosity protostars. 

• The observed luminosity distribution for embedded objects with L int < 1.0 L Q is incon- 
sistent with the simplest picture of star formation wherein mass accretes from the core 
onto the protostar at a constant rate. Combining this result with other studies that 
find clear indications of episodic outflow activity strongly suggests that protostellar 
mass accretion is episodic in nature. 
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We have outlined several avenues of future work that must be pursued now that rela- 
tively complete and unbiased samples of embedded, low-mass protostars are being compiled. 
Only with such future studies can we begin to build a coherent picture of low-mass star 
formation consistent with the growing observational database provided by systematic, large- 
scale surveys of low-mass star forming regions. 
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A. Internal Luminosity Completeness and the Intrinsic Luminosity 

Distribution 

The regions observed by c2d are located at distances ranging from 125 — 500 pc. Using 
the result from §3.11 that our observations are sensitive to embedded protostars with L int 
> 4 x 10~ 3 (<i/140pc) 2 L Q , the effect of these non-uniform distances is that our sensitivity to 
embedded protostars varies between 2.5 x 10~ 3 and 5 x 10~ 2 L . Two important questions 
thus arise: (1) How many objects are lost because of these non-uniform sensitivities (stated 
another way, how many additional objects would we expect to identify if all observed regions 
were located at the distance of the closest regions)?, and (2) What conclusions can be drawn 
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about the intrinsic luminosity distribution? The latter is of particular significance to future 
work aimed at explaining the observed distribution of sources in L bo i-T bo i space through 
evolutionary modeling. 

A.l. Internal Luminosity Completeness 

The regions observed consist of both targeted observations of 82 regions with 95 dense 
cores and unbiased surveys of 5 large, molecular clouds Thus, an analysis of complete- 
ness that uses the area observed at different distances is incorrect since the sample consists 
both of observations biased towards dense cores and observations with no such bias. Instead, 
we start from the total number of dense cores observed at different distances since all em- 
bedded protostars are, by definition, within these cores. The sample consists of 122 cores in 
Perseus (Enoch et al. 2006), 35 cores in Serpens (Enoch et al. 2007), and 48 cores in Ophi- 
uchus (Young et al. 2006b), along with the targeted observations of 95 cores. This yields 
a total of 300 cores; the distribution of the sensitivity to protostars embedded within these 
cores is shown in Figure [T71 We have no information on the number of dense cores in the 
Lupus and Chamaeleon II clouds; any such cores have comparable luminosity sensitivities 
to the rest of the sample (Figure [T7|) . 

To start with a simple analysis, we ask what percentage of the cores are located at 
distances close enough such that sources in each of the lowest three bins of Figure [12] could 
be detected. Statistically, this should be equal to the percentage of sources in that bin we 
are able to detect. All of the cores are located at distances such that sources in the bin -1.2 
< Log(L int /L Q ) < -0.9 could be detected, thus we are complete to sources in this bin. 96% 
of the cores are located at distances such that sources in the bin -1.5 < Log(L int /L Q ) < 
-1.2 could be detected; correcting for completeness raises the total number of sources in 
this bin from 5 to 5.2. 90% of the cores are located at distances such that sources in the 
bin -1.8 < Log(L int /L Q ) < -1.5 could be detected; correcting for completeness raises the 
total number of sources in this bin from 4 to 4.4. We conclude that we miss, at most, 1 
source in the range -1.8 < Log(L int /L Q ) < -0.9; this is actually a strong upper limit since 
the percentage of cores with distances close enough to detect sources in these bins would 
increase if the dense cores in the Lupus and Chamaeleon II clouds were included. However, 
this says nothing about the bins that are above our sensitivity limit for the nearest regions 
(2.5 x 10~ 3 L Q ) in which we detect no sources. To examine the completeness in these bins, 
and confirm the above results, we turn to a Monte Carlo simulation. 

For the Monte Carlo simulation, we create a sample of 10,000 embedded protostars. We 
assign an internal luminosity to each protostar from a distribution specified at the start; in 
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all cases, this distribution has a maximum of Li nt = 0.1 L . We then randomly place each of 
the 10,000 sources in one of the 300 dense cores. If a source has L int less than the sensitivity 
limit for the core in which it is placed, it is considered "missed" by our observations and 
rejected; otherwise, the source is kept. We run the simulation 1000 times, average the results, 
and compare the observed distribution of L int for the 15 VeLLOs with the average simulated 
observed distribution. 

Figure dH] shows the results of a simulation where the 10,000 sources have internal 
luminosities evenly distributed (in linear space) between 10~ 5 and 10 _1 L . The normalized 
observed distributions of L{ nt for the 15 VeLLOs identified by this work and for the simulated 
sources agree given the uncertainties due to small number statistics. Furthermore, since we 
detect no sources with Log(L int ) < —1.8, we can set an upper limit of < 1/15 sources (or 
< 7% of all sources detected) have such luminosities. Out of all the sources detected in the 
Monte Carlo simulation, 6.1% have such luminosities, consistent with the observed upper 
limit of 7%. 

Out of the 10,000 sources considered by this simulation, 1158 (12%) have luminosities 
below the sensitivity limit of the core in which they are placed and are thus "missed" by our 
observations. This means that the 15 VeLLOs identified are only 88% of the total number of 
objects, indicating that approximately 2 sources are missed. The exact details of the input 
luminosity distribution do not affect these results, as long as the simulation is consistent 
with the upper limit on the number of detected sources with Log(L int ) < —1.8 described 
above. Input luminosity distributions that have a much larger fraction of sources at lower 
luminosities, such as a distribution flat in log space, would lead to an increase in the number 
of sources missed. However, such distributions are inconsistent with the above upper limit 
(for example, a simulation with an input luminosity distribution flat in log space between 
10~ 3 and 10" 1 L results in 24% of the total sources detected having Log(L int ) < —1.8, 
clearly inconsistent with the observed upper limit of 7%). Thus, we conclude that ~ 2 
sources are missed due to the non-uniform distances of the observed regions. 



A. 2. Intrinsic Luminosity Distribution 

We now consider what constraints we can place on the low end of the intrinsic luminosity 
distribution. Additional constraints will be provided by a future paper aimed at analyzing 
the shape of the full luminosity distribution, including higher luminosity sources (Evans et 
al. 2008, in preparation). These constraints will serve as guides to future work aimed at 
explaining the full distribution of sources in L^-T^ space. 
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As stated in §6.11 the intrinsic luminosity distribution for the objects identified here 
increases to lower luminosities. To examine whether or not this increase continues to our 
sensitivity limit, Figure [Tj|] shows the distribution of internal luminosities, in linear bins 
of 0.01 L Q , for the 15 objects with L int < 0.1 L . Considering the small sample size, 
this distribution appears uniform down to L int = 0.02 L . From the discussion above, 
approximately two sources are missing from this figure because of incompleteness arising 
from the non-uniform distances to the observed regions. Additionally, a strong upper limit 
of 1 source is lost due to incompleteness for L int > 10~ L8 L (denoted by the dashed line 
in Figure [T9|) . suggesting that the sources missing due to incompleteness likely fall in the 
lowest two bins. Taking into account the small sample size, adding ~ 2 sources to these two 
bins results in a distribution that appears approximately uniform all the way down to the 
sensitivity limit (denoted by the dotted line in Figure [T91 . This suggests that the intrinsic 
luminosity distribution may extend to lower luminosities than probed by these observations. 
However, even if all missed sources fall in the lowest bin of Figure [191 this portion of the 
intrinsic luminosity distribution (Li nt < 0.1 L Q ) would appear to be uniform rather than 
increasing to lower luminosities. We thus conclude that the instrinsic luminosity distribution 
does not continue to increase to our sensitivity limit, although K-S tests should be used to 
quantitatively confirm this result once a larger sample becomes available. 

This result provides an important constraint for future attempts at understanding these 
very low luminosity objects. However, we caution that we can draw no conclusions about 
the possible presence of objects with luminosities below ~ 10~ 3 L . To illustrate this with 
an example, Figure [2171 shows the results of a Monte Carlo simulation where half of the 10,000 
sources have luminosities evenly distributed between 10 -4 and 10 -3 L , while the other half 
have luminosities evenly distributed between 10~ 3 and 10 _1 L . The simulated observed 
distribution matches the observed distribution quite well; the population of objects with 
L int between 10 -4 and 10~ 3 L have no effect on the final result since none of them are 
detected. A large population of embedded objects with such low luminosities would be very 
difficult to explain physically (see discussion of VeLLOs in §6.41) . but we cannot rule out the 
presence of such objects based on these observations. 

B. Errors in Evolutionary Indicators From Incomplete, Finitely Sampled 

Spectral Energy Distributions 

The integrals defined in Equations 0] — M are evaluated using the trapezoid rule to 
integrate over the finitely sampled source SEDs. The contribution Si to the total integral 
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from measurements yi and yi+i, measured at Xi and Xj+i, is 

6 i = (x w -x i )(y^^j . (Bl) 

The total value of the integral is then the sum of each individual <5j. 

Clearly, the result will depend on how well the source SED is sampled, especially near 
the peak. Most of the material in the envelope of embedded objects is cold (10 — 15 K; 
e.g., Evans et al. 2001; Shirley et al. 2002). The emission peaks from blackbodies with 
temperatures of 10 and 20 K are located at approximately 500 and 250 /im, respectively, 
overlapping with the ~ 50 — 300 /im far-infrared regime that cannot be observed from the 
ground. Incomplete sampling near the peak may thus introduce significant errors in the 
calculations of L boU T hoU and L bo i/L smm . 

To quantify these errors, we consider the 13 group 1 candidates in Perseus. Most of these 
sources have well-sampled SEDs, with 2MASS and Spitzer detections from 1.25 — 70 /im, 
Spitzer detections at 160 /im and possibly IRAS detections at 100 /im, SHARC-II detections 
at 350 /xm, SCUBA detections at 450 and 850 /im, and Bolocam detections at 1.1 mm. We 
construct four SEDs for each source: 

1. We include all detections. These "category 1" SEDs are those that are the most 
completely sampled. 

2. We only include detections in the wavelength ranges 1.25 < A < 70 /im and A > 350 
/im. These "category 2" SEDs are meant to simulate sources for which 2MASS and 
Spitzer 1 — 70 /im data, 350 or 450 /im submillimeter data, and millimeter wavelength 
data are available, but neither Spitzer 160 /im nor IRAS 100 /im data are available. 

3. We only include detections in the wavelength ranges 1.25 < A < 160 /im and A > 850 
/im. These "category 3" SEDs are meant to simulate sources for which 2MASS and 
Spitzer 1.25 — 70 /im data, Spitzer 160 /im or IRAS 100 /im data, and millimeter 
wavelength data are available, but no 350 or 450 /im submillimeter data are available. 

4. We only include detections in the wavelength ranges 1.25 < A < 70 /im and A > 850 
/im. These "category 4" SEDs are meant to simulate sources for which only 2MASS 
and Spitzer 1.25 — 70 /im and millimeter wavelength data are available. 

For each of the 13 sources, we calculate L^oh Tboi, and Lboi/L smm for each of the four 
versions of the SED. Comparing the results will give an estimate of the errors introduced by 
incompletely sampling the source SED. 
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B.l. Bolometric Luminosity 

The three panels of figure [21] shows the percent error for each source between Lboi 
calculated from the SEDs in categories 2 — 4 and Lboi calculated from the SED in category 
1. For category n, where n = 2, 3, 4, this percent error (PE) is defined as 

PE = 100 x Lbol ' n - Lho ^ (B2) 

LboLl 



For all three cateogories, the calculated Lboi is an underestimate of the Lboi calculated 
from the more completely sampled Category 1 SEDs. The average of the absolute values 
of PE, ignoring values of which simply indicate that the complete SED for the source 
is not a Category 1 SED, are 31%, 25%, and 54% for categories 2, 3, and 4, respectively. 
As expected, including detections either between 100 — 160 or 350 — 450 /mi significantly 
improves the accuracy of Lb u adding detections in both ranges improves it further. 



B.2. Bolometric Temperature 

The results of the same analysis for Tboi as that performed above for Lboi are shown in 
Figure [22j For category 3 and category 4 SEDs, the calculated Tboi is always higher than 
that calculated from a well-sampled SED. This is expected since both categories leave out 
detections near the peak of the SED. For category 2, however, the direction of the error in 
T bo i depends on the true value of T ho i- Category 2 SEDs of the coldest sources (T bo i < 50 
K) will lead to underestimates of Tboi, whereas they will lead to overestimates for warmer 
sources (T bo i > 50 K). 

We calculate the percent error in Tboi in the same manner as above for Lboi- The average 
of the absolute values of these percent errors are 21%, 23%, and 64% for categories 2, 3, and 
4, respectively. As for Lboi-, category 2 and 3 SEDs yield significantly more accurate values 
of Tboi than category 4 SEDs, but are not as accurate as Category 1 SEDs. 

We note here that our results are in good agreement with Enoch (2007), who found that 
overall uncertainties for measured L bo i and T bo i values are approximately 20 — 50%, depending 
on whether or not 160 fim data are available. We also note that distance uncertainties will 
add additional uncertainties to Lw beyond those considered here. 
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B.3. Bolometric to Submillimeter Luminosity 

We performed a similar anaysis for L bo i/L smm as above for L bo i and T bo \. The average 
values of the percent errors are 31%, 1282%, and 605% for categories 2, 3, and 4, respectively, 
significantly larger than the errors in L boi and T bo i. The value of L bo i/L smm calculated for 
category 4 SEDs overestimates the actual value; both L bol and L smm are underestimated by 
not sampling the SED at all between 70 and 850 //m, but L bo i is underestimated less than 
L smm . The same is true for category 3 SEDs, except the error is even larger because the 
calculation of L bo \ is now less of an underestimate (see above), but the calculation of L smm 
remains the same. On the other hand, L bo i/L smm calculated from category 2 SEDs underes- 
timates the actual value, since L bo i is underestimated but L smm is accurately calculated. The 
errors here are more than an order of magnitude smaller than Categories 3 and 4. Clearly, 
great care must be taken to obtain the most well-sampled SEDs before attempting to draw 
conclusions from the value of L bo i/L smm for a given source. 
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Table 1. Embedded Protostars Observed by Spitzer with Published Models 





Lir 


Lint 


Lbol 




Protostar 


(L ) 


(L©) 


(Lq) 


Reference a 


CB130-3-IRS 


0.05 


0.15 ± 0.05 


0.96 


1 


IRAM 04191+1522 


0.02 


0.08 ± 0.04 


0.28 


2 


IRAS 04191+1523 


0.25 


0.4 ± 0.1 


0.64 


3 


IRAS 18148-0440 


8.19 


13.0 ± 0.5 


12.9 


4 


IRAS 19156+1906 


1.13 


2.6 ± 0.5 


3.4 


4 


L1014-IRS 


0.06 


0.09 ± 0.03 


0.3 


5 


L1148-IRS 


0.08 


0.1 ± 0.05 


0.15 


6 


L1221-IRS1 


1.12 


2.5 ± 0.25 


2.7 


7 


L1221-IRS3 


0.36 


1.0 ± 0.15 


1.9 


7 


L1521F-IRS 


0.01 


0.05 ± 0.02 


0.36 


8 


L673-7-IRS 


0.02 


0.04 ± 0.02 


0.13 


3 



a Rcfcrcncc for the radiative transfer model and determination of 
Li n t and !/(,„;: 

(1) H.J. Kim ct al. (2008, in preparation); (2) Dunham et al. (2006); 
(3) M. Dunham et al. (2008, in preparation); (4) Shirley et al. 
(2002); (5) Young ct al. (2004); (6) Kauffmann ct al. (2005); (7) C. 
Young et al. (2008, in preparation); (8) Bourkc et al. (2006.). 



Table 2. Linear Least-Squares Fits to Plots of Log F„ vs. Log L int . 



Wavelength 


Observations 






Models 




H 


m 


b 


Xr 


m 


b 


Xr 


3.6 


0.57 + 0.42 


-11.07 + 0.32 


296 


0.86 ± 0.05 


-11.87 + 0.03 


3309 


4.5 


0.88 + 0.26 


-10.54 + 0.20 


148 


0.94 ± 0.04 


-11.63 + 0.03 


2846 


5.8 


0.97 + 0.29 


-10.54 + 0.22 


179 


0.99 ± 0.04 


-11.43 + 0.03 


2305 


8.0 


0.83 + 0.27 


-10.72 ± 0.22 


160 


1.02 ± 0.04 


-11.38 + 0.03 


2044 


24 


0.87 + 0.20 


-10.05 ± 0.17 


85 


1.38 ± 0.02 


-10.12 + 0.01 


433 


70 


1.06 + 0.06 


-9.26 + 0.06 


3 


1.06 ± 0.01 


-9.02 + 0.01 


15 



Table 3. Candidate Low-Luminosity Embedded Objects in the Cores 



Source 


Spitzcr 


c2d 


Dist. 


RA 


Dec 


L IR d 




Nearest IRAS 


Number Group a 


Source Namc b 


Core Region c 


(PC) 


(J2000) 


(J2000) 


(L©) 


(") 


Source 



001 


1 


J042156. 


88+152946. 





IRAM04191+1522 


140 ± 10 


04 


21 


56 


.88 


+15 


29 


46. 








.023 


43.6 


IRAS 


04191+1523 


002 


1 


J042200. 


.07+153024. 


8 


IRAM04191+1522 


140 ± 10 


04 


22 


00 


.07 


+15 


30 


24. 


.8 


0. 


033 


17.4 


IRAS 


04191+1523 


003 


1 


J042200. 


.41+153021. 


2 


IRAM04191+1522 


140 ± 10 


04 


22 


00 


.41 


+15 


30 


21. 


2 





.250 


20.8 


IRAS 


04191+1523 


004 


1 


J042838. 


.90+265135. 


6 


L1521F 


140 ± 10 


04 


28 


38 


.90 


+26 


51 


35. 


6 





.015 


380.9 


IRAS 


04260+2642 


005 


1 


J044112. 


.65+254635. 


.4 


TMC1 


140 ± 10 


04 


11 


12 


.65 


+25 


16 


35. 


.4 





.383 


15.7 


IRAS 


04381+2540 


006 


1 


J054429. 


26+090856. 


.8 


B35A 


400 ± 40 


05 


44 


29 


.26 


+09 


08 


56. 


.8 


0. 


141 


8.1 


IRAS 


05417+0907 


007 


6 


J 054443 


.94+090307. 


2 


B35A 


400 ± 40 


05 


44 


43 


.94 


+09 


03 


07. 


2 


0. 


006 


405.4 


IRAS 


05417+0907 


008 


5 


J071357 


.60—482711. 


9 


BHR22 


450 ± 50 


07 


13 


57 


.60 


-48 


27 


11. 


9 





.006 


170.7 


IRAS 


07128-4820 


009 


5 


J080533. 


.05-390924. 


.8 


BHR16 


440 ± 100 


08 


05 


33 


.05 


-39 


09 


24. 


.8 





.081 


241.9 


IRAS 


08036-3904 


010 


5 


J080850 


.47-360652. 


2 


CG30-31 


450 ± 50 


08 


08 


50 


.47 


-36 


06 


52. 


2 





.064 


525.5 


IRAS 


08076-3556 


011 


3 


J081705. 


.26-395415. 


8 


DC2573-25 


440 ± 100 


08 


17 


05 


.26 


-39 


54 


15. 


8 





193 


51.1 


IRAS 


08152-3945 


012 


5 


J081706. 


.34-394635. 


.8 


DC2573— 25 


440 ± 100 


08 


17 


06 


.34 


-39 


16 


35. 


8 


0. 


013 


194.3 


IRAS 


08151-3934 


013 


5 


J081710 


92—395244. 


4 


DC2573— 25 


440 ± 100 


08 


17 


10 


.92 


—39 


52 


14. 


1 





.018 


161.0 


IRAS 


08152-3945 


014 


6 


J092851 


50-513558. 


6 


DC2742— 04 


130 ± 10 


09 


28 


51 


.50 


-51 


35 


58. 


6 





.005 


712.4 


IRAS 


09274-5134 


015 


5 


J123040. 


90-710007. 


6 


Mu8 


150 ± 30 


12 


30 


40 


.90 


-71 


00 


07. 


6 





.001 


413.5 


IRAS 


12288-7047 


016 


2 


J124539. 


.96—552522. 


.4 


DC302. 1+7.4 


300 ± 50 


12 


45 


39 


.96 


-55 


25 


22. 


.4 





.232 


44.5 


IRAS 


12427-5508 


017 


1 


J130736 


89-770009. 


.7 


DC303.8— 14.2 


178 ± 18 


13 


07 


36 


.89 


-77 


00 


09. 


7 





.387 


5.7 


IRAS 


13036-7644 


018 


3 


J154216. 


99—524802. 


2 


DC3272+18 


250 ± 50 


15 


12 


16 


.99 


—52 


48 


02. 


2 





.034 


337.1 


IRAS 


15387-5233 


019 


5 


J155856. 


14-373354. 





DC3391+117— 2 


150 ± 20 


15 


58 


56 


.14 


-37 


33 


54. 








.005 


331.8 


IRAS 


15552-3727 


020 


5 


J163705 


11-353219. 


7 


DC3460+78 


150 ± 20 


16 


37 


05 


.11 


-35 


32 


19. 


7 





.003 


221.3 


IRAS 


16335-3528 


021 


1 


J165719. 


63-160923. 


.4 


CB68 


125 ± 25 


16 


57 


19 


.63 


-16 


09 


23. 


4 





.299 


2.0 


IRAS 


16544-1604 


022 


6 


J171111 


.83-272655. 





B59 


125 ± 25 


17 


11 


11 


.83 


-27 


26 


55. 








.024 


125.9 


IRAS 


17082-2724 


023 


3 


J171122 


.18-272602. 





B59 


125 ± 25 


17 


11 


22 


.18 


-27 


26 


02. 








.343 


86.3 


IRAS 


17081-2721 


024 


1 


J181616 


39-023237. 


7 


CB130-3 


270 ± 50 


18 


16 


16 


.39 


-02 


32 


37. 


.7 





.054 


905.7 


IRAS 


18144-0242 


025 


1 


J181659 


.47-180230. 


.5 


L328 


270 ± 50 


18 


16 


59 


.47 


-18 


02 


30. 


5 





.065 


445.0 


IRAS 


18145-1801 


026 


6 


J191744 


28+191523. 


.7 


L723 


300 ± 100 


19 


17 


14 


.28 


+19 


15 


23. 


7 





.048 


225.4 


IRAS 


19156+1906 


027 


1 


J192025 


32+112217. 


4 


L673 


300 ± 100 


19 


20 


25 


.32 


+11 


22 


17. 


4 


0. 


201 


42.9 


IRAS 


19180+1116 


028 


5 


J192025 


92+112221. 





L673 


300 ± 100 


19 


20 


25 


.92 


+ 11 


22 


21. 





0. 


168 


52.5 


IRAS 


19180+1116 


029 


1 


J192026. 


.16+111949. 


1 


L673 


300 ± 100 


19 


20 


26 


.16 


+11 


19 


49. 


1 


0. 


037 


5.3 


IRAS 


19180+1114 


030 


3 


J192026 


54+112025. 


4 


L673 


300 ± 100 


19 


20 


26 


.54 


+11 


20 


25. 


4 





138 


35.4 


IRAS 


19180+1114 


031 


1 


J192134 


.82+112123. 


4 


L673-7 


300 ± 100 


19 


21 


34 


.82 


+11 


21 


23. 


4 





.017 


427.6 


IRAS 


19189+1109 


032 


3 


J204056 


.66+672304. 


9 


L1148 


325 ± 25 


20 


10 


56 


.66 


+67 


23 


04. 


9 





.081 


222.3 


IRAS 


20410+6710 


033 


6 


J204105 


.95+671820. 


9 


L1148 


325 ± 25 


20 


41 


05 


.95 


+67 


18 


20. 


9 





.003 


250.9 


IRAS 


20410+6710 


034 


6 


J204355. 


.51+673850. 


3 


L1155E 


325 ± 25 


20 


43 


55 


.51 


+67 


38 


50. 


3 





.004 


642.1 


IRAS 


20423+6736 



Table 3 — Continued 



Source 




Spitzer 




c2d 


Dist. 




RA 




Dec 






L IR d 




Nearest IRAS 


Number 


Group a 


Source Name b 




Core Region 


(pc) 


(J2000) 


(J2000) 




(L©) 


(") 




Source 


035 


6 


J204427.17+673835. 


9 


L1155E 


325 ± 25 


20 


41 


27. 


17 


+67 


38 


35 


.9 


0.004 


768.6 


IRAS 


20423+6736 


036 


6 


J205706.72+773656. 


2 


L1228 


200 ± 50 


20 


57 


06. 


72 


+77 


36 


56 


.2 


0.199 


72.7 


IRAS 


20582+7724 


037 


6 


J205707.85+773659. 


.8 


L1228 


200 ± 50 


20 


57 


07. 


85 


+77 


36 


59 


.8 


0.038 


75.4 


IRAS 


20582+7724 


038 


1 


J212407.58+495908. 


9 


L1014 


250 ± 50 


21 


24 


07. 


.58 


+49 


59 


08. 


.9 


0.087 


246.8 


IRAS 


21220+4944 


039 


6 


J220633.22+590232. 


6 


L1165 


300 ± 50 


22 


06 


33. 


22 


+59 


02 


32 


.6 


0.030 


135.6 


IRAS 


22051+5848 


040 


6 


J220637.27+590315. 


.8 


L1165 


300 ± 50 


22 


06 


37. 


27 


+59 


03 


15. 


.8 


0.008 


107.5 


IRAS 


22051+5848 


041 


1 


J222807.42+690038. 


9 


L1221 


250 ± 50 


22 


28 


07. 


42 


+69 


00 


38. 


.9 


0.404 


41.9 


IRAS 


22266+6845 


042 


6 


J222933.39+751316. 





L1251 


300 ± 50 


22 


29 


33. 


39 


+75 


13 


16. 


.0 


0.066 


119.7 


IRAS 


22290+7458 


043 


3 


J222959.52+751403. 


1 


L1251 


300 ± 50 


22 


29 


59. 


52 


+75 


14 


03 


1 


0.272 


13.1 


IRAS 


22290+7458 


044 


1 


J223031. 94+751408. 


9 


L1251 


300 ± 50 


22 


30 


31. 


94 


+75 


14 


08. 


.9 


0.156 


120.5 


IRAS 


22290+7458 


045 


1 


J223105.59+751337. 


2 


L1251 


300 ± 50 


22 


31 


05. 


59 


+75 


13 


37 


.2 


0.077 


252.0 


IRAS 


22290+7458 


046 


6 


J223514.06+751502. 


5 


L1251 


300 ± 50 


22 


35 


14. 


06 


+75 


15 


02 


.5 


0.030 


129.3 


IRAS 


22343+7501 


047 


6 


J223731. 13+751041. 


5 


L1251 


300 ± 50 


22 


37 


31. 


13 


+75 


10 


41. 


.5 


0.004 


295.5 


IRAS 


22376+7455 


048 


1 


J223846.15+751132. 


3 


L1251 


300 ± 50 


22 


38 


46, 


15 


+75 


11 


32 


.3 


0.274 


5.2 


IRAS 


22376+7455 


049 


1 


J223846.44+751128. 





L1251 


300 ± 50 


22 


38 


46. 


44 


+75 


11 


28. 


.0 


0.081 


2.7 


IRAS 


22376+7455 



a Group indicating likelihood of being an embedded protostar (see EI5.2J I, 
b All source names are preceded by the prefix "SSTc2d ". 

c Name of the core region from Evans et al. (2007) in which this source is located. 

integrated luminosity using all available detections between 1.25 (2MASS J-band) and 70 fim. Entries in italics denote that the calculation only 
extends to 24 fim due to no flux information available at 70 fj,m. 
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Table 4. Candidate Low-Luminosity Embedded Objects in Perseus 



Source 
Number 


Group a 


Spitzer 
Source Name b 




RA 
(J2000) 


Dec 
(J2000) 




Lir c 

(Lq) 


(") 


Nearest IRAS 
Source 


050 


5 


J032353.74+310244. 


.5 


03 


23 


53. 


74 


+31 


02 


44 


.5 


0.010 


440.4 


IRAS 


03203+3048 


051 


5 


J032430.79+310651. 


.8 


03 


24 


30. 


79 


+31 


06 


51 


.8 


0.017 


419.3 


IRAS 


03219+3056 


052 


5 


J032433.14+311433. 





03 


24 


33. 


14 


+31 


14 


33 


.0 


0.008 


407.6 


IRAS 


03219+3103 


053 


5 


J032503.14+310322. 


.7 


03 


25 


03. 


14 


+31 


03 


22 


.7 


0.015 


178.7 


IRAS 


03222+3053 


054 


5 


J032515.26+312511. 


3 


03 


25 


15. 


26 


+31 


25 


11 


3 


0.008 


689.4 


IRAS 


03219+3103 


055 


1 


J032536.22+304515. 


8 


03 


25 


36. 


22 


+30 


45 


15 


.8 


0.061 


4.3 


IRAS 


03225+3034 


056 


1 


J032539. 12+304358. 


1 


03 


25 


39. 


12 


+30 


43 


58 


.1 


0.432 


90.4 


IRAS 


03225+3034 


057 


6 


J032621. 65+305726. 


3 


03 


26 


21. 


65 


+30 


57 


26 


.3 


0.006 


281.5 


IRAS 


03232+3051 


058 


5 


J032627.62+293002. 


2 


03 


26 


27. 


62 


+29 


30 


02 


.2 


0.016 


932.8 


IRAS 


03239+2905 


059 


5 


J032629.33+312154. 


4 


03 


26 


29. 


33 


+31 


21 


54 


.4 


0.018 


434.6 


IRAS 


03235+3118 


060 


3 


J032637.46+301528. 


1 


03 


26 


37. 


16 


+30 


15 


28 


.1 


0.462 


5.8 


IRAS 


03235+3004 


061 


5 


J032655.27+312831. 


4 


03 


26 


55. 


27 


+31 


28 


31 


.4 


0.008 


227.9 


IRAS 


03235+3118 


062 


5 


J032705.62+294601. 


9 


03 


27 


05. 


62 


+29 


46 


01 


.9 


0.008 


133.6 


IRAS 


03239+2933 


063 


1 


J032738.26+301358. 


.8 


03 


27 


38. 


26 


+30 


13 


58 


.8 


0.347 


60.0 


IRAS 


03245+3002 


064 


1 


J032832.57+311105. 


3 


03 


28 


32. 


57 


+31 


11 


05 


.3 


0.036 


163.6 


IRAS 


03255+3103 


065 


1 


J032839. 10+310601. 


8 


03 


28 


39. 


10 


+31 


06 


01 


.8 


0.015 


72.7 


IRAS 


03256+3055 


066 


6 


J032843.58+311736. 


2 


03 


28 


43. 


58 


+31 


17 


36 


.2 


0.4U 


15.8 


IRAS 


03256+3107 


067 


5 


J032843.70+291050. 


2 


03 


28 


43. 


70 


+29 


10 


50 


.2 


0.010 


757.9 


IRAS 


03253+2848 


068 


1 


J032845.29+310542. 





03 


28 


45. 


29 


+31 


05 


42 


.0 


0.185 


10.5 


IRAS 


03256+3055 


069 


6 


J032856.33+312227. 


.8 


03 


28 


56. 


33 


+31 


22 


27 


.8 


0.119 


182.5 


IRAS 


03260+3111 


070 


6 


J032856.59+310736. 


.8 


03 


28 


56. 


59 


+31 


07 


36 


.8 


0.014 


194.6 


IRAS 


03256+3055 


071 


1 


J032900.55+311200. 


7 


03 


29 


00. 


55 


+31 


12 


00 


.7 


0.105 


170.3 


IRAS 


03258+3104 


072 


6 


J032906.05+303039. 


2 


03 


29 


06. 


05 


+30 


30 


39 


.2 


0.019 


212.5 


IRAS 


03258+3023 


073 


1 


J032912.07+311301. 


6 


03 


29 


12 


07 


+31 


13 


01 


6 


0.018 


217.2 


IRAS 


03259+3105 


074 


2 


J032913.54+311358. 


1 


03 


29 


13. 


54 


+31 


13 


58 


.1 


0.305 


183.8 


IRAS 


03259+3105 


075 


2 


J032917.16+312746. 


4 


03 


29 


17. 


16 


+31 


27 


46 


.4 


0.183 


184.3 


IRAS 


03262+3114 


076 


6 


J032919.75+311256. 


9 


03 


29 


19. 


75 


+31 


12 


56 


.9 


0.002 


283.7 


IRAS 


03259+3105 


077 


6 


J032920.35+311250. 


.4 


03 


29 


20. 


35 


+31 


12 


50 


.4 


0.020 


293.9 


IRAS 


03259+3105 


078 


2 


J032923.47+313329. 


5 


03 


29 


23. 


47 


+31 


33 


29 


.5 


0.149 


9.3 


IRAS 


03262+3123 


079 


5 


J032950.42+300328. 


.8 


03 


29 


50. 


42 


+30 


03 


28 


.8 


0.010 


922.0 


IRAS 


03279+2950 


080 


3 


J032951. 82+313906. 


1 


03 


29 


51. 


82 


+31 


39 


06 


.1 


0.166 


3.3 


IRAS 


03267+3128 


081 


1 


J033032.69+302626. 


5 


03 


30 


32. 


69 


+30 


26 


26 


.5 


0.033 


147.0 


IRAS 


03273+3018 


082 


6 


J033047.81+313906. 


.8 


03 


30 


47. 


.81 


+31 


39 


06 


.8 


0.009 


402.2 


IRAS 


03279+3123 


083 


6 


J033116.90+302958. 


2 


03 


31 


16. 


90 


+30 


29 


58 


.2 


0.004 


383.7 


IRAS 


03286+3017 


084 


1 


J033120.98+304530. 


2 


03 


31 


20. 


98 


+30 


45 


30 


.2 


0.299 


9.6 


IRAS 


03282+3035 


085 


6 


J033146.63+302440. 





03 


31 


46. 


63 


+30 


24 


40 


.0 


0.015 


209.2 


IRAS 


03286+3017 


086 


6 


J033147.33+302438. 


5 


03 


31 


47. 


33 


+30 


24 


38 


.5 


0.002 


211.3 


IRAS 


03286+3017 


087 


5 


J033149.90+295423. 


8 


03 


31 


49. 


90 


+29 


54 


23 


.8 


0.014 


642.7 


IRAS 


03295+2947 


088 


1 


J033217.95+304947. 


6 


03 


32 


17. 


95 


+30 


49 


47 


.6 


0.135 


6.2 


IRAS 


03292+3039 


089 


5 


J033225.70+300631. 


3 


03 


32 


25. 


70 


+30 


06 


31 


.3 


0.007 


409.5 


IRAS 


03297+2951 


090 


1 


J033229.18+310240. 


9 


03 


32 


29. 


18 


+31 


02 


40 


.9 


0.180 


18.4 


IRAS 


03293+3052 


091 


6 


J033257.84+310608. 


3 


03 


32 


57. 


84 


+31 


06 


08 


.3 


0.012 


258.9 


IRAS 


03301+3057 


092 


1 


J033314.38+310710. 


9 


03 


33 


14. 


38 


+31 


07 


10 


.9 


0.098 


47.3 


IRAS 


03301+3057 


093 


1 


J033316.44+310652. 


6 


03 


33 


16. 


44 


+31 


06 


52 


6 


0.144 


58.4 


IRAS 


03301+3057 


094 


5 


J033357.17+314330. 





03 


33 


57. 


17 


+31 


43 


30 


.0 


0.050 


461.0 


IRAS 


03302+3131 
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Table 4 — Continued 



Source 




Spitzer 






RA 




Dec 






Lir c 




Nearest IRAS 


Number 


Group a 


Source Namc b 




(J2000) 


(J2000) 




(L©) 


(") 




Source 




K 
O 




r: 
o 


03 


35 


14. 


18 


+30 


23 


57. 


K 
D 





.057 




IRAS 


03324+3020 


uyo 


(-; 
O 


J uoo oo^t. <-± 1 -t-ou^t OU 1 . 


1 


03 


35 


54. 


41 


+30 


45 


01. 


1 
1 





.067 




IRAS 


03328+3035 


097 


6 


J033626.54+313634. 


9 


03 


36 


26. 


.54 


+31 


36 


31. 


9 





.026 


381.9 


IRAS 


03335+3120 


098 


6 


J033745.17+305445. 


.7 


03 


37 


45. 


.17 


+30 


51 


45. 


7 





.029 


694.5 


IRAS 


03337+3043 


099 


6 


J033812.17+314020. 


6 


03 


38 


12. 


17 


+31 


40 


20. 


6 





.019 


683.5 


IRAS 


03356+3121 


100 


6 


J033925.54+321708. 


2 


03 


39 


25. 


.54 


+32 


17 


08. 


2 





.292 


28.8 


IRAS 


03363+3207 


101 


5 


J034010.20+304446. 


.7 


03 


40 


10. 


20 


+30 


11 


46. 


7 





.008 


652.3 


IRAS 


03375+3043 


102 


5 


J034015.79+305504. 


8 


03 


10 


15. 


79 


+30 


55 


04. 


8 





.013 


326.0 


IRAS 


03375+3043 


103 


5 


J034227.67+310145. 


1 


03 


42 


27. 


67 


+31 


01 


45. 


1 





.019 


641.4 


IRAS 


03399+3059 


101 


3 


J034351. 02+320307. 


9 


03 


43 


51. 


02 


+32 


03 


07. 


9 





.172 


99.9 


IRAS 


03407+3152 


105 


1 


J034356.52+320052. 


9 


03 


43 


56. 


52 


+32 


00 


52. 


9 





.310 


59.0 


IRAS 


03407+3152 


106 


1 


J034356.83+320304. 


7 


03 


43 


56. 


83 


+32 


03 


04. 


7 





.253 


105.2 


IRAS 


03407+3152 


107 


3 


J034402.40+320204. 


.9 


03 


14 


02. 


.40 


+32 


02 


04. 


9 





.101 


124.0 


IRAS 


03407+3152 


108 


3 


J034402.64+320159. 


.5 


03 


11 


02. 


61 


+32 


01 


59. 


5 


0. 


021 


122.0 


IRAS 


03410+3152 


109 


3 


J034421. 36+315932. 


6 


03 


44 


21. 


36 


+31 


59 


32. 


6 





.200 


171.1 


IRAS 


03410+3152 


110 


6 


J034511. 57+322753. 


3 


03 


45 


11. 


.57 


+32 


27 


53. 


3 





.026 


512.3 


IRAS 


03426+3214 


111 


6 


J034513.82+321210. 


1 


03 


45 


13. 


82 


+32 


12 


10. 


1 





.051 


311.9 


IRAS 


03424+3203 


112 


6 


J034527.98+323635. 


6 


03 


45 


27. 


98 


+32 


36 


35. 


6 





Oil 


185.7 


IRAS 


03421+3229 


113 


5 


J034530.41+324801. 


.8 


03 


45 


30. 


.41 


+32 


18 


01. 


8 





.010 


247.3 


IRAS 


03424+3234 


114 


6 


J034632.45+313443. 


.7 


03 


46 


32. 


.45 


+31 


31 


43. 


7 





.018 


223.8 


IRAS 


03436+3123 


115 


5 


J034635.38+330746. 


6 


03 


46 


35. 


38 


+33 


07 


16. 


6 





.012 


198.8 


IRAS 


03436+3300 


116 


6 


J034656.74+325248. 


.4 


03 


46 


56. 


74 


+32 


52 


48. 


1 





.015 


549.0 


IRAS 


03434+3235 


117 


6 


J034734.34+324640. 


1 


03 


47 


31. 


.34 


+32 


46 


40. 


1 





.006 


314.5 


IRAS 


03445+3242 


118 


5 


J034754.26+331512. 


2 


03 


47 


51. 


26 


+33 


15 


12. 


2 





.013 


229.8 


IRAS 


03448+3302 


119 


6 


J034910.08+323249. 


2 


03 


49 


10. 


08 


+32 


32 


49. 


2 





.010 


594.5 


IRAS 


03454+3230 



a Group indicating likelihood of being an embedded protostar (see £15-21 )- 
b All source names are preceded by the prefix "SSTc2d ". 

integrated luminosity using all available detections between 1.25 (2MASS J-band) and 70 £im. Entries in italics 
denote that the calculation only extends to 24 fim due to no flux information available at 70 fim. 
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Table 5. Candidate Low-Luminosity Embedded Objects in Chamaeleon II 



Source 




Spitzer 






RA 


Dec 






L IR C 


D/iJAS 


Nearest IRAS 


Number 


Group 3 


Source Name b 




(J2000) 


(J2000) 




(L©) 


(") 




Source 




120 


5 


J124657.72- 


-771916. 


.0 


12 


16 


57.72 


-77 


19 


16 


.0 


0.004 


372.4 


IRAS 


12416- 


-7703 


121 


5 


J124816.42- 


-781511. 


.5 


12 


48 


16.42 


-78 


15 


11 


.5 


0.012 


896.9 


IRAS 


12440- 


-7813 


122 


5 


J124940.03- 


-775657. 


.1 


12 


49 


40.03 


-77 


56 


57. 


.1 


0.006 


897.8 


IRAS 


12504- 


-7745 


123 


5 


J125254.34- 


-780913. 


.7 


12 


52 


54.34 


-78 


09 


13. 


.7 


0.010 


483.2 


IRAS 


12504- 


-7745 


124 


3 


J125342.86- 


-771511. 


.5 


12 


53 


42.86 


-77 


15 


11. 


.5 


0.184 


5.7 


IRAS 


12500- 


-7658 


125 


5 


J125507.66- 


-765755. 


.4 


12 


55 


07.66 


-76 


57 


55. 


,4 


0.006 


182.1 


IRAS 


12522- 


-7640 


126 


5 


J125701.58- 


-764834. 


.9 


12 


57 


01.58 


-76 


48 


34. 


.9 


0.050 


10.5 


IRAS 


12533- 


-7632 


127 


5 


J130100.05- 


-770453. 


.0 


13 


01 


00.05 


-77 


04 


53 


.0 


0.003 


347.6 


IRAS 


12571- 


-7654 


128 


5 


J130242.22- 


-764402. 


.0 


13 


02 


42.22 


-76 


44 


02 


.0 


0.003 


396.9 


IRAS 


12584- 


-7621 


129 


5 


J130326.59- 


-764502. 


.9 


13 


03 


26.59 


-76 


45 


02 


.9 


0.020 


355.1 


IRAS 


13005- 


-7633 


130 


5 


J130509.53- 


-770716. 


.3 


13 


05 


09.53 


-77 


07 


16. 


.3 


0.005 


203.0 


IRAS 


13022- 


-7650 


131 


5 


J130515.96 


-771115. 


.7 


13 


05 


15.96 


-77 


11 


15. 


.7 


0.003 


336.8 


IRAS 


13022- 


-7650 


132 


1 


J130737.22- 


-770009. 


.4 


13 


07 


37.22 


-77 


00 


09. 


.4 


0.361 


6.2 


IRAS 


13036- 


-7644 


133 


5 


J130832.38- 


-773226. 


.2 


13 


08 


32.38 


-77 


32 


26 


.2 


0.003 


308.5 


IRAS 


13031- 


-7714 


134 


5 


J131031.92- 


-772423. 


.8 


13 


10 


31.92 


-77 


24 


23 


,8 


0.008 


694.4 


IRAS 


13030- 


-7707 


135 


5 


J131331.97- 


-773438. 


.3 


13 


13 


31.97 


-77 


34 


38. 


3 


0.018 


1271.6 


IRAS 


13031- 


-7714 



a Group indicating likelihood of being an embedded protostar (see i|5.2l l, 
b All source names are preceded by the prefix "SSTc2d ". 

integrated luminosity using all available detections between 1.25 (2MASS J-band) and 70 fim. Entries in italics 
denote that the calculation only extends to 24 [im due to no flux information available at 70 lira. 



Table 6. Candidate Low-Luminosity Embedded Objects in Lupus 



Source 


Spitzer 




RA 


Dec 


Lib. Diras 


Nearest IRAS 


Number Group a 


Source Name b 


Region 


(J2000) 


(J2000) 


(La) (") 


Source 



136 


5 


J153706 


36- 


334102 


4 


Lupus 




15 


37 


06 


36 


-33 


41 


02 


4 





004 


682.1 


IRAS 


15334- 


-3340 


137 


5 


J153751 


12- 


335806 


6 


Lupus 




15 


37 


51 


12 


-33 


58 


06 


6 





004 


680.6 


IRAS 


15337- 


-3350 


138 


5 


J153800 


05- 


340524 


1 


Lupus 




15 


38 


00 


05 


-34 


05 


24 


4 





004 


709.8 


IRAS 


15349- 


-3407 


139 


5 


J153826 


95- 


325628 


7 


Lupus 




15 


38 


26 


95 


-32 


56 


28 


7 





003 


178.5 


IRAS 


15350- 


-3246 


140 


5 


J153833 


84- 


334030 





Lupus 




15 


38 


33 


84 


-33 


40 


30 








004 


511.7 


IRAS 


15349- 


-3337 


1 11 


5 


J153834 


66- 


344819 


4 


Lupus 




15 


38 


34 


66 


-34 


48 


19 


4 





003 


486.1 


IRAS 


15358- 


-3444 


142 


5 


J153914 


59- 


341105 


3 


Lupus 




15 


39 


1 1 


59 


-34 


11 


05 


3 





003 


532.3 


IRAS 


15367- 


-3405 


143 


5 


J153953 


90- 


334821 


6 


Lupus 




15 


39 


53 


90 


-33 


48 


21 


6 





005 


186.7 


IRAS 


15364- 


-3339 


144 


5 


J154017 


64- 


324649 


i 


Lupus 




15 


10 


17 


61 


-32 


46 


49 


4 





006 


926.8 


IRAS 


15364- 


-3250 


145 


lc 


J154051 


62- 


342104 


7 


Lupus 




15 


10 


51 


62 


-34 


21 


04 


7 





011 


11.8 


IRAS 


15376- 


-3411 


146 


5 


J154125 


92- 


344445 


6 


Lupus 




15 


41 


25 


92 


-34 


41 


45 


6 





001 


659.8 


IRAS 


15379- 


-3445 


147 


5 


J154126 


95- 


325204 


1 


Lupus 




15 


41 


26 


95 


-32 


52 


04 


1 





004 


845.4 


IRAS 


15375- 


-3253 


148 


5 


J154138 


54- 


325005 


3 


Lupus 




15 


41 


38 


54 


-32 


50 


05 


3 





005 


733.6 


IRAS 


15387- 


-3228 


149 


5 


J154141 


93- 


324025 


7 


Lupus 




15 


41 


41 


93 


-32 


40 


25 


7 





008 


192.0 


IRAS 


15387- 


-3228 


150 


5 


J154148 


05- 


333536 


6 


Lupus 




15 


41 


48 


05 


-33 


35 


36 


6 





003 


203.4 


IRAS 


15384- 


-3327 


151 


5 


J154302 


28- 


344406 


4 


Lupus 




15 


43 


02 


28 


-34 


41 


06 


4 





007 


559.0 


IRAS 


15399- 


-3425 


152 


5 


J154339 


98- 


335554 


8 


Lupus 




15 


43 


39 


98 


-33 


55 


54 


8 





003 


603.3 


IRAS 


15397- 


-3350 


153 


5 


J154421 


91- 


332303 


8 


Lupus 




15 


44 


21 


91 


-33 


23 


03 


8 





008 


196.2 


IRAS 


15409- 


-3312 


154 


5 


J154507 


97- 


334254 


7 


Lupus 




15 


45 


07 


97 


-33 


42 


54 


7 





006 


375.2 


IRAS 


15417- 


-3327 


155 


5 


J154513 


15- 


335524 


2 


Lupus 




15 


45 


13 


15 


-33 


55 


24 


2 





006 


1110.8 


IRAS 


15417- 


-3327 


156 


5 


J154548 


24- 


340510 


7 


Lupus 




15 


45 


48 


24 


-34 


05 


10 


7 





006 


862.0 


IRAS 


15420- 


-3408 


157 


5 


J154601 


18- 


335816 


7 


Lupus 




15 


16 


01 


18 


-33 


58 


16 


7 





007 


1302.4 


IRAS 


15420- 


-3408 


158 


5 


J154612 


58- 


340250 


3 


Lupus 




15 


16 


12 


58 


-34 


02 


50 


3 





004 


1152.9 


IRAS 


15420- 


-3408 


159 


5 


J160003 


79- 


421024 


2 


Lupus 


IV 


16 


00 


03 


79 


-42 


10 


24 


2 





006 


487.1 


IRAS 


15570- 


-4155 


160 


5 


J160020 


16- 


410333 


8 


Lupus 


IV 


16 


00 


20 


16 


-41 


03 


33 


8 





003 


225.4 


IRAS 


15571- 


-4052 


161 


4a 


J160115 


55- 


415235 


4 


Lupus 


IV 


16 


01 


15 


55 


-41 


52 


35 


4 





049 


390.2 


IRAS 


15573- 


-4147 


162 


5 


J160132 


74- 


423628 


4 


Lupus 


IV 


16 


01 


32 


74 


-42 


36 


28 


4 





006 


637.5 


IRAS 


15571- 


-4230 


163 


5 


J160709 


79- 


392817 





Lupus 




16 


07 


09 


79 


-39 


28 


17 








008 


894.9 


IRAS 


16043- 


-3933 


164 


5 


J160727 


19- 


383842 





Lupus 




16 


07 


27 


19 


-38 


38 


12 








015 


960.9 


IRAS 


16051- 


-3820 


165 


5 


J160729 


66- 


391634 





Lupus 




16 


07 


29 


66 


-39 


16 


34 








009 


401.1 


IRAS 


16042- 


-3901 


166 


5 


J160740 


42- 


384602 


3 


Lupus 




16 


07 


40 


42 


-38 


46 


02 


3 





009 


1081.8 


IRAS 


16036- 


-3854 


167 


5 


J160741 


50- 


383614 





Lupus 




16 


07 


41 


50 


-38 


36 


14 








018 


737.8 


IRAS 


16051- 


-3820 


168 


5 


J160754 


74- 


391544 


6 


Lupus 




16 


07 


54 


71 


-39 


15 


44 


6 





012 


411.1 


IRAS 


16042- 


-3901 


169 


6 


J160851 


43- 


390530 


5 


Lupus 




16 


08 


51 


43 


-39 


05 


30 


5 





035 


52.1 


IRAS 


16054- 


-3857 


170 


1 


J160918 


07- 


390453 


4 


Lupus 




16 


09 


18 


07 


-39 


04 


53 


4 





134 


19.3 


IRAS 


16059- 


-3857 


171 


5 


J160934 


85- 


384654 


1 


Lupus 




16 


09 


34 


85 


-38 


16 


51 


1 





004 


715.3 


IRAS 


16059- 


-3850 


172 


5 


J160950 


09- 


380756 


3 


Lupus 




16 


09 


50 


09 


-38 


07 


56 


3 





014 


312.3 


IRAS 


16063- 


-3755 


173 


6 


J160954 


82- 


391122 


6 


Lupus 




16 


09 


54 


82 


-39 


11 


22 


6 





009 


63.2 


IRAS 


16064- 


-3903 


174 


5 


J161002 


52- 


383025 


9 


Lupus 




16 


10 


02 


52 


-38 


30 


25 


9 





021 


679.7 


IRAS 


16069- 


-3811 


175 


5 


J161022 


30- 


373618 


4 


Lupus 




16 


10 


22 


30 


-37 


36 


18 


4 





009 


611.7 


IRAS 


16072- 


-3738 


176 


5 


J161110 


30- 


372008 


9 


Lupus 




16 


11 


10 


30 


-37 


20 


08 


9 





008 


897.5 


IRAS 


16084- 


-3725 


177 


5 


J161316 


27- 


384827 


7 


Lupus 




16 


13 


16 


27 


-38 


48 


27 


7 





006 


531.3 


IRAS 


16106- 


-3842 



a Group indicating likelihood of being an embedded protostar (see !|5.2> . 

b All source names are preceded by the prefix "SSTc2d " . 

c Name of the Lupus region (I, III, or IV) in which this source is located. 
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integrated luminosity using all available detections between 1.25 (2MASS J-band) and 70 /im. Entries in italics denote that the 
calculation only extends to 24 (im due to no flux information available at 70 (im. 
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Table 7. Candidate Low-Luminosity Embedded Objects in Ophiuchus 



Source 




Spitzer 






RA 


Dec 






Lir c 




Nearest IRAS 


Number 


Group a 


Source Name b 




(J2000) 


(J2000) 




(Le) 


(") 




Source 




178 


5 


J162135.57- 


-224351. 


6 


16 


21 


35.57 


-22 


43 


51. 


6 


0.004 


1467.3 


IRAS 


16200- 


-2251 


179 


5 


J162145.12- 


-234231. 


.7 


16 


21 


45.12 


-23 


42 


31. 


.7 


0.033 


246.9 


IRAS 


16187- 


-2339 


180 


5 


J162458.99- 


-252119. 


1 


16 


24 


58.99 


-25 


21 


19. 


1 


0.012 


6.6 


IRAS 


16219- 


-2514 


181 


3 


J162648.48- 


-242838. 


.6 


16 


26 


48.48 


-24 


28 


38. 


.6 


0.103 


366.6 


IRAS 


16235- 


-2416 


182 


1 


J162705.23- 


-243629. 


.5 


16 


27 


05.23 


-24 


36 


29. 


5 


0.107 


55.0 


IRAS 


16240- 


-2430 


183 


3 


J162715.89- 


-243843. 


1 


16 


27 


15.89 


-24 


38 


43. 


1 


0.440 


163.5 


IRAS 


16240- 


-2430 


184 


3 


J162717.57- 


-242856. 


3 


16 


27 


17.57 


-24 


28 


56. 


3 


0.374 


10.4 


IRAS 


16242- 


-2422 


185 


1 


J162821.60- 


-243623. 


1 


16 


28 


21.60 


-24 


36 


23. 


,4 


0.046 


11.6 


IRAS 


16253- 


-2429 


186 


5 


J162843.68- 


-231918. 


1 


16 


28 


43.68 


-23 


19 


18. 


1 


0.000 


516.2 


IRAS 


16262- 


-2317 


187 


5 


J163122.87- 


-231655. 


2 


16 


31 


22.87 


-23 


16 


55. 


2 


0.003 


1354.8 


IRAS 


16275- 


-2251 


188 


6 


J163131.25- 


-242628. 





16 


31 


31.25 


-24 


26 


28. 





0.013 


70.9 


IRAS 


16285- 


-2421 


189 


6 


J163136.77- 


-240419. 


9 


16 


31 


36.77 


-24 


04 


19. 


9 


0.036 


51.1 


IRAS 


16285- 


-2358 


190 


5 


J163307.27- 


-225041. 


6 


16 


33 


07.27 


-22 


50 


41. 


6 


0.034 


15.0 


IRAS 


16301- 


-2244 


191 


5 


J163502.45- 


-235100. 


.7 


16 


35 


02.45 


-23 


51 


00. 


.7 


0.006 


1074.3 


IRAS 


16318- 


-2402 


192 


5 


J164001.80- 


-244429. 





16 


40 


01.80 


-24 


44 


29. 





0.004 


447.9 


IRAS 


16375- 


-2439 


193 


5 


J164114.35- 


-243758. 


1 


16 


41 


14.35 


-24 


37 


58. 


1 


0.007 


539.7 


IRAS 


16377- 


-2426 


194 


5 


J164147.64- 


-243030. 


2 


16 


41 


47.64 


-24 


30 


30. 


.2 


0.005 


591.5 


IRAS 


16381- 


-2419 


195 


5 


J164206.74- 


-233859. 


3 


16 


42 


06.74 


-23 


38 


59. 


3 


0.005 


538.5 


IRAS 


16384- 


-2334 


196 


5 


J164301.87- 


-245420. 


9 


16 


43 


01.87 


-24 


54 


20. 


9 


0.002 


295.6 


IRAS 


16403- 


-2447 


197 


5 


J164305.88- 


-242357. 


5 


16 


43 


05.88 


-24 


23 


57. 


5 


0.005 


345.0 


IRAS 


16396- 


-2419 


198 


5 


J164433.26- 


-233348. 


2 


16 


14 


33.26 


-23 


33 


48. 


2 


0.005 


329.6 


IRAS 


16411- 


-2329 


199 


5 


J164512.50- 


-233848. 


1 


16 


45 


12.50 


-23 


38 


48. 


1 


0.016 


375.3 


IRAS 


16420- 


-2327 



a Group indicating likelihood of being an embedded protostar (see i|5.2t . 
b All source names are preceded by the prefix "SSTc2d ". 

integrated luminosity using all available detections between 1.25 (2MASS J-band) and 70 (im. Entries in italics 
denote that the calculation only extends to 24 fim due to no flux information available at 70 fim. 
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Table 8. Candidate Low-Luminosity Embedded Objects in Serpens 



Source 




Spitzer 






RA 


Dec 










Nearest IRAS 


Number 


Group a 


Source Name b 




colhead(J2000) 


(J2000) 




(L©) 


(") 




Source 


200 


6 


J182841. 86-000321. 


2 


18 


28 


41.86 


-00 


03 


21 


.2 


0.153 


240.0 


IRAS 


18261-0009 


201 


6 


J182844.02+005338. 





18 


28 


44.02 


+00 


53 


38. 





0.146 


99.1 


IRAS 


18262+0050 


202 


3 


J182844.78+005125. 


9 


18 


28 


44.78 


+00 


51 


25. 


.9 


0.016 


40.9 


IRAS 


18262+0050 


203 


6 


J182844.95+005203. 


4 


18 


28 


44.95 


+00 


52 


03. 


1 


0.209 


14.2 


IRAS 


18262+0050 


204 


6 


J182854.91+001832. 


.8 


18 


28 


54.91 


+00 


18 


32. 


.8 


0.033 


321.8 


IRAS 


18267+0016 


205 


6 


J182902.11+003120. 


6 


18 


29 


02.11 


+00 


31 


20. 


.6 


0.032 


70.9 


IRAS 


18265+0028 


206 


6 


J182902.83+003009. 


7 


18 


29 


02.83 


+00 


30 


09. 


.7 


0.134 


52.0 


IRAS 


18265+0028 


207 


6 


J182902.95+003008. 


3 


18 


29 


02.95 


+00 


30 


08. 


.3 


0.035 


51.3 


IRAS 


18265+0028 


208 


1 


J182909.05+003127. 


.8 


18 


29 


09.05 


+00 


31 


27. 


.8 


0.009 


70.9 


IRAS 


18265+0028 


209 


6 


J182923.18+013454. 


.8 


18 


29 


23.18 


+01 


34 


54. 


.8 


0.040 


891.6 


IRAS 


18276+0124 


210 


2 


J182949.63+011522. 





18 


29 


49.63 


+01 


15 


22. 





0.395 


11.3 


IRAS 


18273+0113 


211 


3 


J182951. 98+011538. 


2 


18 


29 


51.98 


+01 


15 


38. 


2 


0.028 


31.2 


IRAS 


18273+0113 


212 


3 


J182952.51+003611. 


9 


18 


29 


52.51 


+00 


36 


11. 


.9 


0.198 


18.8 


IRAS 


18273+0034 


213 


2 


J182953.04+003606. 


.8 


18 


29 


53.04 


+00 


36 


06. 


.8 


0.169 


10.0 


IRAS 


18273+0034 


214 


2 


J182954.31+003601. 


.4 


18 


29 


54.31 


+00 


36 


01. 


.4 


0.119 


10.0 


IRAS 


18273+0034 


215 


2 


J182957.67+011304. 


.4 


18 


29 


57.67 


+01 


13 


04. 


.4 


0.166 


64.9 


IRAS 


18274+0112 


216 


6 


J183000.55+010304. 





18 


30 


00.55 


+01 


03 


04. 





0.025 


157.8 


IRAS 


18273+0059 


217 


6 


J183004.03+003424. 


2 


18 


30 


04.03 


+00 


34 


24. 


2 


0.151 


144.5 


IRAS 


18273+0030 


218 


6 


J183005.26+004104. 


6 


18 


30 


05.26 


+00 


41 


04. 


6 


0.056 


92.5 


IRAS 


18275+0040 



a Group indicating likelihood of being an embedded protostar (see ^5.21 . 
b All source names are preceded by the prefix "SSTc2d ". 

integrated luminosity using all available detections between 1.25 (2MASS J-band) and 70 ^m. Entries in italics denote 
that the calculation only extends to 24 fim due to no flux information available at 70 (im. 



Table 9. Groups Showing Likelihood Each Candidate is an Embedded Protostar 





Associated with 


Embedded in 


Associated with 








Region of High 


Region of High 


Region of High 


4.5 (im 


Number 


Group 


Volume Dcnsity a 


Volume Density a 


Column Density a 


Jets/Nebulosity a 


of Candidates' 3 


1 


Yes 


Yes 


N/A 


N/A 


10 


2 


Yes 


No/Unknown 


N/A 


Yes 


8 


3 


Yes 


No/Unknown 


N/A 


No/Unknown 


19 


4a 


Unknown 


No/Unknown 


Yes 


Yes 


1 


4b 


Unknown 


No/Unknown 


Yes 


No /Unknown 





4c 


Unknown 


No/Unknown 


No/Unknown 


Yes 


1 


5 


Unknown 


No/Unknown 


No/Unknown 


No /Unknown 


95 


6 


No 


N/A 


N/A 


N/A 


54 



a "Yes" indicates data exists and satisfies this condition. "No" indicates data exists and shows that this 
condition is not satisfied. "Unknown" indicates no data exists to evaluate this condition. As discussed in the 
text, "No" and "Unknown" are considered equivalent for all but "Associated with Region of High Volume 
Density" since negative results do not prove a source is not an embedded protostar. 

b Total number of the 218 candidates in this group. 
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Table 10. Candidates Divided Into Groups 



Group 


Source 
Number 


Core/Cloud 


(L©) 


MIPS3 b 


SIMBAD 


1 


001 


IRAM 04191+1522 


0.023 


N/A 


NAME IRAM 04191-IRS 


1 


002 


IRAM 04191+1522 


0. 033 


N/A 


IRAS F04191+1523 


1 


003 


IRAM 04191+1522 


0.250 


N/A 


IRAS F04191+1523 


1 


004 


L1521F 


0.015 


N/A 


NAME LDN 1521F-IRS 


1 


005 


TMC1 


0.383 


N/A 


IRAS F04381+2540 


1 


006 


B35A 


0.141 


N/A 


Barnard 35A 


1 


017 


DC303. 8-14.2 


0.387 


N/A 


BHR 86 


1 


021 


CB68 


0.299 


N/A 


[CB88] 68 


1 


024 


CB130-3 


0.054 


N/A 


[CB88] 130 


1 


025 


L328 


0.065 


N/A 


[VRC2001] L328 SMM 1 


1 


027 


L673 


0.201 


N/A 


[VRC2001] L673 SMM 1 


1 


029 


L673 


0.037 


N/A 


[VRC2001] L673 SMM 2 


1 


031 


L673-7 


0.017 


N/A 


NAME LDN 673 7 


1 


038 


L1014 


0.087 


N/A 


NAME LDN 1014-IRS 


1 


041 


L1221 


0.404 


N/A 


[LH2005] MM 2 


1 


044 


L1251 


0.156 


N/A 


LDN 1251A 


1 


045 


L1251 


0.077 


N/A 


[KP93] 3-6 


1 


048 


L1251 


0.274 


N/A 


LDN 1251B 


1 


049 


L1251 


0.081 


N/A 


LDN 1251B 


1 


055 


Perseus 


0.061 


C/S 


[BC86b] LDN 1448 IRS 3B 


1 


056 


Perseus 


0.432 


c/s 


NAME LDN 1448* 


1 


063 


Perseus 


0.347 


Y 


2MASS J03273825+3013585 


1 


064 


Perseus 


0.036 


Y 


[JJK2007] 9 


1 


065 


Perseus 


0.015 


Y 


[JJK2007] 12 


1 


068 


Perseus 


0.185 


Y 


IRAS 03256+3055 


1 


071 


Perseus 


0.105 


C/S 


[JCC87] IRAS 4B1 


1 


073 


Perseus 


0.018 


C/S 


[JCC87] IRAS 4B 


1 


081 


Perseus 


0.033 


Y 


SSTc2d J033032. 7+302626 


1 


084 


Perseus 


0.299 


Y 


[WBK2005] SMM J033135+30455 


1 


088 


Perseus 


0.135 


Y 


IRAS 03292+3039 


1 


090 


Perseus 


0.180 


Y 


IRAS 03293+3052 


1 


092 


Perseus 


0.098 


C/S 


[WBK2005] SMM J033328+31078 


1 


093 


Perseus 


0.144 


Y 


[JJK2007] 35 


1 


105 


Perseus 


0.310 


C/S 


HH 211 


1 


106 


Perseus 


0.253 


C/S 


NAME IC 348 MMS 


1 


132 c 


Chamaeleon II 


0.361 


Y 


BHR 86 


1 


170 


Lupus III 


0.134 


Y 


NAME LUPUS 3 MMS 


1 


182 


Ophiuchus 


0.107 


C/S 


CRBR 2403.7-2948 


1 


185 


Ophiuchus 


0.046 


Y 


[SSG2006] MMS126 


1 


208 


Serpens 


0.009 


C/S 


[DAB2006] VLA-7 


2 


016 


DC302. 1+7.4 


0.232 


N/A 


CG 19 


2 


074 


Perseus 


0.305 


C/S 


[JCC87] IRAS 4B2 


2 


075 


Perseus 


0.183 


C/S 


[JJK2007] 27 


2 


078 


Perseus 


0.149 


Y 


IRAS 03262+3123 


2 


210 


Serpens 


0.395 


C/S 


[SB86] S68 lb 



-61 - 



Table 10 — Continued 



Source LiR a 



Group 


Number 


Core/Cloud 




MIPS3 b 


SIMBAD 


2 


213 


Serpens 


0.169 


c/s 


IRAS 18273+0034 


2 


214 


Serpens 


0.119 


Y 


IRAS 18273+0034 


2 


215 


Serpens 


0.166 


C/S 


GCNM 104 


3 


Oil 


DC2573-25 


0.193 


N/A 


NAME DCLD 257.3-2.5 1 


3 


018 


DC3272+18 


0.034 


N/A 


[LM99] 135 


3 


023 


B59 


0.343 


N/A 


[BHB2007] 10 


3 


030 


L673 


0.138 


N/A 


[VRC2001] L673 SMM 2 


3 


032 


L1148 


0.081 


N/A 


IRAS F20404+6712 


3 


043 


L1251 


0.272 


N/A 


IRAS 22290+7458 


3 


060 


Perseus 


0.462 


Y 


[JJK2007] 4 


3 


080 


Perseus 


0.166 


Y 


IRAS 03267+3218 


3 


104 


Perseus 


0.172 


C/S 


[JJK2007] 42 


3 


107 


Perseus 


0.101 


C/S 


[JJK2007] 47 


3 


108 


Perseus 


0.021 


C/S 


CI* IC 348 LRL 54460 


3 


109 


Perseus 


0.200 


C/S 


CI* IC 348 LRL 1889 


3 


124 


Chamaeleon II 


0.184 


Y 


DENIS-P J125338. 9-771553 


3 


181 


Ophiuchus 


0.103 


C/S 


BBRCG 7 


3 


183 


Ophiuchus 


0.440 


C/S 


NAME WL 20W 


3 


184 


Ophiuchus 


0.374 


C/S 


BBRCG 36 


3 


202 


Serpens 


0.016 


N 


None 


3 


211 


Serpens 


0.028 


C/S 


NAME Serpens SMM 10 IR 


3 


212 


Serpens 


0.198 


C/S 


IRAS 18273+0034 


4a 


161 


Lupus IV 


0.049 


Y 


2MASS J16011549-4152351 


4c 


145 


Lupus I 


0.011 


Y 


IRAS 15376-3411 


5 


008 


BHR22 


0.006 


N/A 


None 


5 


009 


BHR16 


0.081 


N/A 


None 


5 


010 


CG30-31 


0.064 


N/A 


None 


5 


012 


DC2573-25 


0.013 


N/A 


None 


5 


013 


DC2573-25 


0.018 


N/A 


Dcld 257.3-02.5 


5 


015 


Mu8 


0.001 


N/A 


None 


5 


019 


DC3391+117-2 


0.005 


N/A 


None 


5 


020 


DC3460+78 


0.003 


N/A 


None 


5 


028 


L673 


0.168 


N/A 


None 


5 


050 


Perseus 


0.010 


N 


None 


5 


051 


Perseus 


0.017 


N 


None 


5 


052 


Perseus 


0.008 


N 


None 


5 


053 


Perseus 


0.015 


N 


None 


5 


054 


Perseus 


0.008 


N 


None 


5 


058 


Perseus 


0.016 


N 


None 


5 


059 


Perseus 


0.018 


N 


None 


5 


061 


Perseus 


0.008 


N 


None 


5 


062 


Perseus 


0.008 


N 


None 


5 


067 


Perseus 


0.010 


U 


None 


5 


079 


Perseus 


0.010 


N 


None 


5 


087 


Perseus 


0.014 


N 


None 



-62 - 



Table 10 — Continued 



Group 


Source 
Number 


Core/Cloud 


(L©) 


MIPS3 b 


SIMBAD 


5 


089 


Perseus 


0.007 


N 


None 


5 


094 


Perseus 


0.050 


N 


None 


5 


095 


Perseus 


0.057 


I J 


None 


5 


101 


Perseus 


0.008 


U 


None 


5 


102 


Perseus 


0.013 


N 


None 


5 


103 


Perseus 


0.019 


U 


None 


5 


113 


Perseus 


0.010 


N 


None 


5 


115 


Perseus 


0.012 


N 


None 


5 


118 


Perseus 


0.013 


N 


None 


5 


120 


Chamacleon II 


0.004 


U 


None 


5 


121 


Chamaeleon II 


0.012 


U 


None 


5 


122 


Chamaeleon II 


0.006 


u 


None 


5 


123 


Chamaeleon II 


0.010 


N 


None 


5 


125 


Chamaeleon II 


0.006 


N 


None 


5 


126 


Chamaeleon II 


0.050 


N 


IKAs 12533—7632 


5 


128 


Chamaeleon II 


0.003 


N 


None 


5 


129 


Chamaeleon II 


0.020 


N 


None 


5 


130 


Chamaeleon II 


0.005 


N 


None 


5 


131 


Chamacleon II 


0.003 


N 


None 


5 


133 


Chamaeleon II 


0.003 


N 


None 


5 


134 


Chamaeleon II 


0.008 


N 


None 


5 


135 


Chamacleon II 


0.018 


N 


None 


5 


136 


Lupus I 


0.004 


N 


None 


5 


137 


Lupus I 


0.004 


N 


None 


5 


138 


Lupus I 


0.004 


N 


None 


5 


139 


Lupus I 


0.003 


N 


None 


5 


140 


Lupus I 


0.004 


N 


None 


5 


141 


Lupus I 


0.003 


N 


None 


5 


142 


Lupus I 


0.003 


N 


None 


5 


143 


Lupus I 


0.005 


IN 


None 


5 


144 


Lupus I 


0.006 


u 


None 


5 


146 


Lupus I 


0.001 


N 


None 


5 


147 


Lupus I 


0.004 


N 


None 


5 


148 


Lupus I 


0.005 


U 


None 


5 


149 


Lupus I 


0.008 


U 


None 


5 


150 


Lupus I 


0.003 


N 


None 


5 


151 


Lupus I 


0.007 


N 


None 


5 


152 


Lupus I 


0.003 


N 


None 


5 


153 


Lupus I 


0.008 


U 


None 


5 


154 


Lupus I 


0.006 


U 


None 


5 


155 


Lupus I 


0.006 


u 


None 


5 


156 


Lupus I 


0.006 


N 


None 


5 


157 


Lupus I 


0.007 


N 


None 


5 


158 


Lupus I 


0.004 


N 


None 


5 


159 


Lupus IV 


0.006 


N 


None 
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Group 


Source 
Number 


Core/Cloud 


(L©) 


MIPS3 b 


SIMBAD 


5 


160 


Lupus IV 


0.003 


I J 


None 


5 


162 


Lupus IV 


0.006 


N 


None 


5 


163 


Lupus III 


0.008 


N 


None 


5 


164 


Lupus III 


0.015 


N 


None 


5 


165 


Lupus III 


0.009 


N 


None 


5 


166 


Lupus III 


0.009 


N 


None 


5 


167 


Lupus III 


0.018 


N 


None 


5 


168 


Lupus III 


0.012 


N 


ZIVLAbb JloU75475 — o91544o 


5 


171 


Lupus III 


0.004 


N 


None 


5 


172 


Lupus III 


0.014 


N 


None 


5 


174 


Lupus III 


0.021 


N 


None 


5 


175 


Lupus III 


0.009 


N 


None 


5 


176 


Lupus III 


0.008 


N 


None 


5 


177 


Lupus III 


0.006 


N 


None 


5 


178 


Ophiuchus 


0.004 


N 


None 


5 


179 


Ophiuchus 


0.033 


Y 


[L89J H7 


5 


180 


Ophiuchus 


0.012 


N 


lHAb Z16219— 2514 


5 


186 


Ophiuchus 


0.000 


N 


None 


5 


187 


Ophiuchus 


0.003 


N 


None 


5 


190 


Ophiuchus 


0.034 


N 


TTi AO 1 Pom O A A 

IKAs 16301 — 2244 


5 


191 


Ophiuchus 


0.006 


N 


None 


5 


192 


Ophiuchus 


0.004 


N 


None 


5 


193 


Ophiuchus 


0.007 


U 


None 


5 


194 


Ophiuchus 


0.005 


N 


None 


5 


195 


Ophiuchus 


0.005 


Y 


None 


5 


196 


Ophiuchus 


0.002 


N 


None 


5 


197 


Ophiuchus 


0.005 


N 


None 


5 


198 


Ophiuchus 


0.005 


N 


None 


5 


199 


Ophiuchus 


0.016 


Y 


None 


6 


007 


B35A 


0.006 


N/A 


None 


6 


014 


DC2742-04 


0.005 


N/A 


None 


6 


022 


B59 


0.024 


N/A 


2MASS J17111182-2726547 


6 


026 


L723 


0.048 


N/A 


None 


6 


033 


L1148 


0.003 


N/A 


None 


6 


034 


L1155E 


0.004 


N/A 


None 


6 


035 


L1155E 


0.004 


N/A 


None 


6 


036 


L1228 


0.199 


N/A 


2MASX J20570656+7736557 


6 


037 


L1228 


0.038 


N/A 


None 


6 


039 


L1165 


0.030 


N/A 


None 


6 


040 


L1165 


0.008 


N/A 


None 


6 


042 


L1251 


0.066 


N/A 


None 


6 


046 


L1251 


0.030 


N/A 


None 


6 


047 


L1251 


0.004 


N/A 


None 


6 


057 


Perseus 


0.006 


N 


None 


6 


066 


Perseus 


0.4U 


Y 


HBC 341 
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Group 


Source 
Number 


Core/Cloud 


(L©) 


MIPS3 b 


SIMBAD 


6 


069 


Perseus 


0.119 


c/s 


2MASS J03285630+3122279 


6 


070 


Perseus 


0.014 


N 


None 


6 


072 


Perseus 


0.019 


N 


None 


6 


076 


Perseus 


0.002 


C/S 


None 


6 


077 


Perseus 


0.020 


C/S 


HH5 


6 


082 


Perseus 


0.009 


N 


None 


6 


083 


Perseus 


0.004 


N 


None 


6 


085 


Perseus 


0.015 


N 


None 


6 


086 


Perseus 


0.002 


N 


None 


6 


091 


Perseus 


0.012 


N 


None 


6 


096 


Perseus 


0.067 


N 


IRAS 03328+3035 


6 


097 


Perseus 


0.026 


N 


None 


6 


098 


Perseus 


0.029 


N 


None 


6 


099 


Perseus 


0.019 


N 


None 


6 


100 


Perseus 


0.292 


Y 


None 


6 


110 


Perseus 


0.026 


N 


None 


6 


111 


Perseus 


0.051 


N 


CI* IC 348 LRL 904 


6 


112 


Perseus 


0.011 


N 


None 


6 


114 


Perseus 


0.018 


N 


None 


6 


116 


Perseus 


0.015 


N 


None 


6 


117 


Perseus 


0.006 


N 


None 


6 


119 


Perseus 


0.010 


N 


None 


6 


127 


Chamacleon II 


0.003 


N 


None 


g 


169 


T 11 T"\1 IC III 

Ijlipilb 111 


0.035 


C /S 


V* VI 1 9? Srn 
v v iiyi oco 


6 


173 


Lupus III 


0.009 


N 


None 


6 


188 


Ophiuchus 


0.013 


N 


None 


6 


189 


Ophiuchus 


0.036 


Y 


None 


6 


200 


Serpens 


0.153 


N 


None 


6 


201 


Serpens 


0.146 


N 


None 


6 


203 


Serpens 


0.209 


Y 


IRAS 18262+0050 


6 


204 


Serpens 


0.033 


Y 


None 


6 


205 


Serpens 


0.032 


C/S 


SSTc2d J192902.1+003120 


6 


206 


Serpens 


0.134 


C/S 


2MASS J18290283+0030092 


6 


207 


Serpens 


0.035 


C/S 


None 


6 


209 


Serpens 


0.040 


N 


None 


6 


216 


Serpens 


0.025 


C/S 


None 


6 


217 


Serpens 


0.151 


N 


None 


6 


218 


Serpens 


0.056 


N 


None 



integrated luminosity using all available detections between 1.25 and 70 fim. Italics 
denote that the calculation only extends to 24 fj,m due to no flux information available at 70 
(im. 

b Indicates whether each candidate is associated with a 160 (im source. Y: Yes; N: No; U: 
Unknown - located off the edge of the 160 fim map; C/S: Unknown - located in a region that 
is confused and/or saturated at 160 /im; N/A: Not applicable for the dense cores. 

c Sourcc 132 is the same as Source 017. 
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Table 11. SEDs of Groups 1-3 



Source 


Wavelength 


Flux Density 


<7 


Apcrture a 




Evolutionary 


Number 


(H 


(mjy) 


(mjy) 


\ ) 


b 

Reference 


Indicators 


001 


3.6 


0.22 


0.02 




1 


L IR + = 0.04 




4.5 


0.78 


0.06 




1 


L 70 , = 0.05 L 




5.8 


0.94 


0.08 




1 


L hn , = 0.12 ± 0.02 L P 

OOL v — ^ w 




8.0 


0.63 


0.06 




1 


J-/1-.S-.I / J-JGWivn O ~i~ 1 

OOL I a til i il ^ — -■- 




24 


16.4 


1.73 




1 


T hn , = 27 ± 3 K 

OOL 




60 


500 


100 


60 


2 






70 


875 


120 




1 






90 


800 


160 


60 


2 






160 


6000 


1200 


60 


2 






200 


10000 


2000 


60 


2 






350 


5000 


800 


40 


3 






450 


10000 


2000 


60 


2 






850 


2500 


500 


60 


2 






1300 


650 


65 


60 


2 




003 


1.25 


0.3 


0.03 




4 


L IR t = 0.42 L m 




1.65 


1.7 


0.2 




4 


L 70 . = 0.23 L« 




2.17 


8.3 


0.9 




4 


L bol = 0.50 ± 0.06 L 




3.6 


48.5 


3.42 




1 


LiV.r.l / /j^mm — 197 ZL 23 




4.5 


93.0 


6.60 




1 


Tk„, = 149 ±8K 




5.8 


135 


9.36 




1 






8.0 


179 


12.5 




1 






24 


985 


105 




1 






25 


1300 


117 




5 






60 


5910 


591 




5 






70 


4320 


591 




1 






100 


14400 


1730 




5 






450 


1400 


800 


40 


6 






850 


600 


100 


40 


6 






1300 


400 




60 


7 




004 


3.6 


0.42 


0.06 




8 


*W = °- 03 L © 




4.5 


0.24 


0.03 




1 






5.8 


0.35 


0.05 




1 


L 6o! = 0.13 ± 0.02 L Q 




8.0 


0.62 


0.06 




1 


Lbol/Lsmm — 7 lb 1 




24 


24.0 


2.52 




1 


T 6oi = 20 ± 3 K 




70 


509 


86.5 




1 






350 


8900 


1400 


40 


9 






450 


7000 


1800 


40 


10 






850 


1500 


500 


40 


6 






1200 


600 


150 


40 


11 




005 


1.25 


0.22 


0.05 




4 


= °- 64 L © 




1.65 


2.27 


0.17 




4 


iJSt = °- 36 L © 




2.17 


16.1 


0.88 




4 


L bol = 0.73 ± 0.07 Lq 




4.5 


139 


10.3 




1 


Lbol/Lsmm = 89 ± 9 




8.0 


346 


24.7 




1 


T bol = 165 ± 10 K 




12 


396 


32 




5 
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Table 11 — Continued 



Smirrp 


Wq VP 1 pTI 0"T Vl 


T^llly T~)pTT31TV 




Li LIL L 




Pivnl 1 1 1 i on a rv 


Wi im hpi* 

1 ^ LI 1 1 1 \_I L L 




I'm Tvl 




C"1 


R of or pti ro 

± l^lLl LULL 


Tti H i pa tor 1 ? 




24 


1720 


182 




1 






25 


2520 


151 




5 






60 


9810 


490 




5 






70 


6740 


919 




1 






100 


12600 


1260 




5 






450 


6410 


2800 


120 


12 






850 


1350 


110 


120 


12 






1200 


353 


71 


40 


11 




017 


3.6 


0.98 


0.07 




1 


L>« = 0.65 L Q 




4.5 


12.8 


0.89 




1 


L i°t = °- 46 L 




5.8 


40.1 


2.76 




1 


L bol = 0.91 ± 0.09 Lq 




8.0 


82.7 


5.75 




1 


L bot /L smm = 2502 ± 257 




24 


995 


105 




1 


T bol = 69 ± 5 K 




25 


1050 


83.9 




5 






60 


4100 


200 


43.5 


13 






70 


5450 


574 




1 






100 


8100 


1100 


43.5 


13 






200 


29000 


3400 


89.4 


13 






1200 


1260 


140 


120 


14 






1300 


225 


49.0 




15 




021 


3.6 


0.36 


0.03 




1 


L l"t = °- 50 L o 




4.5 


1.38 


0.12 




1 


L I°t = °- 54 l O 




5.8 


1.90 


0.15 




1 


L bol = 0.94 ± 0.16 L Q 




8.0 


1.12 


0.09 




1 


Lbol/ Lsmm — 16 i 3 




24 


659 


69.3 




1 


T bol = 52 ± 4 K 




25 


1480 


133 




5 






60 


19500 


1950 




5 






70 


13100 


1780 




1 






100 


33700 


8080 




5 






350 


19100 


2900 


40 


3 






450 


27800 


14900 


40 


6 






850 


9400 


2300 


120 


6 




024 


3.6 


1.44 


0.11 




1 


L t J « = 0.09 L Q 




4.5 


4.23 


0.30 




1 


LjO t = 0.07 L Q 




5.8 


5.85 


0.41 




1 


L bol = 0.20 ± 0.04 L Q 




8.0 


6.76 


0.47 




1 


Lbol / Lsmm — 9 i 2 




24 


54.1 


5.68 




1 


T bol = 55 ± 10 K 




70 


312 


35.0 




1 






350 


2800 


700 


40 


3 






450 


1100 


400 


40 


3 






850 


1480 


296 


120 


16 




025 


3.6 


1.55 


0.14 




1 


LlSt = 0-11 Lq 




4.5 


3.79 


0.30 




1 


Lj° t = 0.07 L Q 




5.8 


3.34 


0.28 




1 


L bol = 0.18 ± 0.04 L Q 




8.0 


1.75 


0.27 




1 


Lbol/ Lsmm — 9 i 2 
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Table 11 — Continued 



Source 
Number 


Wavelength 
H 


Flux Density 
(mjy) 


a 

(mjy) 


Aperture 3 
(") 


Reference 1 " 


Evolutionary 
Indicators 




24 
70 
350 


87.0 
337 
2000 


9.14 
109 
500 


20 


1 
1 

3 


T bol = 62 ± 9 K 


031 


4.5 
24 
70 
450 
1200 


0.11 
1.94 
148 
1100 

358 


0.01 
0.28 
23.8 
400 
72 


40 
80 


1 
1 
1 

3 

11 


Lift = 0-03 L 
Lj° t = 0.04 L 
L bol = 0.09 ± 0.03 L 

Lbol/Lsmm = 10 ± 3 

T bol = 24 ± 6 K 


038 


3.6 
4.5 
5.8 
8.0 
24 
70 
350 
450 
850 
1200 


4.02 
12.2 
21.2 
26.2 
89.5 
503 
3200 
21500 
1800 
630 


0.34 
0.91 
1.53 
1.85 
9.39 
70.4 
500 
16100 
400 
126 


40 
120 
120 
120 


1 
1 
1 
1 
1 
1 
3 

17 
17 
17 


L\St = 0-15 L 
Lj° t = 0.09 L 
L bol = 0.34 ± 0.11 L 

Lbol/Lsmm - 3 ± 1 

T bol = 66 ± 21 K 


041 


3.6 
4.5 
5.8 
8.0 
24 
70 
350 
450 
850 


0.57 

3.10 

4.96 

3.84 

47.5 

5080 

6400 

11300 

7900 


0.04 
0.24 
0.25 
0.35 
5.00 
691 
1000 
6100 
1900 


40 
40 
120 


1 
1 
1 
1 
1 
1 
3 
6 
6 


L\nt = 0-67 L 

L Tnt = °- 81 L 

L bol = 0.93 ± 0.16 L 
L bo i/L srnrn = 8 i 1 
T bol = 42 ± 4 K 


044 


3.6 
4.5 
5.8 
8.0 
24 
70 
350 
1200 


0.13 
0.42 
0.32 
0.16 
4.97 
1400 
9400 
1010 


0.02 
0.05 
0.04 
0.03 
0.54 
192 
1400 
202 


40 
80 


1 
1 
1 
1 
1 
1 
3 
11 


Ll"t = 0-26 L 
Lj°t = 0.34 L 
L bol = 0.78 ± 0.11 L 

Lbol/Lsmm = 8 ± 1 

T bol = 24 ± 3 K 


045 


3.6 
4.5 
5.8 
8.0 
24 
70 
350 
1200 


0.43 
1.07 
0.75 
0.37 
1.90 
688 
7500 
906 


0.04 
0.10 
0.06 
0.04 
0.23 
94.5 
1100 
181 


40 
80 


1 
1 
1 
1 
1 
1 
3 
11 


L« = 0.13 L 
Lj°t = 0.17 L 
L bol = 0.55 ± 0.08 L 

Lbol/Lsmm - 7 i 1 

T bol = 21 ± 2 K 


063 


1.25 
1.65 


0.24 
7.99 


0.05 
0.32 




4 
4 


LlZt = 0-58 L 
L\ nt = 0-20 L 
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\A/n Amnion o+h 


Hjiiy Tlr>Ti<3i't"'\r 

L 1UA J_y(_-ilOl Li^y 








i_j vui u iiuiiai y 


lMi i tn Vior 


(//m) 


Cm TV) 


I'm Tvl 


CM 




Tri H i na+or^ 




2.17 


36.4 


0.84 




4 


L bol = 1.1 ± 0.2 L Q 




3.6 


62.5 


4.43 




1 


Lbol/Lsmm — 36 i 7 




4.5 


72.0 


5.13 




1 


T bot = 199 ± 33 K 




5.8 


79.1 


5.46 




1 






8.0 


93.1 


6.57 




1 






24 


343 


36.2 




1 






70 


1120 


160 




1 






160 


16200 


4600 




1 






350 


9800 


1500 


40 


3 






450 


1600 




T 


18, 19 






850 


720 




T 


18, 19 






850 


260 




T 


20 






1100 


270 


30 


T 


21 




064 


1.25 


0.17 


0.01 




4 


L 'int = °- 06 L © 




1.65 


0.69 


0.03 




4 


Lj° t = 0.03 L© 




2.17 


1.07 


0.11 




4 


L bo i = 0.20 ± 0.05 L Q 




3.6 


1.01 


0.09 




1 


LboljLsmm — 12 lb 3 




4.5 


2.28 


0.17 




1 


T bol = 65 ± 12 K 




5.8 


2.57 


0.18 




1 






8.0 


4.01 


0.28 




1 






24 


55.3 


5.83 




1 






70 


179 


26.5 




1 






160 


3340 


996 




1 






350 


3100 


500 


40 


9 






850 


240 




T 


18, 19 






1100 


210 


70 


120 


21 




065 


3.6 


0.13 


0.01 




1 


L lSt = °- 03 L © 




4.5 


0.27 


0.02 




1 


Ljo = o,02 L 




5.8 


0.43 


0.04 




1 


L bol = 0.22 ± 0.06 Lq 




8.0 


0.94 


0.07 




1 


L ol/Lsmm = 15 ± 4 




24 


24.4 


2.58 




1 


T bol = 29 ± 3 K 




70 


94.0 


14.5 




1 






160 


4930 


1420 




1 






350 


1900 


300 


40 


9 






850 


750 




T 


18, 19 






850 


830 




T 


20 






1100 


290 


20 


T 


21 




068 


1.25 


0.06 


0.01 




4 


= 0.31 L 




1.65 


0.46 


0.02 




4 


Lj° t = 0.26 Lq 




2.17 


0.72 


0.10 




4 


L bol = 0.45 ± 0.09 L Q 




3.6 


1.22 


0.09 




1 


L bo i/L amm = 193 ± 38 




4.5 


2.73 


0.19 




1 


T bol = 64 ± 8 K 




5.8 


2.67 


0.19 




1 






8.0 


2.69 


0.19 




1 






12 


20 






22, 23 
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tJHJ Lll 


\A/si\rr>l r>vi ol n 


TlIiiy Tir i TiQi't"\r 

L ILIA J_y t-lltM Lj' 




1 l_JJi_-l L Lll l_> 




i_J VU1 LI L1L_>11<X1 y 


lMi i tyi Vipr 

i. > LI 1 1 1 1 




I'm Tvl 


I'm Tv"t 


(") 


pf PT'Pl 1 PP 
± V. 1 V. 1 L^ 1 1 L^ \_ 


TnHI pa+nr^ 

111 LI ILUliWl o 




24 


213 


22.6 




1 






25 


150 






22, 23 






60 


1430 


129 




5 






70 


1520 


214 




1 






100 


6000 






22, 23 






160 


4960 


1620 




1 






450 


552 


140 


T 


23 






850 


102 


12 


T 


23 




071 


4.5 


0.09 


0.01 




1 


L\% = 0.18 L 




5.8 


0.11 


0.03 




1 


Lj° t = 0.23 L 




24 


15.4 


1.63 




1 


L bol = 0.24 ± 0.09 L 




70 


1320 


183 




1 


Lbol/Lsmm = 24 ± 4 




450 


2100 




T 


18, 19 


T bol = 41 ± 4 K 




850 


720 




T 


18, 19 






850 


130 




T 


20 






1100 


370 


50 


80 


21 




081 


3.6 


0.06 


0.01 




1 


L\* = 0.06 L 




4.5 


0.20 


0.01 




1 


L i°t = °- 06 L 




5.8 


0.19 


0.04 




1 


L bol = 0.18 ± 0.04 L 




8.0 


0.16 


0.03 




1 


Lbol i Lsmm — 7 i 2 




24 


30.5 


3.21 




1 


T bol = 33 ± 4 K 




70 


310 


43.1 




1 






160 


2560 


865 




1 






350 


3200 


500 


40 


9 






1100 


800 


70 


120 


21 




084 


3.6 


0.13 


0.01 




1 


L\K = 0.50 L 




4.5 


0.85 


0.07 




1 


L?° t = 0.63 L 




5.8 


0.83 


0.07 




1 


L bol = 1.1 ± 0.2 L 




8.0 


0.62 


0.06 




1 


Lbol 1 ' Lsmm — 22 lb 4 




24 


14.7 


1.55 




1 


T bol = 35 ± 3 K 




60 


2332 


210 




5 






70 


3880 


533 




1 






100 


13900 


1520 




5 






160 


16400 


5680 




1 






350 


12300 


2100 


40 


3 






450 


3300 




T 


18, 19 






850 


2260 




T 


18, 19 






850 


910 




T 


20 






1100 


760 


30 


T 


21 




088 


3.6 


0.39 


0.08 




1 


Lint = 0- 2 3 L© 




4.5 


0.42 


0.14 




1 


L™t = °-29 L 




8.0 


1.29 


0.11 




1 


L bol = 1.2 ± 0.3 L 




24 


13.0 


1.38 




1 


Lbol/ Lsmm - 7 ± 2 




60 


1590 


127 




5 


T bol = 28 ± 2 K 




70 


1720 


183 




1 
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finnrrp 


Wa vol PT1 P"t" ll 


T^lnx l~)pTi<^i1~v 




A npi"tnrp 

lipLl Lr Ul l, 




Pivnl 1 1 1 i on a rv 


lMi i tyi hpr 






I'm Tv"t 


(") 

\ ) 


pf PT'Pl 1 PP 


Tti A \ ratnri 




100 


9087 


1270 




1 






160 


24200 


8380 




1 






350 


17500 


2800 


40 


3 






450 


21500 




T 


18, 19 






850 


5600 




T 


18, 19 






850 


2130 




T 


20 






1100 


2310 


60 


120 


21 




090 


1.25 


0.71 


0.06 




4 


L lnt = °- 30 L 




1.65 


2.69 


0.16 




4 


L J°t = °- 20 L 




2.17 


5.81 


0.21 




4 


L bol = 0.57 ± 0.11 L 




3.6 


11.1 


0.86 




1 


Lboll 'Lsmm = 30 i 6 




4.5 


17.2 


1.25 




1 


T bol = 114 ± 17 K 




5.8 


23.3 


1.62 




1 






8.0 


32.8 


2.28 




1 






24 


181 


19.0 




1 






60 


1350 


122 




5 






70 


1130 


157 




1 






160 


8250 


2380 




1 






350 


2900 


400 


40 


3 






1100 


290 


50 


80 


21 




092 


3.6 


0.16 


0.01 




1 


L\%. = 0.16 L 




4.5 


0.94 


0.07 




1 


L I°t = o- 14 L © 




5.8 


2.05 


0.15 




1 


L bol = 0.29 ± 0.04 L 




8.0 


4.21 


0.30 




1 


Lbol/ Lsmm = 8 ± 1 




24 


111 


11.8 




1 


T bol = 47 ± 5 K 




70 


809 


116 




1 






350 


3700 


600 


20 


9 






1100 


1840 


40 


T 


21 




093 


4.5 


0.26 


0.04 




1 


L 1H . = 0.24 L n 

%rtt <J 




5.8 


0.44 


0.05 




1 


L? n t = 0-31 L 




24 


19.5 


2.08 




1 


L bol = 1.1 ± 0.2 L 




70 


1820 


256 




1 


L bo l/ Lsmm = 4 ± 1 




160 


13600 


3930 




1 


T bol = 24 ± 2 K 




350 


13300 


500 


40 


9 






450 


50400 




T 


18, 19 






850 


5370 




T 


18, 19 






850 


4960 




T 


20 






1100 


1840 


40 


T 


21 




105 


3.6 


0.06 


0.02 




1 


Hnt = °- 52 L 




5.8 


0.28 


0.06 




1 


Lj°t = 0-66 L 




24 


2.20 


0.34 




1 


L bol = 1.8 ± 0.3 L 




70 


4080 


565 




1 


L bo ll Lsmm = 4 ± 1 




350 


25200 


4100 


40 


3 


T bol = 21 ± 3 K 




450 


83300 




T 


18, 19 






850 


8190 




T 


18, 19 
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Source 
Number 


Wavelength 

(H 


Flux Density 
(mjy) 


a 

(mjy) 


Aperture a 
(") 


Reference 15 


Evolutionary 
Indicators 




850 
1100 


3140 
2070 


50 


T 

80 


20 
21 




106 


3.6 
4.5 
5.8 
8.0 
24 
70 
450 
850 
850 
1100 


0.05 

0.20 

0.52 

0.60 

10.8 

3280 

17600 

3590 

2270 

1420 


0.01 
0.06 
0.06 
0.12 
1.17 
454 

30 


T 
T 
T 
T 


1 
1 
1 
1 
1 
1 

18, 19 
18, 19 

20 

21 


HZt = 0-42 L 
L™ = 0.54 L 
L bol = 1.1 ± 0.2 L 

Lbol/Lsmm — 15 i 3 

T bol = 25 ± 4 K 


170 


3.6 
4.5 
5.8 
8.0 
24 
60 
70 
160 
1200 


0.26 
1.00 
0.99 
0.55 
32.4 
1670 
2610 
8710 
990 


0.03 
0.08 
0.09 
0.06 
3.44 
217 
363 
2510 
130 


56 


1 
1 
1 
1 
1 
5 
1 
1 

24 


Hnt = °- 22 L © 
Lj° t = 0-29 L 
L bol = 0.37 ± 0.08 L 
Lhol/Lsmm = 275 ± 62 
T bol = 39 ± 4 K 


182 


1.25 
1.65 
2.17 
3.6 
4.5 
5.8 
8.0 
24 
70 
850 
1200 


0.06 
0.28 
1.11 
6.11 
10.7 
16.6 
29.6 
459 
3490 
440 
433 


0.01 
0.01 
0.09 
0.46 
0.80 
0.82 
1.17 
48.3 
526 


T 
T 


4 
4 
4 
1 
1 
1 
1 
1 
1 

25 
26 


L\« = 0.18 L 
L™ = 0.15 L 
L bol = 0.15 ± 0.02 L 
Lhol/Lsmm = 335 ± 46 
T bo , = 105 ±2K 


185 


3.6 

4.5 

5.8 

8.0 

24 

60 

70 

100 

160 

350 

1100 

1200 


0.59 
2.55 
4.49 
3.13 
90.1 
2911 
2000 
16300 
18700 
45400 
1140 
411 


0.05 
0.20 
0.37 
0.22 
9.47 
378 
281 
2780 
5560 
6800 
100 


40 
120 
T 


1 
1 
1 
1 
1 
5 
1 
5 
1 
9 
27 
26 


Hnt = °- 08 L © 

Lj° t = 0.09 L Q 
L bol = 0.45 ± 0.08 L 
L bo i/L srnrn = 7 i 1 
T bol = 30 ± 2 K 


016 


3.6 


0.13 


0.01 




1 


Hnt = 0-39 L 
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Source 
Number 


Wavelength 

(H 


Flux Density 
(mjy) 


a 

(mjy) 


Aperture a 
(") 


Reference 15 


Evolutionary 
Indicators 




4.5 
5.8 
8.0 
24 
70 
1200 


1.10 
2.04 
2.39 
87.7 
1750 
370 


0.08 
0.15 
0.16 
9.20 
240 
50 


120 


1 
1 
1 
1 
1 

14 


= °- 42 L © 

L bol = 0.37 ± 0.05 L Q 

Lbol/Lsmm = 321 ± 43 

T bol = 66 ± 1 K 


074 


3.6 
5.8 
24 
70 
450 
850 


0.03 
0.25 
28.4 
3910 
20900 
2920 


0.01 
0.04 
3.00 

538 


T 
T 


1 
1 
1 
1 

18, 19 
18, 19 


L 1H — 51 L^ 
L™t = 0-63 Lq 
L bo i = 1-3 ± 0.2 L 

Lbol/L srnrn = 16 i 3 

T bol = 26 ± 4 K 


075 


3.6 
4.5 
5.8 
8.0 
24 
70 
450 
850 
850 
1100 


0.42 
1.22 
1.24 
0.82 
111 
1950 
2500 
1420 
470 
550 


0.04 
0.12 
0.10 
0.11 
11.7 
292 

50 


T 
T 
T 

80 


1 
1 
1 
1 
1 
\ 

18, 19 
18, 19 

20 

21 


m = 0.31 L 
Lj° t = 0.33 L 
L bol = 0.37 ± 0.06 L 
Lbol/Lsmm — 25 i 4 
T bol = 53 ± 4 K 


078 


3.6 
4.5 
5.8 
8.0 
24 
60 
70 
160 
850 
850 
1100 


0.76 
2.26 
2.47 
3.03 
191 
1070 
1150 
4380 
1560 
390 
510 


0.08 
0.17 
0.24 
0.21 
20.1 
96.3 
159 
1270 

40 


T 
T 

80 


1 
1 
1 
1 
1 
5 

1 
1 

18, 19 
20 
21 


Hut = 0-25 L 
Lj° t = 0.20 L 
L bol = 0.38 ± 0.07 L 
Lbol/Lsmm = 134 ± 26 
T bol = 59 ± 7 K 


Oil 


1.25 
1.65 
2.17 
3.6 
4.5 
5.8 
8.0 
24 
70 
1200 


0.75 
1.87 
2.42 
2.17 
2.12 
1.81 
2.43 
56.4 
492 
100 


0.07 
0.14 
0.14 
0.15 
0.15 
0.13 
0.17 
5.92 
69.5 
20 


40 


4 
4 
4 
1 
1 
1 
1 
1 
1 

14 


L lnt = O- 32 L© 
Lj°t = 0-26 L 
Lboi = 0.28 ± 0.04 L 
Lbol/Lsmm = 416 ± 54 
Tboi = 206 ± 5 K 


018 


3.6 
4.5 


0.47 
2.03 


0.04 
0.16 




1 
1 


Hnt = 0- 06 L © 

Lj° t = 0.04 L 
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Source 


Wavelength 


Flux Density 


cr 


Aperture a 


Evolutionary 


Number 


Om) 


(mjy) 


(mjy) 


///\ t» r 1~> 

(") Reference" 5 


Indicators 




5.8 


4.16 


0.30 


1 


L bo i = 0.06 ± 0.01 L Q 




8.0 


4.94 


0.43 


1 


Li. i/L.™™ = 90 ± 12 




24 


42.9 


4.50 


1 


71 , — ins + a K 




70 


219 


32.9 


1 






1200 


290 


40 


80 14 




023 


3.6 


3.18 


0.23 


\ 


T 1H — n T 
n int — U - ' ^0 




4.5 


11.9 


0.84 


\ 






5.8 


22.1 


1.53 


\ 


7",, , > n 34 T,^ 




8.0 


42.2 


2.91 


\ 


T u , IT — 




24 


1870 


199 


\ 


J bol _^ lx 




70 
( u 


i 0900 








UOU 


O.O 


nil 

U. 11 


u.uo 


i 
± 


r I H — n 90 t „ 




4.5 


0.48 


0.05 


1 


7"/ 70 — 30 T,^-, 




5.8 


0.51 


0.07 




7",, . > 1 4 T,^ 




24 


7.61 


0.85 








70 


1230 


170 








1 25 


n oq 


n ni 

U.U1 


4 


j IR _ n 1 4 T ^ 




l.OO 


u.oy 


u.uo 


fi 


T 70 — n HQ T ^ 




9 17 


n &7 

u.o / 


u.uo 


4 


t , 7 — n 1 1 -t- n t ^ 

^bol — U.IO It U.UZ 1j0 




O.D 


1 91 


u.uy 


i 
i 


T 1 , 1 T — 907 A- 98 

^bol/^smm — I It Zo 




*±.o 


1 

i.oy 


n 1 n 
u. 1U 


i 
i 


7 1 , , — 1 A S -4- 1 Tv" 




O.O 


U.oo 


U.U / 


i 
i 






o.u 


1 25 


n no 
u.uy 








24 


70 ^ 

/U.O 


7 41 

/ .4:1 








25 


126 


15.1 


5, 28 






60 


321 


48.2 


5, 28 






70 


306 


41.8 








1200 


170 


35 


80 11 




043 


1.25 


0.33 


0.02 


4 






1.65 


1.76 


0.13 


4 


r70 — n 91 T,^ 




2.17 


6.43 


0.23 


4 


L bo i = 0.38 ± 0.04 L Q 




3.6 


15.8 


1.22 


1 


Lbol/Lsmm = 173 ± 19 




4.5 


21.8 


1.64 


1 


T 6o! = 236 ± 3 K 




5.8 


25.2 


1.82 


1 






8.0 


27.3 


1.93 


1 






24 


272 


28.5 


1 






25 


375 


48.8 


5 






60 


734 


58.7 


5 






70 


848 


117 


1 






1200 


705 


141 


80 11 




060 


1.25 


0.18 


0.01 


4 


Lj« = 0.77 L 




1.65 


0.50 


0.03 


4 


Lj° t = 0.69 L 




2.17 


1.19 


0.08 


4 


Lboi = 1-1 ± 0.2 L Q 




3.6 


3.62 


0.29 


1 


Lbol/Lsmm = 27 ± 4 




4.5 


9.27 


0.65 


1 


T 6oi = 64 ± 6 K 
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Rnnrrp 


Wq vp 1 pTi t?\ n 

V V <JjVI_ j H_.±1^Li1± 


± 1 LL-A. \_J \_ llol Li V 




A nprtnro 3. 




Pivnl 1 1 1 i on a rv 


TvJi i yti Vior 


(/zm) 


Cm TV) 


I'm Tvl 


(") 


R pf Pt'PTl PP 


Tn H i ra+or^ 




5.8 


11.1 


0.77 




1 






8.0 


12.0 


0.82 




1 






24 


396 


41.7 




1 






25 


545 


76.3 




5 






60 


4180 


376 




5 






70 


4300 


591 




1 






100 


6970 


697 




5 






160 


8970 


2550 




1 






450 


12100 




T 


18, 19 






850 


860 




T 


18, 19 






850 


170 




T 


20 






1100 


400 


60 


120 


21 




080 


4.5 


0.02 


0.005 




1 


L{* = 0.28 L 




5.8 


0.10 


0.03 




1 


L Jnt = °- 34 L © 




8.0 


0.37 


0.04 




1 


L bol = 0.64 ± 0.18 L Q 




24 


48.9 


5.15 




1 


Lhol/Lsmm — 11 i 3 




60 


1880 


151 




5 


T hol = 34 ± 7 K 




70 


1990 


276 




1 






160 


5960 


1730 




1 






450 


11900 


6400 


40 


6 






850 


4900 


1200 


120 


6 






850 


410 




T 


20 






1100 


520 


40 


80 


21 




104 


3.6 


6.81 


0.59 




1 


L\* = 0.29 L 




4.5 


3.99 


0.63 




1 


L™ t = 0.33 L Q 




5.8 


12.8 


0.93 




1 


L bol = 0.32 ± 0.05 L© 




8.0 


2.15 


0.18 




1 


Lhol/Lsmm = 56 ± 8 




24 


45.4 


4.82 




1 


T bol = 63 ± 3 K 




70 


1960 


276 




1 






850 


1460 




T 


20 






1100 


1410 


40 


80 


21 




107 


3.6 


0.73 


0.05 




1 


HZt = 0-17 L 




4.5 


1.93 


0.15 




1 


L?nt = 0-15 L Q 




5.8 


3.82 


0.27 




1 


L bol = 0.18 ± 0.03 L 




8.0 


6.41 


0.45 




1 


Lbol/Lsmm = 77 ± 11 




24 


107 


11.7 




1 


T bol = 76 ± 4 K 




70 


847 


124 




1 






850 


940 




T 


20 






1100 


500 


20 


T 


21 




109 


1.25 


0.33 


0.02 




4 


Lint = 0-33 L 




1.65 


1.17 


0.09 




4 


Lj° t = O.H L Q 




2.17 


6.63 


0.18 




4 


L bol = 0.26 ± 0.03 L 




3.6 


31.0 


2.14 




1 


Lbol/Lsmm = 195 ± 20 




4.5 


50.4 


3.45 




1 


T bol = 345 ±4K 




5.8 


61.7 


4.21 




1 
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Smirrp 


W?) Vpl PT1 P"t" ll 


T^liix T~)pvi<>i1~v 




A nprtnrp a 




Th 1 . vol i ii" i on pi rv 


INJi i \Y] \)nr 
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4 
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4 
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Source 


Wavclc n gt h 


Flux Density 


(7 


Apcrturc a 




Evolutionary 


Number 


(yum) 


(mjy) 


(mjy) 


( ) 


Reference 13 


Indicators 




25 


6560 


1120 




5 






70 


12200 


1720 




1 






100 


22000 


4840 




5 






1200 


375 




T 


26 




161 


1.25 


0.42 


0.06 




4 


Hnt = °- 08 L © 




1.65 


1.67 


0.10 




4 


Lj° t = 0.08 L 




2.17 


4.83 


0.20 




4 


L bot > 0.11 L 




3.6 


8.36 


0.61 




1 


^bol/ Lsmm — ••• 




4.5 


9.92 


0.70 




1 


Uoi < 126 K 




5.8 


8.98 


0.63 




1 






8.0 


7.70 


0.53 




1 






24 


75.9 


7.99 




1 






70 


1220 
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1 






160 


5390 


1560 




1 




145 


5.8 


0.28 


0.09 




1 


Hnt = °- 02 L 




8.0 


0.70 


0.19 




1 


Lj°t = 0-03 L 




24 


8.78 


1.00 




1 


L bol > 0.02 L 




60 


468 


42.1 




5 


^boll Lsmm — ••■ 




70 


372 


54.9 




1 


T bol < 57 K 




160 


705 


225 




1 





a Diameter of the aperture used for photometry at A > 100 ^tm, where the sources are likely to be 
extended, when specified. A "T" indicates the photometry is an attempt to derive the total source 
flux density rather than the flux density in a fixed aperture; see individual references for descriptions 
of the methods used to derive this quantity. 

k References . — (1) Spitzer observations from this work; (2) Andre et al. (1999); (3) Wu et al. 
(2007); (4) 2MASS; (5) IRAS; (6) Young et al. (2006a); (7) Motte & Andre (2001); (8) Bourke et 
al. (2006); (9) M. Dunham et al. (2008), in preparation; (10) Kirk ct al. (2005); (11) Kauffmann 
et al. (2008); (12) Young ct al. (2003); (13) Lehtinen et al. (2005); (14) K. Brede et al. (2008), in 
preparation; (15) Hcnning ct al. (1993); (16) H.J. Kim et al. (2008, in preparation); (17) Young et al. 
(2004); (18) Hatchcll ct al. (2005); (19) Hatchell et al. (2007a); (20) Kirk ct al. (2006); (21) Enoch ct 
al. (2006); (22) Clark (1991); (23) Hodapp et al. (2005); (24) Tachihara et al. (2007); (25) Johnstone 
ct al. (2000); (26) Stankc ct al. (2006); (27) Young ct al. (2006b); (28) Kun (1998); (29) Andre & 
Montmerle (1994). 
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Fig. 1. — (a): Log L int vs. Log for the objects listed in Table [TJ The solid line shows the 
results of a linear least-squares fit in log-log space; it has a slope of 0.88 and a y-intercept 
of 0.32. (b): Li n t/ Li r vs. Log Lm for the same objects. The solid line shows the average 
ratio of 1.7, weighted by the uncertainties in Li nt . 
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Fig. 2. — Histogram of the un-normalized "probability" of being a galaxy, P ga i, using the 
method of Harvey et al. (2007b) for the 851 sources assigned such a probability in the 
ensemble of 82 regions with dense cores observed by c2d. The dotted line shows Log(P ga i) = 
— 1.47 while the dashed line shows Log(P ga i) = —1.25. 
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Fig. 3. — (normalized to 140 pc) vs. L int for the 11 embedded protostars listed in Table 
CD From left-to-right, top-to-bottom: IRAC1 (3.6 //m), IRAC2 (4.5 /mi), IRAC3 (5.8 fim), 
IRAC4 (8.0 /im), MIPS1 (24 /iin), and MIPS2 (70 /mri). The lines represent the results of 
linear least-squares fits at each wavelength. The parameters of the fits are given in Table [2J 
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Fig. 4. — Same as Figure [3], except for the 1460 model SEDs calculated from the grid of 2-D 
radiative transfer models rather than the 11 embedded protostars listed in Table [U 
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Fig. 5. — Histogram of the 70 /im 3a point source sensitivity for the 82 regions with dense 
cores observed by c2d. The distribution has a mean and median of 38.6 mJy and 31.9 mJy, 
respectively. 
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Fig. 6. — Log of the luminosity sensitivity limit for embedded protostars in the 82 regions 
with dense cores observed by c2d (L) vs. the distance to each core (d), calculated by 
translating the 70 /im 3a point source sensitivity for each core into a luminosity sensitivity 
using the correlation found between F 70 and L int and then scaling from 140 pc to the distance 
to the core. The solid line shows the relation L int = 4 x 10~ 3 (rf/140pc) 2 L . 
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Fig. 7. — SED of SSTc2d J032856. 64+311835. 6, a source in Perseus that is representative of 
candidates that pass all 7 selection criteria described in §3.11 but feature SEDs inconsistent 
with being embedded sources. This particular source is SSS 108, a previously known pre- 
main sequence object with a 2 /zm excess indicative of a circumstellar disk. Objects with 
these types of SEDs are rejected (see §3.41 for further information). 
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Fig. 8. — Average SEDs, weighted by l/L^ou for the embedded protostars identified by 
Enoch (2007) and Enoch et al. (2008, in preparation) in Perseus, Serpens, and Ophiuchus 
based on a comparison between 1.1 mm Bolocam dust continuum emission maps and the 
Spitzer c2d maps. The average Class SED is shown in black, the average SED for Class I 
sources with rising (or flat) fluxes from 24 to 70 fim is shown in blue, the average SED for 
Class I sources with decreasing fluxes from 24 to 70 /zm is shown in green, and the average 
SED for the Class II objects identified by Enoch et al. is shown in red. 
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Fig. 9. — Distribution of luminosities for the 851 objects in the 82 regions with dense cores 
observed by c2d classified as either candidate YSOs or candidate galaxies by the classification 
method of Harvey et al. (2007b). The solid line shows the distribution for the sources 
classified as candidate YSOs while the dashed line shows the same for the sources classified 
as candidate galaxies. Of the 604 out of 851 sources with Ljr < 0.05 L , 518 (~ 86%) are 
classified as candidate galaxies. 
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Fig. 10. — Top Left: Three-color image of L673-7 comprised of IRAC 1 (3.6 /mi; blue), 
IRAC 2 (4.5 /mi; green), and IRAC 4 (8.0 /mi; red). The white cross marks the position of 
Source 031 from this work, which is not detected by IRAC; it is clearly seen to fall within 
a dark core at 8.0 /mi. Top Right: Same as top left, except the red now shows the MIPS1 
(24 jum) image. A very red source is seen at the position of Source 031. Bottom: Radial 
profiles of the background intensity at 8 and 24 /im, centered at the position of Source 031. 
The background intensities are calculated in concentric annuli, each with a radius of 2". The 
radius on the x-axis is the mid-point of the annulus (for example, at a radius of 10" the 
background intensity was calculated in an annulus with its inner edge located 9" from the 
source and its outer edge located 11" from the source). The calculations start at radii of 5 
and 7" respectively, for 8.0 and 24 jum. 
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Fig. 11. — L(nt, the internal luminosity estimated from the calculated value of L IR} vs. Lj® t , 
the internal luminosity estimated from the 70 /im flux scaled to its value at 140 pc, for the 
50 objects listed in Table [TT1 The dashed line shows the 1:1 line, while the dotted lines show 
the 2:1 and 1:2 lines. 
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Fig. 12. — Distribution of internal luminosities of the 50 embedded, low-luminosity protostars 
listed in Table (TTJ The internal luminosity of each protostar was estimated from the observed 
70 /j,m flux and the correlation between the two derived in §|3j The dashed lines shows 
our sensitivity limit to embedded protostars for the closest regions observed by c2d (L int 
= 2.5 x 10~ 3 L ) and the most distant regions observed by c2d (L int = 5 x 10~ 2 L ). 
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Fig. 13. — Distribution of internal luminosities, in linear bins with a size of 0.1 L Q , for all 
sources listed in Table [TT1 with L int < 0.5 L . 
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Fig. 14. — Figure 19 from Young & Evans (2005), a Lu vs. T^oi BLT, with the 50 objects 
listed in Table [11] overlaid as large, colored points. The meaning of the colors and symbols 
are described in the text. The uncertainties in and T^i, which are not shown on this 
plot, are dominated by the 20-60% errors introduced by incomplete, finite sampling of the 
source SEDs (see Appendix 2). The thick lines show the evolutionary tracks for the three 
models considered by Young & Evans (2005), which differ in their initial envelope mass. The 
thin lines show the evolutionary tracks for three models considered by Myers et al. (1998), 
which also differ in their initial envelope mass. The small circles are observations from the 
literature compiled by Young & Evans (2005). The vertical dashed lines show the Class 0/1 
and Class I /II Tboi boundaries from Chen et al. (1995). 
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Fig. 15. — Lb i/L srnm vs. Tboi for the 50 objects listed in Table dU The meaning of the 
colors / symbols are the same as for Figure [TH The Class divisions in Tb i are from Chen et 
al. (1995), while the Class divisions in L bo i/L smm are from Young & Evans (2005). 
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Fig. 16. — Same as Figure [HI except only showing the 15 VeLLOs (objects with L int < 0.1 
L ). 
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Fig. 17. — Distribution of the sensitivity to protostars embedded within the 300 dense cores 
in Perseus, Serpens, Ophiuchus, and the sample of 82 regions with dense cores targeted for 
observations. The dashed lines show, from left to right, the sensitivity for cores in Lupus I 
and IV, Chamaeleon II, and Lupus III. 
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Fig. 18. — Monte Carlo simulation of the observation of 10,000 embedded protostars with 
internal luminosities distributed evenly between 10~ 5 and lO^ 1 L and randomly placed in 
one of the 300 dense cores. Top: Normalized distribution of L int for the 10,000 sources. 
Bottom: Results of the simulation. The solid line shows the observed distribution of L int 
(normalized) for the 15 VeLLOs identified by this work, while the dashed line shows the 
simulated observed distribution of L int (normalized) of the simulated sources. 
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Fig. 19. — Distribution of internal luminosities, in linear bins of 0.01 L , for the 15 embedded 
protostars with L int < 0.1 L Q . The vertical dashed line shows L int = 10 -18 L 0.016 L . 
The vertical dotted line shows L int = 0.0025 L . 
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Fig. 20. — Same as Figure [IH except with the internal luminosities of 5,000 of the 10,000 
sources evenly distributed between 10~ 4 and 10~ 3 L Q and the internal luminosities of the 
other 5,000 sources evenly distributed between 10~ 3 and 10 _1 L©. 
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Fig. 21. — Percent errors between L& / calculated from the category 2 (left; category 3, 
middle; category 4, right) SED and the category 1 SED, as defined in the text. The dashed 
line shows the average value of this error for each category. 
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Same as Figure I2TI except for T bol rather than L hoi . 



